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ABSTRACT 


This paper presents a method of evaluating the 
performance of a solid propellant rocket motor of 
fixed geometry and given propellant characteristics 
using the constant (ratio of propellant burning- 

surface area to nozzle-throat area) process and 
group transformation method. This method does not 
require a prior knowledge of the K^. values, but re- 
quires only that the parameters to be evaluated be 
selected at constant regressed distances normal to 
the original propellant surface. Two computer pro- 
grams utilizing this method are presented. The first 
performs the evaluation for general performance pa- 
rameters over selected time intervals, and the second 
performs the evaluation for instantaneous performance 
versus time. 
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SOLID PROPELLANT ROCKET MOTOR PERFORMANCE COMPUTER 

PROGRAMS USING THE GROUP TRANSFORMATION METHOD 

By B. J. Lee and P. B. Burchfield 
Manned Spacecraft Center 

SUMMARY 


A solid propellant rocket motor of fixed geometry and given propellant 
will yield different performance with various prefire propellant temperatures. 
This paper presents two computer programs that transform the test data from 
as many as three different prefire propellant temperature groups to a like 
number of common prefire propellant temperatures of interest. Statistical 
calculations are made on the test data after they have been transformed. A 
detailed description of these computer programs is presented. 


INTRODUCTION 


The internal ballistic performance of a solid propellant rocket motor is 
a function of geometry, propellant temperature, and propellant composition. 
Therefore, a solid propellant rocket motor of fixed geometry (figs. 1(a) 
through 1(c)) and given propellant characteristics (fig. 2) will yield different 
performance with various propellant temperatures (fig. 3 and refs. 1 and 2). 
This phenomenon is due to the variation of propellant burning rate with propel- 
lant temperature (fig. 2) and will be referred to as the thermal sensitivity of 
the propellant. 

By using adequate experimental test data, the variation of motor per- 
formance with propellant temperature can be mathematically defined. With 
proper utilization of this mathematical definition the performance test data 
(from motor test firings conducted at different propellant temperatures) can 
be transformed to any common temperature of interest, within the extremities 
of the test data. 

There is generally a limited number of motors allotted to determine the 
performance variations (performance versus time and product variance) of a 



given motor design. The transformation of test data from different prefire 
propellant temperatures to given temperatures of interest is utilized in this 
report to obtain greater confidence in the statistical analyses performed on 
the test data. 

This report describes two solid propellant rocket motor performance 
computer programs that were written at the Manned Spacecraft Center. These 
programs transform experimental test data from as many as three different 
prefire propellant temperature groups to a like number of propellant tempera- 
tures of interest. The computer program presented in appendix A performs 
the transformations and calculates general performance parameters over 
selected time intervals. The computer program presented in appendix B per- 
forms the transformation for instantaneous performance versus time. 

Statistical analyses are performed on the test data after they have been 
transformed. The statistical analyses consist of calculating means, standard 
deviations, one-sided tolerance limits, and two-sided tolerance limits. 


SYMBOLS 


A area 

a burning-rate equation coefficient, r = aP n 

b KfJ 

e base of natural logarithm, 2. 71828. . . 

F longitudinal thrust 

K tolerance factor 

Kjj ratio of propellant burning-surface area to nozzle-throat area 
In natural logarithm 

m K n equation exponent, K N = bP m 

N sample size 


equation coefficient, K^. = bP 
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n burning-rate equation exponent, r = aP n 
P chamber pressure ; percentage of population 

P a ambient pressure 

r burning rate 

s estimated standard deviation 

T prefire propellant temperature 

t motor operating time interval 

u tolerance limit 

w web thickness 

X experimentally determined values that make up a sample 
y probability or confidence level 

7Tj£ temperature sensitivity coefficient of chamber pressure 
CTj ^ temperature sensitivity coefficient of burning rate 

0 nozzle cant angle 

Subscripts : 
d desired 

e exit 

f value of the parameter that corresponds to as determined from a 
second-order least-squares curve fit of the parameter versus T 

g prefire propellant temperature group whose data are to be transformed 
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h value of the parameter that the explicit motor would be expected to 
experience if the prefire propellant temperature had been 

i value of the parameter experimentally acquired from the explicit motor 
whose data are to be transformed 

j value of the parameter that corresponds to T as determined from 

gp 

a second-order least-squares curve fit of the parameter versus T 
1 lower 

t throat 

u upper 

v vacuum 

1 one-sided 

2 two-sided 

Operator : 

— average (such as F) 

RELATED SOLID PROPELLANT ROCKET MOTOR THEORY 


Ratio of Propellant Burning- Surface Area to Nozzle- Throat Area 

Solid propellants burn in parallel layers and regress normal to the pro- 
pellant surface (fig. 1(c)). Therefore, at any increasing distance normal to 
the original propellant surface, the exposed surface area can be predicted 
(fig. 4). Assuming that nozzle-throat erosion is reproducible, a geometric 
relationship can be determined between the ratio of propellant burning- surface 
area to nozzle-throat area and the distance regressed normal to the 

original propellant surface (fig. 5). 

The chamber pressure of a solid propellant rocket motor, with fixed 
geometry and given propellant characteristics at a known propellant 
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temperature, is a function of . The relationship between K^. , chamber 

pressure, and propellant temperature is generally determined experimentally 
and can, in some cases, be quite complex (fig. 2). However, the relationship 
between and motor chamber pressure for limited chamber-pressure 

ranges can be approximated for a given propellant at a known propellant tem- 
perature by the following empirically determined relationship (fig. 6) 


k n = bpm 


( 1 ) 


where 

b = K^j. equation coefficient. 

P - chamber pressure. 
m= K^. equation exponent. 


Burning Rate 

Burning rate r is the rate at which a solid propellant is consumed. It 
is measured in a direction normal to the propellant surface and expressed by 
in. /sec. The burning rate of a specific propellant is a function of chamber 
pressure and propellant temperature. The complex relationship between 
burning rate, chamber pressure, and propellant temperature (fig. 2) can be 
approximated for a given propellant at a known propellant temperature for 
limited chamber pressure ranges by the following empirically determined 
relationship (fig. 6) 


r 


= aP 


,n 


( 2 ) 


where 


a = burning-rate equation coefficient, 
n = burning-rate equation exponent. 
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PROPELLANT THERMAL SENSITIVITY 


Data obtained from test firings, conducted at different propellant tem- 
peratures, can be transformed to any common temperature of interest within 
the extremities of the available test data by any one of three processes. 
These processes are constant pressure, constant burn rate, and constant K^. 

Compensations must be made for factors such as erosive burning and pressure 
losses along the length of the grain perforation during motor operation. The 
constant process presented in this report is based on a constant geo- 
metric configuration at specifically regressed distances normal to the original 
propellant surface that compensates for the above mentioned factors. The 
actual regressed distance need not be known as long as the transformations 
are conducted at a constant regressed distance. This transformation method 
is also valid for the average geometric configuration during a time interval 
corresponding to a specific regressed distance (such as burn time, action 
time, and tail-off time). 


Chamber-Pressure Transformation- Equation Derivation 

The chamber-pressure transformation equation can be derived by first 
taking the natural logarithm of equation (1) giving 


In = In b + (m) In P (3) 


Taking the partial derivative of equation (3) with respect to the propellant 
temperature T at a constant , and assuming m is independent of tem- 
perature, yields 


3T ( ln S*) = 0 ST ^ ln b ^K N + ^ ST ^ ln P ^K n ^ 
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Therefore, 


3T P ^K n 


aT (In b)^ 


m 


(5) 


The expression - 5 = (In P)«. is defined as the temperature sensitivity coef- 
0T 

ficient of chamber pressure (ref. 1 ) and is expressed as . 


Performing the indicated differentiation of equation (5) yields 



Integrating equation ( 6 ) yields 



( 6 ) 


(7) 


where 


T = mean propellant temperature of the motors in the propellant tem- 
g perature group whose data are to be transformed. 

T . = desired propellant temperature or the propellant temperature to 
which the experimental test data are to be transformed. 

P. = chamber pressure that corresponds to T as obtained from a 

J D 

second-order least- squares curve fit of the experimentally 


7 



determined chamber pressures, at the specific regressed distance 
of interest normal to the original propellant surface, versus pro- 
pellant temperature (fig. 7). 

Pj = chamber pressure that corresponds to as obtained from a 

second-order least- squares curve fit of the experimentally deter- 
mined chamber pressures, at the specific regressed distance of 
interest normal to the original propellant surface, versus propel- 
lant temperature (fig. 7). 


o 

All least- squares curve fits are of the form In x = a + by + cy 
where x is the dependent parameter, y is the independent parameter, 
and a, b, and c are constants. 

The mean value theorem permits writing equation (7) as 


* f T<i 

J 


dT = 


g 


’ r p f 
/ *p 

J p 


( 8 ) 


P. 

J 




where 77^ is the average value of ir^. over the path of constant from 


Tg to . This quantity 77^ will hereafter be referred to as n 
Performing the indicated integration of equation (8) yields 


K' 


77 


K 


( T a - \) - [ 


In P f - In P. 


] 


Kn 


(9) 
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Therefore, 



By rearranging equation (10), the chamber pressure for each motor in 
the prefire propellant temperature group T g (experimentally determined at 

the specific regressed distance for which ir^r was calculated) is transformed 


to the chamber pressure corresponding to the desired prefire propellant tem- 
perature 



where 


( 11 ) 


P^ = expected chamber pressure (at the specific regressed distance of 
interest normal to the original propellant surface) of the explicit 
motor if the propellant temperature had been T^. 

P. = chamber pressure experimentally determined at the spec if ic 
regressed distance of interest normal to the original propellant 
surface of the explicit motor whose data are to be transformed. 

T. = actual propellant temperature of the explicit motor whose data are 
1 to be transformed. 
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Burning- Rate Transformation- Equation Derivation 

Reference 1 presents a relationship between m and n, m = 1 - n, 
thus equations (1) and (2) may be combined and rearranged to give 



Then, substituting equation (12) in equation (2) yields 



Taking the natural logarithm of equation (13) gives 

lnr=lna-^ r 5- 5 )lnb+( r ^- s ) In K N (14) 

Partial differentiation of equation (14) with respect to the propellant tempera- 
ture T at a constant K^., and assuming n is independent of temperature, 

yields 

ST ^ ln r ^K M = ST ^ ln a ^K M " 1 - n ST ^ ln b ^K + 1 - n ir ( ln K n) K ^ 
N N N v • N 


The 


expression J^(ln r)„ 
ST K n 


is defined as the temperature sensitivity coef- 


ficient of burning rate (ref. 1) and is expressed as cr K . 
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Performing the indicated differentiation of equation (15) yields 


°K = 



( 16 ) 


Integrating equation (16) yields 



(17) 


where 


r. 

] 


burning rate that corresponds to 



as obtained from a second- 


order least-squares curve fit of the experimentally determined 
burning rates, at the specific regressed distance of interest 
normal to the original propellant surface, versus propellant 
temperature (fig. 8). 


r^ = burning rate corresponding to as obtained from a second -order 

least-squares curve fit of the experimentally determined burning 
rates, at the specific regressed distance of interest normal to the 
original propellant surface, versus propellant temperature (fig. 8). 


The mean value theorem permits writing equation (17) as 


* 


a 


K 




(18) 
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where a R * is the average value of a K over the path of constant from 

jfe 

Tg to T d ' This ^ nantit y a K wi ^ hereafter be referred to as a K . Per- 
forming the indicated integration of equation ( 18 ) yields 


Therefore, 


K 


T , - T = In r. - 
d gj [_ f 





(19) 


( 20 ) 


By rearranging equation (20), the burning rate for each motor in the 
propellant temperature group T (experimentally determined at the specific 

regressed distance for which cr K was calculated) is transformed to the burn- 
ing rate corresponding to the desired propellant temperature 



r. e 
1 



( 21 ) 


where 


r^ = expected burning rate (at the specific regressed distance of interest 
normal to the original propellant surface) of the explicit motor if 
the propellant temperature had been T^. 

r. = burning rate experimentally determined at the specific regressed 
distance of interest normal to the original propellant surface of the 
esqjlicit motor whose data are to be transformed. 
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Time Transformation Equation Derivation 
Since burning rate is the rate at which a solid propellant is consumed, 


r 



( 22 ) 


where 


r = average burning rate during time interval t. 

w = propellant thickness consumed (measured normal to the 
original propellant surface) during time interval t. 


Assuming that the thermal expansion of the propellant is negligible and 
substituting equation (22) in equation (20) 



(23) 


where 

t. = estimated motor operating time interval at temperature T 
(fig. 9). 

tj = estimated motor operating time interval at temperature T^ 
(fig. 9). 
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By rearranging equation (23), the measured operating time interval for 
each motor in the propellant temperature group T is transformed to the 

operating time interval corresponding to the desired propellant temperature 



(24) 


where 


t. = measured operating time interval of interest of the explicit motor 
1 whose data are to be transformed. 

t, = estimated operating time interval of each motor at propellant tem- 
perature TVj. 


Thrust Related Calculations and Transformation Equation 

In order to transform longitudinal thrust in the simplest maimer and also 
to provide the capability of handling motors that utilize as many as four noz- 
zles (with equal or unequal expansion ratios) and with a nozzle cant angle (such 
as the Apollo launch escape motor), the thrust measurements are first cor- 
rected to vacuum pressure altitude 


F . = F. + P A cos 4 > 
vi i a e ^ 


(25) 


where 


F . = measured average longitudinal thrust of each motor during the 
V1 time interval of interest corrected to vacuum pressure altitude. 

F . = measured average longitudinal thrust of each motor during the 
1 time interval of interest. 


14 



P = ambient pressure experienced by each motor during the time inter- 
a val of interest. 

A = total nozzle- exit area, 
e 

c p = nozzle cant angle from the longitudinal centerline of the motor. 


Equation (25) can also be used to correct thrust data to any pressure 
altitude of interest; however, it should be noted that there has been no attempt 
to evaluate the performance data during the time interval of nozzle flow sepa- 
ration. 

The average longitudinal thrust for each motor in the propellant temper- 
ature group T (experimentally determined at, the specific regressed distance 

for which was calculated) is transformed to the average longitudinal 

Ja. 

thrust corresponding to the desired propellant temperature 



( 26 ) 


where 


F . = average expected thrust of the explicit motor if the propellant 
v temperature had been T^. 


STATISTICAL ANALYSIS 


When parameters are experimentally determined, it is desirable for 
design, performance, and reliability evaluations to establish limits or bounds 
which contain a desired percentage of a specific parameters population, with 
a confidence or probability that the intended condition is satisfied. The bounds 
thus established are called tolerance limits. 

To determine these tolerance limits, a frequency distribution must be 
assumed. This paper assumes normal distribution since experience has 
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shown that the majority of experimental data is approximately normally dis- 
tributed. Although a check for normality may be performed, moderate 
departure from the assumed distribution will not seriously affect the tolerance 
limits computed. 

Given a sample of data, the mean is estimated by 


X 


2X 

N 


( 27 ) 


and the standard deviation is estimated by 


s 


J nzx 2 - (ZX) 2 

V N(N - 1) 


( 28 ) 


where 


N = the number of points in the sample. 

X = experimentally determined values that make up the sample. 

Now that estimates of the mean and standard deviation are known, tolerance 
limits may be determined. The tolerance limits are the bounds which with 
probability y contain at least P percent of the population. 

The tolerance limits are of the form 


u = X ± Ks 


( 29 ) 


where 


K = tolerance factor. 
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Note: 


Kj = tolerance factor for one-sided tolerance limits. 
Kg = tolerance factor for two-sided tolerance limits. 


This equation may be applied in two different ways. It can be used to deter- 
mine the two-sided tolerance limits (fig. 10), or for determining either the 
upper or the lower one-sided tolerance limit (fig. 11). 


Two-Sided Tolerance Limits 

The two-sided tolerance limits are the upper and the lower bounds that 
will enclose at least the desired percentage of the expected population with a 
preselected confidence. The tolerance factors required in equation (29) are 
given in table I with confidence y for P percent of the expected population 
and N data points. 

Example: A sample of 25 points with X = 10.02 and s = 0. 13. 

Two-sided tolerance limits for 90 percent of the population with 95 percent 
confidence are desired. From table I, for P = 0.90 and y = 0.95, K = 2.208. 


u = X ± KgS = 10. 02 ± (2. 208) (0. 13) = 10. 02 ± 0. 29 


Therefore, with 95 percent confidence, at least 90 percent of the population 
lie between 9. 73 and 10. 31. 


One-Sided Tolerance Limits 

The one-sided tolerance limit is the upper or the lower bound below 
which or above which at least the desired percentage of the population can be 
expected to lie with a preselected confidence. The tolerance factors required 
in equation (29) are given in table II with confidence y for P 
(where P = 1 - a) percentage of the expected population and N data points. 
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Example: Find the one-sided tolerance limit for the preceding example. 
From table H, K= 1.838 


u = X ± = 10. 02 ± (1. 838) (0. 13) = 10. 02 ± 0. 24 

Therefore, with 95 percent confidence, at least 90 percent of the population 
lies either above 9. 78 or below 10. 26. 

The computer program presented in appendix A calculates estimated 
standard deviation, one-sided tolerance limits, and two-sided tolerance 
limits for both physical data and transformed general performance parameters 
over selected time intervals. The computer program presented in appendix B 
calculates estimated standard deviation and one-sided tolerance limits 
or two-sided tolerance limits for transformed instantaneous performance 
data versus time. 


GENERAL SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
COMPUTER PROGRAM USING THE GROUP 
TRANSFORMATION METHOD 


General Description 

This computer program was written at the Manned Spacecraft Center 
in Fortran IV for the IBM 7094 computer with 32K storage. The program 
transforms general performance parameters for selected time intervals 
acquired from as many as three different propellant temperature groups to a 
like number of propellant temperatures of interest for a solid propellant 
rocket motor of fixed geometry and a given propellant. The transformations 
are performed for the average geometric configuration during the time inter - 
vals that correspond to the specific regressed distances of interest (such as 
burn time, action time, and tail -off time). 

The actual regressed distances need not be known as long as the trans- 
formations are conducted at the average geometric configuration for the spe- 
cific regressed distances. Statistical calculations are performed on both 
physical data and general performance parameters after the parameters have 
been transformed to the specific common temperatures of interest. The 
input format for the program is presented in table HI. 
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Computer Deck Setup 


This program exceeds 32K core storage, and the overlay feature of 
Fortran IV is employed to allow the program to be run as a single input job. 

Six decks constitute the program in its entirety and include a calling 
program and five subroutines. 

The actual deck setup for the IBM 7094 is illustrated in figure 12. 

1. Setup Cards 

$JOB card 
$LBJOB card 

2. Program 

DRIVER 

STATS 

$ORIGIN ALPHA card 

MOTORS 

GLS1 

$ORIGIN ALPHA card 

MOT 2 

MOT3 

3. Data 

$DATA 
DATA deck 
7/ 8 card 

DRIVER . - DRIVER is the executive routine required in an overlayed 
program. This routine never leaves the machine storage area and calls the 
various overlayed subroutines. 

Subroutine STATS . - Subroutine STATS computes the mean and the stan- 
dard deviation of an input array of given size. Due to the computer arithmetic 
methods, standard deviations were sometimes found to exist with a small non- 
zero value when all elements of the input array were identical. To avoid this 
difficulty, standard deviation is set to zero if its ratio to the mean is less 
than 0. 00025. Since STATS is used by all overlayed sections of the program, 
it has been placed in the permanent storage area with DRIVER. 
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Subroutine MOTORS. - Subroutine MOTORS reads in all the data. It 
sorts the data for each motor into the proper propellant temperature group 
and does performance calculations on the input data. It calculates the trans- 
formation factors for pressures and times (tt 

transformations on these data. 

Subroutine GLS1. - Subroutine GLS1 is used to calculate a second-order 
least- squares curve fit of the natural logarithms of the chamber pressures 
and the operation time intervals versus propellant temperature. 

Subroutine MOT2 . - Subroutine MOT2 calculates transformed thrusts, 
impulses and specific impulses, and the tolerance limits at each desired 
pressure altitude. 

Subroutine MOT3 . - Subroutine MOT3 calculates transformed burning 
rate, characteristic velocity, and their tolerance limits. 


K and cr^ and performs the 


Program Restrictions 

This program was written to perform the desired transformations for a 
maximum of 15 motors per experimental prefire propellant temperature group, 
3 experimental prefire propellant temperature groups, and 9 time intervals 
such as ignition delay time, thrust rise time, burn time, total time, and so 
forth. 


The input restrictions are that the number of propellant temperatures 
(to which the data are to be transformed) must equal the number of experi- 
mental temperature groups, the first ambient pressure (EXPA (1)) must 
equal zero psia (vacuum), and the desired temperatures must be in ascending 
order. 

Care must be exercised in the evaluation of the program results. For 
example, outputs of specific impulse, total impulse, and chamber pressure 
integral are meaningless if instantaneous data are inputed. The evaluation 
of data acquired during nozzle flow separation is beyond the scope of this pro- 
gram ; therefore, action time, tail-off time, and total time data must be care- 
fully handled. It should be noted that the output values for characteristic 
velocity (CSTAR) are valid for total time only. 


I 
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Output Formats 


The output formats are basically of two types : transformed data and 
nontransformed data. Table IV presents a typical page of nontransformed 
data output. The type of output (time, pressure, thrust, or other parameters) 
is identified in the upper left corner of the page. The nine columns represent 
the motor operating time intervals of interest (such as ignition delay time, 
thrust rise time, burn time, total time, and so forth). Nine columns will 
always be printed out even though fewer time intervals are inputed. The non- 
used columns will be printed out as zeros. At the left of the page, each motor 
is identified by its alphanumeric identification code, and the motors are sepa- 
rated into propellant temperature groups (up to three). Means and standard 
deviations are printed out for the motors in each temperature group. In three 
instances (thrust, impulse, and specific impulse) identical formats are used 
for different values of pressure altitude. The pressure altitude is printed out 
at the top of the page. 

Table V presents a typical output page of transformed data. Each out- 
put page represents one time interval such as ignition delay time, thrust rise 
time, burn time, total time, and so forth. The time interval is printed out 
in numerical form (time no. 1, time no. 2, time no. 3, and up to time no. 9) 
in the upper left corner of the page and must be correlated with the time defi- 
nitions in the input data (time no. 1 = ignition delay time, time no. 2 = thrust 
rise time, time no. 3 = burn time, or up to time no. 9 = total time). The type 
of output (time, pressure, thrust, and other outputs) is printed out at the top 
of the page. The three columns represent the three propellant temperatures 
to which the input data have been transformed. At the left of the page each 
motor is identified by its alphanumeric identification code, and the motors are 
separated into the original propellant temperature groups (up to three). Means 
and standard deviations are printed out for motors in each temperature group. 
Means, standard deviations, and tolerance limits are printed out for the three 
groups as a whole. In three instances (thrust, impulse, and specific impulse) 
identical formats are used for different values of pressure altitude. The pres- 
sure altitude is printed out under the type of output. 
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SOLID PROPELLANT ROCKET MOTOR PERFORMANCE VERSUS 
TIME COMPUTER PROGRAM USING THE GROUP 
TRANSFORMATION METHOD 


General Description 

This computer program was written at the Manned Spacecraft Center 
in Fortran IV for the Univac 1107 computer with 65K storage. The program 
effectively transforms experimentally obtained thrust-and-chamber pressure 
versus time data (acquired for as many as three different propellant tempera- 
ture groups) to a common propeUant temperature of interest for a solid 
propellant rocket motor of fixed geometry and a given propellant. The 
transformations are performed at specific regressed distances normal to 
the original propeUant surface. The actual regressed distance need not be 
known, as long as the transformations are conducted at a constant regressed 
distance (such as percents of burn time and percents of tail-off time). See 
figure 13. 

Statistical calculations are performed on the test data after they have 
been transformed to the specific common temperature of interest. The input 
format for the program is presented in table VI. 


Computer Deck Setup 

Twelve decks constitute the program in its entirety, including a main 
program and 12 subroutines. The program for the Univac 1107 is illustrated 
in figure 14. 

1. Setup Card 
$JOB card 

2. Program 

MAIN 

SRCH1 

LSTSQ 

CURVE 

STATS 

GLS1 

ACCEND 

NORMLZ 

FEXIT 
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STDEV 

CF2F1 

CF2F2 

QUTKMV 

3. Data 

XQT card 
DATA deck 
EOF card 
FIN card 

MAIN program . - The MAIN program reads in all data and sorts the data 
for each motor into its proper prefire propellant temperature group. It then 
calls the CURVE FIT routines to curve fit both input pressure and input thrust 
versus input time. It subdivides burn time and tail-off time into the required 
number of subintervals (fig. 13). Finally, it calls subroutine LSTSQ to com- 
pute the transformed data. 

Subroutine SRCH1 . - Subroutine SRCH1 computes time increments in 
terms of percent web time and percent tail-off time (fig. 13). The subroutine 
then selects from the CURVE FIT routines the proper cubic coefficients for 
the particular time increment being considered. After the proper coefficients 
have been selected, the subroutine uses them to compute the values of the 
ordinate (pressure or thrust) at each time increment and returns these values 
to the main program. 

Subroutine LSTSQ . - Subroutine LSTSQ computes the transformation 
factors for pressures and times and ct k ) and performs the transforma- 
tions on these and thrust data at each percent web time and at each percent 
tail-off time and prints them out as final answers. Graphs of the final answers 
are also produced (see subroutine QUIKMV). 

Subroutine STATS . - Subroutine STATS computes the mean and standard 
deviation of an input array of given size. Due to computer ar ithm etic 
methods, standard deviations were sometimes found to exist with a small 
nonzero value when all elements of the input array were identical. To avoid 
this difficulty, standard deviation is set to zero if its ratio to the mean is less 
than 0. 00025. 

Subroutine GLS1 . - Subroutine GLS1 is used to calculate a second-order 
least- squares curve fit of the natural logarithms of chamber pressures, and 
motor operation times (at each experimentally determined percent web time 
and percent tail-off time) versus propellant temperature. 
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CURVE FIT routines. - The CURVE FIT routines consist of subroutines 
CURVE, ACCEND, NORMLZ, FOOT, STDEV, CF2F1, and CF2F2. These 
subroutines fit an input ordinate array versus an input abscissa array using a 
piecewise cubic least- squares curve fit. Because of the piecewise nature of 
the fit there are NPTS-ICON +1 sets of cubic coefficients, where NPTS is the 
size of the input arrays, and ICON points are fit in each piece. Each coeffi- 
cient set is valid over a limited interval. The coefficients may be printed out 
by setting input parameter IPRNT to 1. 

Subroutine QUIKMV. - Subroutine QUIKMV is a MSC general plot routine 
for the SC-4020. It has not been included in appendix B since most facilities 
use plotting routines individually suited to the available plotting equipment. 

The call to this routine may be deleted, or a dummy routine substituted, with- 
out affecting the overall program. 


Program Restrictions 

This program was written to perform the desired transformation for as 
many as 15 motors per experimental propellant temperature group, 3 experi- 
mental propellant temperature groups, and 200 selected inputs of time with 
corresponding chamber pressure and thrust for each motor. 

Since the program is limited to 200 inputs of time with corresponding 
chamber pressure and thrust for each motor, care must be exercised in 
choosing the input points in order that they best define the variation of chamber 
pressure and thrust versus time for each motor (fig. 15). Care must also be 
exercised in choosing the proper output points that best define the variations 
of chamber pressure and thrust versus time. This can be accomplished since 
burn time can be subdivided into as many as five subintervals, and tail-off 
time can be subdivided into as many as four subintervals. As many time 
increments (percents of web time and/or percents of tail- off time) as desired 
can be generated to define the variation of chamber pressure and thrust ver- 
sus time for each subinterval as long as the total does not exceed 200 points 
(fig. 13). 


Output Formats 

The outputs consist of the transformed results in computer printout 
form, such as presented in tables VII and VIH, and SC-4020 graphs as shown 
in figures 16(a) and 16(b). 
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The transformed results present the statistics (mean and tolerance 
limits) of the transformed time, chamber pressure, and thrust at the desired 
percent of burn time and percent tail-off time. Transformed thrust statistics 
can be presented for as many as three desired pressure altitudes. 


CONCLUDING REMARKS 


The method and computer programs presented in this report were de- 
veloped at the Manned Spacecraft Center specifically to evaluate the perform- 
ance of the solid rockets used in the Apollo Launch Escape System. From the 
conception, however, the method of solution, including the machine programs, 
was intended to be generally applicable for any rocket motor utilizing a single 
propellant grain with no throttling capability. 

The programs were designed for two types of evaluations. Often, in the 
requirements, motors are constrained to operate within a given performance 
regime during some portion of the firing, while delivering a specified nominal 
performance over the total operation. The program presented in appendix A 
was intended to evaluate the more stringent constraints placed upon the launch 
escape motor where both the maximum thrust level and the minimum impulse 
were specified for an initial phase of firing, and an overall minimum perform- 
ance required for the abort mission. 

The program presented in appendix B compares the transformed thrust- 
time and pressure-time relationships of all motors tested to define the ex- 
pected limits of the performance at any temperature. This program has been 
used to describe the general thrust and pressure characteristics for documen- 
tation and product variance purposes. 

It is anticipated that these programs will be valuable, particularly for 
contracting agencies and prime contractors, in the evaluation of solid rocket 
test data. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, June 29, 1966 
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TABLE I. - TWO-SIDED TOLERANCE FACTORS FOR NORMAL DISTRIBUTIONS 11 

Factors K such that the probability is y that at least a proportion P of the 
distribution will be included between X ± Ks, where X and s are estimates 
of the mean and the standard deviation computed from a sample of N 
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From Techniques of Statistical Analysis by Eisenhart, Hastay, and Wallis. 
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TABLE I. - TWO-SIDED TOLERANCE FACTORS 
FOR NORMAL DISTRIBUTIONS 21 - Concluded 


Factors K such that the probability is y_ that at least ajproportion P of the 
distribution will be included between X ± Ks, where X and s are estimates 
of the mean and the standard deviation computed from a sample of N 
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3.364 

4.297 

63 

1.251 

1.789 

2.131 

2.801 

3.577 

1.334 

1.907 

2.272 

2.985 

3.813 

1.387 

1.984 

2.363 

3.106 

3.967 

1.498 

2.142 

2.552 

3.354 

4.284 

64 

1.250 

1.787 

2.129 

2.798 

3.574 

1.331 

1.904 

2.268 

2.981 

3.807 

1.384 

1.980 

2.359 

3.100 

3.959 

1.494 

2.136 

2.545 

3.344 

4.272 

65 

1.249 

1.785 

2.127 

2.795 

3.571 

1.329 

1.901 

2.266 

2.976 

3.801 

1.382 

1.976 

2.354 

3.094 

3.951 

1.490 

2.130 

2.538 

3.335 

4.260 

66 - 

a. 247 

1.784 

2.125 

2.793 

3.567 

1.327 

1.898 

2.261 

2.972 

3.796 

1.379 

1.972 

2.350 

3.088 

3.944 

1.486 

2.124 

2.531 

3.326 

4.249 

67 

1.246 

1.782 

2.123 

2.790 

3.564 

1.325 

1.895 

2.258 

2.967 

3.790 

1.377 

1.968 

2.345 

3.082 

3.937' 

1.482 

2.119 

2.524 

3.318 

4.238 

68 

1.245 

1.781 

2.122 

2.788 

3.661 

1.323 

1.892 

2.255 

2.963 

3.785 

1.374 

1.965 

2.341 

3.076 

3.930 

1.478 

2.113 

2.518 

3.309 

4.227 

66 

1.244 

1.779 

2.120 

2.786 

3.558 

1.322 

1.890 

2.252 

2.959 

3.779 

1.372 

1.961 

2.337 

3.071 

3.923 

1.474 

2.108 

2.512 

3.301 

4.216 

60 

1.243 

1.778 

2.118 

2.784 

3.556 

1.320 

1.887 

2.248 

2.955 

3.774 

1.369 

1.958 

2.333 

3.066 

3.910 

1.471 

2.103 

2.506, 

3.293 

4.206 

61 

1.242 

1.776 

2.117 

2.781 

3.553 

1.318 

1.885 

2.245 

2.951 

3.769 

1.367 

1.955 

2.329 

3.061 

3.909 

1.467 

2.098 

2.500 

3.285 

4.196 

es 

1.141 

l.m 

6,115 

6.776 

5.550 

1.110 

1.882 

3,248 

3,947 

8.785 

1.885 

1.961 

9,826 

8.066 

1.903 

1.464 

2,098 

3.494 

8,278 

4.187 

66 

1.840 

1.774 

6.116 

6.777 

6.548 

1.115 

1.880 

1.240 

3.944 

8.760 

1.888 

1.948 

9.829 

8.081 

8.897 

1.461 

2 089 

1480 

s.sn 

4.178 

64 

1.240 

1.772 

6.112 

6.775 

5 . 54 ^ 

1.818 

1 

1.878 

3.237 

2.940 

8.755 

1.861 

1.946 

2.318 

8.046 

8.891 

1.468 

2.064 

3.483 

8.884 

4. 108 


a From Techniques of Statistical Analysis by Eisenhart, Hastay, and Wallis. 
Copyright 1947. McGraw-Hill Book Company. Used by permission. 



TABLE n. - ONE-SIDED TOLERANCE FACTORS FOR NORMAL DISTRIBUTIONS ' 1 

Factors K such that the probability's 7 that at least a proportion 1 - a of 
tae distribution will be less than X ± Ks (or greater than X - Ks), where 
X and s are estimates of the mean and standard deviation computed from 
_ a sample of size N 


'X a 
N \ 

0.25 

7= 0.75 
0.10 0,05 

0. 01 

0. 001 

0. 25 

.y = 0.90 
0. 10 0. 05 

0. 01 

0. 001 

0. 25 

y 

0.10 

= 0.95 
0. 05 

0. 01 

0. 001 

0. 25 

7 

0.10 

= 0.99 
0.05 

0. 01 

0. 001 

3 

1.464 

2.501 

3.152 

4.396 

5.805 

2.602 

4.258 

5.310 

7.340 

9.651 

3.804 

6.158 

7.655 

10.552 

13.857 






4 

1.256 

2.134 

2.680 

3.726 

4.910 

1.972 

3.187 

3.957 

5.437 

7.128 

2.619 

4.163 

5.145 

7.042 

9,215 






5 

1.152 

1.961 

2.463 

3.421 

4.507 

1.698 

2.742 

3.400 

4.666 

6.112 

2.149 

3.407 

4. 202 

5.741 

7.501 






6 

1.087 

1.860 

2.336 

3.243 

4.273 

1.540 

2.494 

3.091 

4. 242 

5.556 

1.895 

3.006 

3.707 

5.062 

6.612 

2.849 

4.408 

5.409 

7.334 

9.540 

7 

1.043 

1.791 

2.250 

3.126 

4.118 

1.435 

2.333 

2.894 

3.972 

5.201 

1.732 

2.755 

3.399 

4.641 

6.061 

2.490 

3.856 

4.730 

6.411 

8.348 

8 

1.010 

1.740 

2.190 

3.042 

4. 008 

1.360 

2.219 

2.755 

3.783 

4.955 

1.617 

2.582 

3.188 

4.353 

5.686 

2.252 

3.496 

4.287 

5. 811 

7.566 

9 

0. 984 

1.702 

2.141 

2.977 

3.924 

1.302 

2. 133 

2.649 

3.641 

4.772 

1.532 

2.454 

3.031 

4. 143 

5.414 

2.085 

3.242 

3.971 

5.389 

7. 014 

10 

0. 964 

1.671 

2.103 

2.927 

3.858 

1.257 

2. 065 

2.568 

3.532 

4.629 

1.465 

2.355 

2.911 

3,981 

5.203 

1.954 

3.048 

3.739 

5.075 

6. 603 

11 

0. 947 

1.646 

2.073 

2.885 

3.804 

1.219 

2.012 

2.503 

3.444 

4.515 

1.411 

2.275 

2.815 

3.852 

5.036 

1.854 

2.897 

3.557 

4. 828 

6. 284 

12 

0. 933 

1.624 

2.048 

2.851 

3.760 

1.188 

1.966 

2.448 

3.371 

4.420 

1.366 

2.210 

2.736 

3.747 

4. 900 

1.771 

2.773 

3.410 

4. 633 

6. 032 

, 13 

0.919 

1.606 

2.026 

2.822 

3.722 

1.162 

1.928 

2.403 

3.310 

4.341 

1.329 

2.155 

2.670 

3.659 

4.787 

1.702 

2. 677 

3.290 

4. 472 

5. 826 

14 

0. 909 

1.591 

2.007 

2.796 

3.690 

1.139 

1.895 

2.363 

3.257 

4.274 

1.296 

2.108 

2.614 

3.585 

4. 690 

1.645 

2.592 

3.189 

4.336 

5.651 

15 

0. 899 

1.577 

1.991 

2.776 

3.661 

1.119 

1.866 

2.329 

3.212 

4.215 

1.268 

2.068 

2.566 

3. 520 

4. 607 

1.596 

2. 521 

3.102 

4.224 

5.507 

16 

0. 891 

1.566 

1.977 

2.756 

3.637 

1.101 

1. 842 

2.299 

3.172 

4. 164 

1.242 

2.032 

2. 523 

3.463 

4. 534 

1.553 

2.458 

3.028 

4. 124 

5.374 

17 

0. 883 

1.554 

1.964 

2. 739 

3.615 

1.085 

1.820 

2.272 

3.136 

4.118 

1.220 

2.001 

2.486 

3.415 

4. 471 

1.514 

2.405 

2. 962 

4. 038 

5.268 

18 

0. 876 

1.544 

1.951 

2.723 

3.595 

1.071 

1.800 

2.249 

3.106 

4. 078 

1.200 

1.974 

2.453 

3.370 

4.415 

1.481 

2.357 

2.906 

3.961 

5.167 

19 

0. 870 

1.536 

1.942 

2.710 

3.577 

1.058 

1.781 

2.228 

3.078 

4. 041 

1.183 

1.949 

2. 423 

3.331 

4.364 

1.450 

2.315 

2.855 

3.893 

5.078 

20 

0. 865 

1.528 

1.933 

2.697 

3.561 

1.046 

1.765 

2.208 

3.052 

4. 009 

1.167 

1.926 

2.396 

3.295 

4.319 

1.424 

2. 275 

2.807 

3.832 

5.003 

21 

0. 859 

1.520 

1.923 

2.686 

3.545 

1.035 

1.750 

2.190 

3.028 

3.979 

1.152 

1.905 

2. 371 

3.262 

4.276 

1.397 

2. 241 

2.768 

3.776 

4. 932 

22 

0. 854 

1.514 

1.916 

2,675 

3.532 

1.025 

1.736 

2.174 

3.007 

3.952 

1.138 

1.887 

2.350 

3.233 

4. 238 

1.376 

2. 208 

2.729 

3.727 

4. 866 

23 

0. 849 

1.508 

1.907 

2.665 

3.520 

1.016 

1.724 

2.159 

2.987 

3.927 

1.126 

1.869 

2.329 

3.206 

4.204 

1.355 

2. 179 

2.693 

3.680 

4. 806 

24 

0. 845 

1.502 

1.901 

2. 656 

3.509 

1.007 

1.712 

2. 145 

2.969 

3.904 

1.114 

1.853 

2.309 

3.181 

4.171 

1.336 

2.154 

2.663 

3.638 

4.755 

25 

0.842 

1.496 

1.895 

2, 647 

3.497 

0.999 

1.702 

2.132 

2.952 

3.882 

1.103 

1.838 

2.292 

3.158 

4.143 

1.319 

2.129 

2.632 

3.601 

4. 706 

30 

0..825 

1.475 

1.869 

2.613 

3.454 

0. 966 

1.657 

2.080 

2.884 

3.794 

1.059 

1.778 

2.220 

3.064 

4. 022 

1.249 

2.029 

2.516 

3.446 

4. 508 

35 

0. 812 

1.458 

1.849 

2.588 

3.421 

0. 942 

1.623 

2.041 

2. 833 

3.730 

1. 025 

1.732 

2.166 

2.994 

3.934 

1.195 

1.957 

2.431 

3.334 

4.364 

40 

0. 803 

1.445 

1.834 

2.568 

3.395 

0. 923 

1.598 

2.010 

2. 793 

3.679 

0. 999 

1.697 

2.126 

2.941 

3.866 

1.154 

1. 902 

2.365 

3.250 

4.255 

45 

0. 795 

1.435 

1.821 

2.552 

3.375 

0. 908 

1.577 

1.986 

2.762 

3.638 

0. 978 

1.669 

2.092 

2.897 

3.811 

1.122 

1.857 

2.313 

3.181 

4.168 

50 

0.788 

1.426 

j 

1.811 

2.538 

3.358 

0. 894 

1.560 

1.965 

2. 735 

3.604 

0. 961 

1.646 

2.065 

2, 863 

3.766 

1.096 

1.821 

2.296 

3,124 

4.096 


a From Engineering Statistics by Albert H. Bowker and Gerald J. Lieberman. Copyright 1959. Prentice -Hall, Inc. Used by permission. 



TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM 


Card 

number 

(a) 

Variable 

name 

Format 

Columns 

Description 

Units 

1 

CASE 

A6 

1-6 

Case number being 
evaluated 

- 

2 

HEADING 

80H 

2-80 

Designation of the 
motor being evaluated 

- 

3 

NTIMES 

12 

11-12 

Integer designating the 
number of time inter- 
vals being evaluated 
(maximum of 9) 

- 

4 

Description of time 
interval no. 1 

8 OH 

2-80 

Definition of regressed 
distance no. 1 (such as 
ignition delay time) 

- 

5 

Description of time 
interval no. 2 

80H 

2-80 

Definition of regressed 
distance no. 2 (such as 
thrust rise time) 

- 

6 

Description of time 
interval no. 3 

80H 

2-80 

Definition of regressed 
distance no. 3 (such as 
ignition time) 

- 

7 

Description of time 
interval no. 4 

8 OH 

2-80 

Definition of regressed 
distance no. 4 (such as 
burn time) 

- 

8 

Description of time 
interval no. 5 

80H 

2-80 

Definition of regressed 
distance no. 5 (such as 
action time) 

- 

9 

Description of time 
interval no. 6 

80H 

2-80 

Definition of regressed 
distance no. 6 (such as 
tail -off time) 

- 

10 

Description of time 
interval no. 7 

80H 

2-80 

Definition of regressed 
distance no. 7 (such as 
total time) 

- 

11 

Description of time 
interval no. 8 

8 OH 

2-80 

Definition of regressed 
distance no. 8 

- 

12 

Description of time 
interval no. 9 

8 OH 

2-80 

Definition of regressed 
distance no. 9 

- 


a Cards 4 through 12 are shown; however, only NTIMES cards are required. 



TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 


■ ■■IIIIIM 111 


ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued 


Card 

number 

Variable 

name 

Format 

Columns 

Description 

Units 

13 

NP 

112 

12 

Integer designating 
the number of pressure 
altitudes at which the 
performance data are to 
be reported (up to 4) 



EXPA(l) b 

F12.4 

13-24 

First pressure altitude 
at which the perform- 
ance data are to be 
reported 

psia 


EXPA(2) 

F12.4 

25-36 

Continue for each indi- 
vidual pressure altitude 
at which the perform- 
ance data are to be 
reported (up to 4) 

psia 


EXPA(3) 

F12.4 

37-48 


psia 


EXPA(4) 

F12.4 

49-60 


psia 

14 

NT C 

112 

12 

Integer designating the 
number of prefire pro- 
pellant temperatures to 
which the performance 
data are to be trans- 
formed (up to 3) 



EXTEM(l) d 

F12.4 

13-24 

First prefire propel- 
lant temperature to 
which the performance 
data are to be trans- 
formed 

°F 

i 

EXTEM(2) d 

F12.4 

25-36 

Continue for each individ- 
ual prefire propellant 
temperature to which the 
performance data are to 
be transformed (up to 3) 

°F 


EXTEM(3) d 

F12.4 

37-48 


°F 


b The first pressure altitude (at which the performance data are to be reported) must be 
zero psia (vacuum). 


c The number of prefire propellant temperatures (to which the data are to be trans- 
formed) must equal the number of experimental temperature groups; therefore, NT on 
card 14 must equal NT on card 15. 

d The value of EXTEM(3) must be greater than EXTEM(2) and EXTEM(2) must be greater 
than EXTEM(l). 
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TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued 


Card 

number 

Variable 

name 

Format 

Columns 

Description 

Units 

15 

NT C 

112 

12 

Integer designating the 
number of prefire pro- 
pellant temperature 
group (up to 3) 



TGRP(lf 

F12. 4 

13-24 

First prefire propellant 
temperature about which 
the motors are grouped 
(prefire propellant tem- 
perature group no. 1) 

°F 


TGRP(2) e 

F12.4 

25-36 

Continue for each indi- 
vidual prefire propel- 
lant temperature about 
which the motors are 
grouped (up to 3) 

°F 


TGRP(3) e 

F12.4 

37-48 


°F 

16 

G 

F12. 5 

1-12 

Acceleration of gravity 

ft/sec^ 


PHI 

F12.5 

13-24 

Nozzle cant (from the 
longitudinal center line 
of the motor) 

rad 


NTOT 

12 

25-26 

Integer designating the 
time interval that is 
total time 

— 


NBURN 

12 

27-28 

Integer designating the 
time interval that is 
burn time 

- 


NFIRST* 

12 

29-30 

. 

Integer designating the 
first time interval that 
all data are to be evalu- 
ated (see note) 



c The number of prefire propellant temperatures (to which the data are to be trans- 
formed) must equal the number of experimental temperature groups; therefore, NT on 
card 14 must equal NT on card 15. 


e The value of TGRP(3) must be greater than TGRP(2), and TGRP(2) must be greater than 
TGRP(l). 

*In some cases (as in ignition transient data, for example) it is desirable to transform 
only the time data. This is accomplished in the program by inputing ’’time only” data in 
the first few time intervals, and setting the value of NFERST to indicate the first time 
interval that all data are to be evaluated. 



TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued 


Card 

number 

fe) 

Variable 

name 

Format 

Columns 

Description 

Units 

17 

CT1 

F 12. 4 

1-12 

One-sided tolerance 
factor for transformed 
data 

- 


CT2 

F12.4 

13-24 

Two-sided tolerance 
factor for transformed 
data 

- 


CPI 

F12.4 

25-36 

One-sided tolerance 
factor for nontrans- 
formed data 

- 


CP2 

F12.4 

36-48 

Two-sided tolerance 
factor for nontrans- 
formed data 

- 

18 

XMOT 

A6 

1-6 

Motor number 

- 


FIRETP 

F12. 4 

7-18 

Prefire propellant 
temperature 

°F 


ENDMOT 

12 

19-2 

Integer designating that 
the last motor to be pro- 
cessed has been reached 
(1 for last motor, 0 for 
all other motors) 


19 

ATI 

E12. 5 

1-12 

Throat area of nozzle 
no. 1 

. 2 
in. 


AT2 

E12. 5 

13-24 

Continue for each indi- 
vidual nozzle (up to 4) 

, 2 
in. 


AT3 

E12. 5 

25-36 


. 2 
in. 


AT4 

E12. 5 

37-48 


. 2 
in. 


®Cards 18 through 26 are required for each individual motor. 
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TABLE HI. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued. 


Card 

number 

(g) 

Variable 

name 

Format 

Columns 

Description 

Units 

20 

AE1 

E12.5 

1-12 

Exit area of nozzle 
no. 1 

. 2 
in. 


AE2 

E12.5 

13-24 

Continue for each indi- 
vidual nozzle (up to 4) 

. 2 
in. 


AE3 

E12. 5 

25-36 


2 

in. 


AE4 

E12. 5 

37-48 


. 2 
in. 

21 

WP 

E12. 5 

1-12 

Total propellant weight 

lbf 


WEB 

E12. 5 

13-24 

Web thickness 

in. 


DENS 

E12. 5 

25-36 

Propellant density 

CO 

g 


PAF 

E12.5 

37-48 

Ambient pressure at the 
time of motor firing 

psia 

22 

Time 1 

F9.4 

1-9 

Time interval required 
to achieve regress dis- 
tance no. 1 

sec 


Time 2 

F9. 4 

10-18 

Continue for each indi- 
vidual time interval 
being evaluated (up to 
NTIMES) 

sec 


Time 3 

F9.4 

19-27 


sec 


Time 4 

F9.4 

28-36 


sec 


Time 5 

F9.4 

37-45 


sec 


Time 6 

F9.4 

46-54 


sec 


Time 7 

F9.4 

55-63 


sec 


Time 8 

F9.4 

64-72 


sec 


Time 9 

F8.4 

73-80 


sec 


'’Cards 18 through 26 are required for each individual motor. 
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TABLE El. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued 


Card 

number 

(g) 

Variable 

name 

(h) 

Format 

Colurms 

Description 

Units 

23 

INTP1 


F9.4 

1-9 

Measured chamber pres- 
sure integral during 
Time interval 1 

psia- 

sec 


INTP2 


F9.4 

10-18 

Continue for each indi- 
vidual time interval 
being evaluated (up to 
NTIMES) 

psia- 

sec 


INTP3 


F9.4 

19-27 


psia- 







sec 


INTP4 


F9.4 

28-36 


psia- 

sec 


INTP5 


F9.4 

37-45 


psia- 

sec 


INTP6 


F9.4 

46-54 


' psia- 
sec 


INTP7 


F9.4 

55-63 


psia- 







sec 


INTP8 


F9.4 

64-72 


psia- 







sec 


INTP9 


F8. 4 

73-80 


psia- 







sec 

24 

PI 


F9.4 

1-9 

Measured average cham- 
ber pressure during 
Time interval 1 

psia 

, 

P2 


F9.4 

L 

10-18 

Continue for each indi- 
vidual time interval 
being evaluated (up to 
NTIMES) 

psia 


£ 

& Cards 18 through 26 are required for each individual motor. 

Average or integral data (but not both) are required for a particular time interval; 
however, cards 23, 24, 25, and 26 are required because average and integral data may be 
intermingled with the various time intervals. 
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TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Continued 


Card 

number 

(g) 

Variable 

name 

00 

Format 

Columns 

* 

Description 

Units 

24 

P3 

F9.4 

19-27 

Continue for each indi- 
vidual time interval 
being evaluated ( up to 
NTIMES) 

psia 


P4 

F9.4 

28-36 


psia 


P5 

F9.4 

37-45 


psia 


P6 

F9.4 

46-54 


psia 


P7 

F9. 4 

55-63 


psia 


P8 

F9. 4 

64-72 


psia 


P9 

F8. 4 

73-80 


psia 

25 

INTF1 

F9. 4 

1-9 

Measured impulse during 
Time interval 1 

lbf- 

sec 


INTF2 

F9.4 

10-18 

Continue for each indivi- 
dual time interval being 
evaluated (up to NTIMES) 

lbf- 

sec 


INTF3 

F9. 4 

19-27 


lbf- 

sec 


INTF4 

! 

F9. 4 

28-36 


lbf- 

sec 


INTF5 

F9.4 

37-45 


lbf- 

sec 


INTF6 

F9. 4 

46-54 


lbf- 

sec 


INTF7 

F9.4 

55-63 


Ibf- 

sec 


^Cards 18 through 26 axe required for each individual motor. 

^Average or integral data (but not both) are required for a particular time interval; 
however, cards 23, 24, 25, and 26 are required because average and integral data may be 
intermingled with the various time intervals. 
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TABLE m. - INPUT FORMAT FOR THE GENERAL SOLID PROPELLANT 
ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM - Concluded 


Card 

number 

(g) 

Variable 

name 

(h) 

Format 

Columns 

Description 

Units 

25 

INTF8 

F9.4 

64-72 

Continue for each indi- 
vidual time interval 
being evaluated (up to 
NTIMES) 

Ibf- 

sec 


INTF9 

F8.4 

73-80 


lbf- 

sec 

26 

FI 

F9.4 

1-9 

Measured average thrust 
during Time interval 1 

lbf 


F2 

F9.4 

10-18 

Continue for each indi- 
vidual time interval 
being evaluated (up to 
(NTIMES) 

lbf 


F3 

F9.4 

19-27 


lbf 


F4 

F9.4 

28-36 


lbf 


F5 

F9.4 

37-45 


lbf 


F6 

F9. 4 

46-54 


lbf 


F7 

F9.4 

55-63 


lbf 


F8 

F9.4 

64-72 


lbf 


F9 

F8. 4 

73-80 


lbf 


g Cards 18 through 26 are required for each individual motor. 

^Average or integral data (but not both) are required for a particular time interval; 
however, cards 23, 24, 25, and 26 are required because average and integral data may be 
intermingled with the various time intervals. 
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TABLE IV. - TYPICAL NONTRANSFORMED DATA OUTPUT FORMAT FOR THE GENERAL SOLID 


PROPELLANT ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM 


IMPULSE 


PA = 000 







MOTOR NO. 

IT1 

IT2 IT3 

IT4 

IT 5 

IT6 

IT7 

IT8 

IT9 


TE MP GRO UP .. X 

-37A .00 . tO P 479. .00 1.51 Q ..Q7 I75$,fl8 1 78 1 .6 1 271.11 *i!0 *m 


70 

.00 

.00 

730.06 

1515.86 

1728.77 

1752.70 

236.41 

.00 

.00 

44A 

.00 

.00 

_630i.9_6.__ 

I438i^5_. 

1757,28 

_ 1780. Q2„ 

280.95 

__.iO o . 

.00 

45A 

.00 

.00 

544.38 

1535.37 

1772.00 

1795.93 

260.17 

.00 

.00 

43A 

.00 

• 00 

^49.86 

15g5.3fl 

1771.66 

1798.37 

272. 6fl .... _ 

• 00 

.00. 

74A 

.00 

.00 

688.55 

1556.49” 

1768.54 

1794,93 

238.10 

.00 

.00 

MEAN 

• 00 

.00 

603.80 

1523.64 

1759.02 

1783.93 

259.89 

.00 

.00 


xOO tUQ 95. 79 — . .. 20*4-4. 16*4J .. 12 * 15 . 1B_*75- .00 .00 


TEMP GROUP 2 


77 

.00 

• 00 

552.44 

14Q6.92 

1754.25 

1779.85 

282. 5q 

.on 

.00_ 

50A 

.00 

• 00 

780.61 

1532,28 

1771.82 

1797.06 

264,50 

• 00 

.00 

38A _ 

i-QO 



742.15 

1470.46 

1723.1Q . . 

_ 1751,85 

-.^281.05- 

.00 . 

*00 

59A 

.00 

.00 

708.81 

1544,61 

1765,84 

1787.62 

242.64 

• 00 

.00 

MEAN 

.00 

.00 

696.00 

1511.07 

1753,75 

1779.09 

267,69 

.00 

*00 

STAND, qEV. 

.on 

.on 

inn. li 

33 T 78 

21.70 

19.48 

18^6£_ 

_ .no 

t pp 

TEMP ftRQUP -3 










69 

iQO 

.00 

848.42 

1485j_66_ _ 

. .1755*91 . 

1778.68 

292,68 

.00 

.00 

64A 

.00 

• 00 

620.51 

1537.45 

1797.06 

1820.64 

282.93 

.00 

.00 

58A 

.00 

.00 

860.31 

lbn8.21 

1801 . 22 

1823.41 

314.78 

*jQ0 

.00 

65A 

.00 

• 00 

705.46 

1521.27 

1795.95 

1811.99 

290.27 

.00 

.00 

4Q_A 

.00 

iM. ... 

.92,60. _ 

__1505*_43 

17.7*. 3.5._ . 

1800,10 . . 

„ 294*36 

• 00 

.00 

78 

.00 

.00 

682.98 

1513.31 

1747.85 

1772.93 

259.16 

.00 

• 00 

. 60A 

*00 

.P_Q 

796*89 _ 

1503,47 . 

. 1768,25 

1794.46 

290.66 

.00 

.00 

MEAN 

.00 

.00 

729.60 

1Rin . fefl 

1777.23 

1800.32 

289.26 

QQ_ ... 

__ juOQ 

STAND. DEV. 

.00 

• 00 

107.36 

16.06 

21.31 

19.73 

16.52 

• 00 

• 00 





TABLE V.- TYPICAL TRANSFORMED DATA OUTPUT FORMAT FOR THE GENERAL 
SOLID PROPELLANT ROCKET MOTOR PERFORMANCE COMPUTER PROGRAM 


[Transformed impulse] 


time no. 5 


MOTOR NO. 
TEMP GROUP 1 


PA = 
TEMPI 

-♦00 

TEMP2 

TEMP3 

TEMP4 

37A 

70 

44A - 

4bA 

43A- 


1717.24 

1598.79 
1687.36 
1708.90 

1691.79 

1808.10 

1683.38 

1776.63 

1798.20 

1780.72 

1997.65 
1860.04 
1963.08 __ 
1986.00 
1966.69 

.00 

.00 

.00 

.00 

.^on 

74A 


1690.42 

1779.28 

1965.10 

0 

°. 

MEAN 

' ■ 

1662,33 

1771.05 

1956.46 

.00 

.standard dev.. 


. 42*56. 

_ 44.6.9 - 

49.20 _ 

.00 


_JEMP GROUP,. 


77 

1617 ♦ 59 

1704.09 

1884.33 

.00 

50 A 

1725.21 

1816.26 

2006,52 

.00 

. 3BA_ 

_ 1660*79 

1748.16 

1930.84 _ 

.00 

59A 

1707.16 

1796.08 

1982.41 

.00 

MEAN 

1677.69 

1766,15 

1951.03 

.00 

STANDARD DEV. 

48.39 

50*28 

54*53 _ 

.00 

TEMP GROUP .3 

— 

- 



69. 

1636.31 

1723.54 

1905.41 

.00 

64A 

1672.84 

1762.01 

1947.95 

.00 

58A __ 

1746.90 

ia40.02 

2034.19 _ 

*.0Q 

65A 

1676.67 

1765.46 

1950.87 

.00 

40A 

1681.40 

1769.57 

1954.06 

.00 

78 

1633.62 

1718.44 

1896.29 

.00 

6QA_ ____ „ 

1678.60 

1765.76 

1948.50 

.00 

MEAN 

1675.15 

1763.54 

1948.18 ... 

.08 

STANDARD DEV. 

37.51 

39.86 

44.65 

.00 

TOTAL MEAN 

1678.30 

1766.81 

1951.77 

.00 

TOTAi STANDARD DEV. 

39.28 

41.30 

45.56 

. jOQ 

CONFIDENCE ON.. NORMAL 

DISTRIBUTION 




(ONE MIN 

1537.47 

1618.76 

1788.44 

.00 

.SIDED. MAX 

1819,13 

1914.86 

2115.11 

,00 

(TWO MIN 

1522.94 

1603.48 

.1271* 5A 

«_0Q 

SIDED MAX 

1833,66 

1930.14 

2131.97 

.00 
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TABLE VI- INPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER PROGRAM 


Card 

number 

(a) 

Variable 

name 

Format 

Columns 

Description 

Units 

1 

NXNW 

16 

1-6 

Integer designating the number of 
subintervals into which web time 
is divided (up to 5) 



NXNT 

16 

7-12 

Integer designating the number of 
subintervals into which tail- off 
time is divided (up to 4) 

- 


TD 

F12.4 

13-24 

Prefire propellant temperature to 
which the performance data are 
to be transformed 

°F 


PHI 

F12.4 

25-36 

Nozzle cant angle (from the longi- 
tudinal center line of the motor) 

rad 


CT2 

F12.4 

37-48 

Tolerance factor (one- or two- 
sided) 

- 


NSW b 

16 

49-54 

Integer designating the method by 
which TO is to be determined 
(0 if TO is inputed, 1 if TO is to 
be determined by the program as 
a function of chamber pressure) 


1A 

PSIA b 

F10. 0 

1-10 

Value of chamber pressure at 
which TO is determined (card 
not included if NSW = 0) 

psia 

2 

PCTW1 

F10. 4 

1-10 

Upper limit of first web- time 
subinterval 

percent 


NW1 

16 

11-16 

Integer designating the number of 
equally divided per cent -web -times 
in the first subinterval that time, 
chamber pressure, and thrust 
transformations are to be per- 
formed 



PCTW2 

F10. 4 

17-26 

Continue for each subinterval 
up to NXNW (last subinterval 
must equal 100 percent web 
time) 

percent 


a Cards 1 through 10 apply for all motors and are read in only once; however, cards 
NPTS 

11 through 13 + — ~ — are required for each motor and must be entered in sequence, 
h " 

D Card 1A can be omitted if TO is an input value (NSW = 0). 
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TABLE VL - INPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Continued 


Card 

number 

(a) 

Variable 

name 

Format 

Columns 

Description 

Units 

2 

NW2 

16 

27-32 

Continue for each subinterval 
up to NXNW (last subinterval 
must equal 100 percent web 
time) 



PCTW3 

F10. 4 

33-42 


percent 


NW3 

16 

43-48 


- 


PCTW4 

F10. 4 

49-58 


percent 


NW4 

16 

59-64 


- 


PCTW5 

F10. 4 

65-74 


percent 


NW5 

16 

75-80 


- 

3 

PCTT1 

F10. 4 

1-10 

Upper limit of first tail -off time 
subinterval 

percent 


NTl 

16 

11-16 

Integer designating the number of 
equally divided percent tail -off 
times in the first subinterval 
that time, chamber pressure, and 
thrust transformations are to be 
performed 



PCTT2 

F10. 4 

17-26 

Continue for each subinterval 
up to NXNT (last subinterval 
must equal 100 percent tail-off 
time) 

percent 


NT2 

16 

27-32 


- 


PCTT3 

F10. 4 

33-42 


percent 


NT3 

16 

43-48 


- 


a Cards 1 through 10 apply for all motors and are read in only once; however, cards 11 
NPTS 

through 13 + 9 - are required for each motor and must be entered in sequence. 
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TABLE VI. - INPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 


PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Continued 


Card 

number 

(a) 

Variable 

name 

Format 

Columns 

Description 

Units 

3 

PCTT4 

F10. 4 

49-58 

Continue for each subinterval up 
to NXNT (last subinterval must 
equal 100 percent tail -off time) 

percent 


NT 4 

16 

59-64 


- 

4 

NA 

112 

1-12 

Integer designating the number 
of prefire propellant tempera- 
tures about which the motors are 
grouped (up to 3) 



TGRP1 C 

F12. 4 

13-24 

First prefire propellant tempera- 
ture about which the motors are 
grouped (prefire propellant tem- 
perature group no. 1) 

°F 


TGRP2 C 

F12. 4 

25-36 

Continue for each prefire propel- 
lant temperature about which the 
motors are grouped (up to 3) 

°F 


TGRP3 C 

F12. 4 

37-48 


°F 

5 

NP 

112 

1-12 

Integer designating the number of 
pressure altitudes at which thrust 
is to be reported (up to 3) 

- 


EXPA1 

F12. 4 

13-24 

The first pressure altitude at 
which thrust is to be reported 

psia 


EXPA2 

F12. 4 

25-36 

Continue for each pressure alti- 
tude at which thrust is to be 
reported (up to 3) 

psia 


EXPA3 

F12. 4 

37-48 


psia 

6 

HEADX 

12A6 

1-72 

Heading for the abscissa of all 
graphs (time axis) 

sec 

7 

HEAD1 

12A6 

1-72 

Ordinate heading on plot of cham- 
ber pressure versus time 

psia 


Cards 1 through 10 apply for all motors and are read in only once; however, cards 11 
NPTS 

through 13 + — ~ — are required for each motor and must be entered in sequence. 


C The value of TGRP(3) must be greater than TGRP(2), and TGRP(2) must be 
greater than TGRP(l). 



TABLE VI. - INPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER* PROGRAM - Continued 


Card 

number 

(a),(d) 

Variable 

name 

Format 

Columns 

Description 

Units 

8 

HEAD2 

12A6 

1-72 

Ordinate heading on plot of thrust 
at pressure altitude no. 1 
(EXPA 1) versus time 

lbf 

9 

HEAD3 

12A6 

1-72 

Ordinate heading on plot of thrust 
at pressure altitude no. 2 
(EXPA 2) versus time 

lbf 

10 

HEAD4 

12A6 

1-72 

Ordinate heading on plot of thrust 
at pressure altitude no. 3 
(EXPA 3) versus time 

lbf 

11 

XMOT 

A6 

1-6 

Motor number 

- 


FIRETP 

F12.4 

7-18 

Prefire propellant temperature 

°F 


PAX 

F12.4 

19-30 

Ambient pressure at the time of 
motor firing 

psia 


TB 

F12. 4 

31-42 

Burn-time for the motor 

sec 


TT 

F12. 4 

43-54 

Total time for the motor 

sec 


TO 

F12. 4 

55-66 

Time offset (time that is zero 
percent web time) 

sec 


PCP 

F12.4 

67-78 

Prefire chamber pressure (used 
to convert input chamber pres- 
sures from psig to psia) 

psia 


NDMOT 

11 

80 

Integer designating that the last 
motor to be processed has been 
reached (1 for last motor, 0 for 
all other motors) 


12 

AEl 

E12.5 

1-12 

Exit area of nozzle no. 1 

' in. 2 


AE2 

E12. 5 

13-24 

Continue for each nozzle (up to 4) 

. 2 
in. 


a Cards 1 through 10 apply for all motors and are read in only once; however, cards 11 
through 13 + ~ are required for each motor and must be entered in sequence. 

d If fewer than three ambient pressures are used, cards 9 and 10 must be present, even 
as blank cards. 



TABLE VI. - INPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Concluded 


Card 

number 

(a) 

Variable 

name 

Format 

Columns 

Description 

Units 

12 

AE3 

^12. 5 

25-36 

Continue for each nozzle (up to 4) 

. 2 
in. 


AE4 

E12. 5 

37-48 


. 2 
in. 

13 

NPTS 

16 

1-6 

Integer designating number of 
points to be read in to define 
motor performance versus time 
(up to 200) 

' 


ICON 

16 

7-12 

Integer designating the number of 
input points in piecewise curve 
fit (6 is a good choice) 



IPRNT 

16 

13-18 

Integer designating that the curve 
fit coefficients are to be printed 
out (0 for no coefficients, 1 for 
all coefficients) 


14 

TREP1 

E16. 8 

1-16 

First time input in the perform- 
ance input array 

sec 


PCI 

E12. 5 

17-28 

Chamber pressure that corre- 
sponds to the time input in the 
performance input array 

psig 


FI 

E12. 5 

29-40 

Thrust that corresponds to time 
input in the performance input 
array 

lbf 


TREP2 

E16. 5 

41-56 

Continue for each chamber pres- 
sure and thrust at the proper 
time (up to NPTS) at 2 times, 
pressures and thrusts inputs per 
card 

sec 


PC2 

E12. 5 

57-68 


psig 


F2 

E12.5 

69-80 


lbf 


a Cards 1 through 10 apply for all motors and are read in only once; however, cards 11 
NPTS 

through 13 + 0 are required for each motor and must be entered in sequence. 
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TABLE VH. - TRANSFORMED TIME AND CHAMBER- PRESSURE STATISTICAL 


PCI * 

WEB 
TIME 
.00 
_ *81 
1.61 
2.42 
3.23 
4.03 

4.64 

5.65 
6.45 

SI .26 
8.06 
_£U_&7 
9.68 
10*48 ■ 

11.29 
12 * 1-0 
12.90 
13-.71 

14.52 

15.32 

16.13 

16.94 

17.74 
16.55 
19.35' 
20.16 
20.97 
21.77 
22.58 
23.39 
24.19 
25.00 
25.81 
26.61 
27.42 
28.23 
29.03 
29.84 
30.65 
31. 45 
32.26 
33.06 
33.87 
34.68 
35.46 

36.29 
37.10 
37.90 
36.71 

39.52 

40.32 

41.13 

41.94 

42.74 


DATA OUTPUT FORMAT FOR THE SOLID PROPELLANT ROCKET 
MOTOR PERFORMANCE VERSUS TIME COMPUTER PROGRAM 


TRANSFORMED TIMES TRANSFORMED CHAMBER PRESSURE 



MEANS 

WITH TWO 

sided tolerance limits 


MEAN 

MIN. 

MAX. 

MEAN 

MIN. 

MAX. 

.0000 

-.0000 

.0000 

100.0727 

96.2117 

103.9337 

.0050 

.0045. 

.0054 

186.2759 

47*4130 

325.1388 

.0099 

.0091 

.0108 

288.4526 

76.4223 

500.4828 

.0149 

.0136 

.0162 

400.0431 

154.3467 

645.7395 

.0199 

.0182 

.0215 

514.5706 

243.7581 

785.3831 

.0248 

.0227 

.0269 . 

624*0311 

338*039.6.. 

910*0225 

.0298 

.0273 

.0323 

731.4412 

443.1322 

1019.7502 

.0348 

.0318 

.0377 . 

834.9045 

. 548*8521 

1120.9568 

.0397 

.0364 

.0431 

928.8627 

649.3971 

1208.3282 

.0447 

.0409 

.0485 

1011.5578 

735.0849 

1288.0307 

.0497 

.0455 

.0533 

1084.9299 

807.9524 

1361.9073 

.0546 

.0500 

-0592 

1150.4757 

871.0561 

1429.8953 

.0596 

.0545 

.0646 

1208.8554 

926.5884 

1491.1225 

.0645 

.0591 

.0700- 

1281*3881 

974.8762 

— 4548*8998- 

.0695 

.0636 

.0754 

1308.2346 

1017.1144 

1599.3549 

.0 745 

.0682 

.0808 

1350.0422 

1056.0642 

1644.0203 

.0794 

.0727 

.0862 

1387.5762 

1091.8444 

1683.3080 

.0844 

.0773 

.0915 

1421.3233 

1125.666 2 

1716.9605 

.0894 

.0818 

.0969 

1451.6023 

1157.9969 

1745.2076 

.0943 

.0864 

.1023 

1478.7972 

- 1188*0852- 

1769*4291- 

.0993 

.0909 

.1077 

1503.5798 

1216.8781 

1790.2816 

.1043 

.0955 

.1131 

1526.5542 

1244.7761 

1808.3324 

.1092 

.1000 

.1185 

1546.9323 

1270.6313 

1823.2334 

.1142 

.1045 

.1238 

1564.9614 

1294.1378 

1835.7851 

.1192 

.1091 

.1292 

1581.5182 

1315.8552 

1847.1812 

.1241 

.1136 

.1346 

1596*28 7 9.__ 

._ 1335*8879 — 

1856.6880 

.1291 

.1182 

.1400 

1609.3506 

1354.5360 

1864.1652 

.1341 

.1227 

.1454 

1620.4328 

1371.0981 

1869.7676 

.1390 

.1273 

.1508 

1630.3001 

1386.5419 

1874.0583 

.1440 

.1316 

.1562 

1639.5942 

1401.9310 

1877.2574 

.1490 

.1364 

.1615 

1647.5151 

1414.6645 

1880.3656 

.1539 

.1409 

.1669 

1654.6235 

1425*0915- 

- 1883.1555 

.1589 

.1455 

.1723 

1661.2788 

1437.3437 

1885.2139 

.1638 

.1500 

.1777 

1666.9022 

1447.2263 

1886.5780 

.1688 

.1545 

.1831 

1672.2225 

1456.2725 

1888.1725 

.1738 

.1591 

.1885 

1676.4756 

1464.7989 

1888.1524 

.1787 

.1636 

.1938 

1680.2541 

1471.3614 

1889.1469 

.1837 

.1682 

.1992 

1683.6993 

_ 14.76*.9883^ 

.. 1890.4102 

.1887 

.1727 

.2046 

1686.7656 

1482.4679 

1891.0632 

.1936 

.1773 

.2100 

1689.6254 

1485.9445 

1893.3062 

.1986 

.1818 

.2154 

1691.9939 

1489.3389 

1894.6489 

.2036 

.1864 

.2208 

1694.2330 

-1492.0268 

1896,4392 

.2085 

.1909 

.2261 

1696.0794 

1493.4425 

1898.7162 

.2135 

.1955 

.2315 

1697.3031 

.1495. 5069 _ 

L8.99*_0_993_ 

.2185 

.2000 

.2369 

1698.3392 

1497.4662 

1899.2123 

.2234 

.2045 

.2423 

1698.8275 

1498.6623 

1898.9928 

.2234 

.2091 

.2477 

1699.3327 

1499.6715 

1898.9940 

.2334 

.2136 

.2531 

1699.4336 

1500.4301 

1898.4371 

.2383 

.2182 

.2585 

1699.4032 

1501.1038 

1897.7027 

.2433 

.2227 

.2633 

1699.6042 

1501.03.58 .. 

1 09 8_._1J72.6_ 

.2483 

.2273 

.2692 

1699.4354 

1500.4852 

1898.3857 

.2532 

.2318 

.2746 

1698.9325 

1498.8431 

1899.0219 

.2582 

.2364 

.2300 

1698.1861 

1497.3153 

1899.0568 

.2631 

.2409 

.2854 

1697.0608 

1495.6663 

1898.4554 
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TABLE TO. - TRANSFORMED TIME AND CHAMBER- PRESSURE STATISTICAL 
DATA OUTPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Continued 


43.55 

.2681 

.2455 

.2908 

1695.5937 

1493.6824 

1897.5050 

44.35 

.2731 

.2500 

.2961 

1694.0745 

1491.7364 

1896* 4126 

Hb. ib 

.2780 

• 2545 

.3015 

1692.4855 

1489.7363 

1895.2347 

45.97 

.2830 

.2591 

.3069 

1690.9648 

1487.6661 

1894*2635 

46.77 

.2680 

.2636 

.3123 

1689.2429 

1486.1593 

1892.3265 

47 . 58 

.2929 

.2662 

.3177 

1687.9059 

- 1484.8948 

1890.917a 

48.59 

.2979 

.2727 

.3231 

1686.0139 

1482.4461 

1889.5818 

49. i9 

.3029 

.2773 

.3285 

1683.7669 

1480.0048 

1887* 529a 

50.00 

.3078 

.2618 

.3338 

1681.6899 

1477.6165 

1085.7634 

5<j.81 

.3128 

.2864 

.3392 

1679.0422 

1474.1835 

1883.9009^ 

51 .bl 

.3178 

.2909 

,3446 

1676.3945 

1470.2831 

1882.5059 

52. 42 

.3227 

.2955 

.3500 

1673.4878 

1466.3678 

1880.6078 

53.23 

.3277 

.3000 

.3554 

1670.3387 

1463.3387 

1877.3388 

54.03 

.3327 

.3045 

.3603 

1667.2517 

1460.9975 

1873. sasa 

54.84 

.3376 

.3091 

.3661 

1663.6708 

1457.1346 

1870.2070 

55.65 

.3426 

.3136 

.3715 

1659.7595 

1453.2724 

1866*2467- 

56.45 

.3476 

.3182 

.3769 

1655,1865 

1449.5914 

1860.7816 

57.26 

.3523 

.3227 

.3323 

1650.3546 

1444.5683 

1856. 1418 

5b, 06 

.3575 

.3273 

.3877 

1645.4045 

1439.3457 

1851.4635 

58.87 

.3624 

.3318 

.3931 

1639.4879 

1433.0053 

1845.9706 

59.68 

• 3o74 

• 33b4 

.3985 

1633.2585 

1426.7380 

1839.7791 

60.48 

.3724 

. 3409 

.4038 

1626.4598 

1420.3194 

1832.6002 

61.29 

.3773 

.3455 

.4092 

1619.3300 

1412.7360 

1825.9241 

62.10 

.3823 

.3500 

.4146 

1612.2121 

1405.1977 

1819.2267 

&2.90 

.3873 

• 3545 

.4200 

1604.6806 

1399.5986 

1809.7626 

63.71 

.3922 

.3591 

.4254 

1596.7055 

1392.7685 

1800.6426 

64.52 

.3972 

.3636 

.4308 

1588.9195 

1386.0671 

1791.7719 

65.32 

.4022 

.3b02 

.4361 

1581.3076 

1377.7398 

1784.8756 

66.13 

.4071 

.3727 

.4415 

1573.1360 

1369.0930 

1777.1790 

66.94 

.4121 

.3773 

.4469 

1564.8889 

1361.5914- 

1768*1865 

67.74 

.4171 

.3818 

.4523 

1556.3409 

1351.5890 

1761.0929 

68. o5 

.4220 

.3864 

.4577 

1547.9214 

1344.4722 

1751*3706 

69.35 

.4270 

.3909 

.4631 

1538.9636 

1336.8781 

1741.0492 

70.16 

.4320 

.3955 

.4685 

1529.6991 

1328.6521 

1730*7462 

70,97 

.4369 

.4000 

.4733 

1520.4931 

1320.7040 

1720.2822 

71.77 

.4419 

.4046 

.4792 

1511.7381 

1313.5326 

1709*9436 

72.58 

.4469 

.4091 

.4846 

1503.2595 

1306.7792 

1699.7398 

73.39 

,4518 

,4136 

.4900 

1494.4600 

1299.1186 

1689.8014 

74.19 

.4568 

.4182 

.4954 

1485.9684 

1292.4371 

1679.4997 

7b. QQ 

,4617 

.4227 

.5008 

1477.3739 

1285.7598 

1668.9881 

75.81 

• 4o67 

.4273 

.5061 

1467.8084 

1277.2995 

1658.3173 

76.61 

.4717 

.4318 

.5115 

1457.9984 

126B • 0385 

1642.9583- 

77.42 

.4766 

.4364 

.5169 

1447.9451 

1258*8627 

1637.0274 

78,23 

,4816 

,4409 

.5223 

1438.8756 

1252.0673 

1625.6840 

79.03 

♦ 4866 

.4455 

.5277 

1430.2524 

1245.2208 

1615,2841 

79.84 

,4915 

,4500 

.5331 

1421.2985 

1237.7059 

16Q4.8912 

80.65 

.4965 

.4546 

.5384 

1411.9651 

1229.8191 

1594.1112 

61.45 

.5015 

.4591. 

.5438 . 

1402,8683.. 

. 1222.0734 

1 583.6631 

82.26 

.5064 

.4636 

.5492 

1393.4044 

1213.5523 

1573.2565 

83.06 

.5114 

.4682 

.5546 

1384.1218 

1205.4579 

1 562.7856 

83.87 

.5164 

.4727 

.5600 

1374.9182 

1197.3819 

1552.4545 

84.68 

.5213 

.4773 

.5654 

1365.6267 

1189.2108 

1542.0426 

85.48 

.5263 

.4818 

.5708 

1356.5140 

1181.2270 

1531.8011 

86.29 

.5313 

.4864 

.5761 

1347.228a 

1 172.4927 

1521.9639 

87.10 

.5362 

.4909 

.5815 

1337.7622 

1163.4899 

1512.0345 

87.90 

,5412 

.4955 

.5869 

1328.4356 

1154.5029 

1502.3684 

88.71 

.5462 

.5000 

.5923 

1318.8686 

1145.1311 

1492.6062 

89.52 

.5511 

.5046 

.5977 

1308.9783 

1136.0403 

1481*9163. 
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TABLE VC. - TRANSFORMED TIME AND CHAMBER- PRESSURE STATISTICAL 
DATA OUTPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 
PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Continued 


90,32 

.5561 

.5091 

.6031 

1299.0389 

1126.5444 

1471.5334 

9 JU 13 

.5610 

.5136 

.6084 

1288.9000 

1116.8996 

1460.9004 

91,94 

• 5660 

.5182 

.6138 

1278.3531 

1106.6177 

1450.0885 

9.2,74 

.5710 

.5227 

.6192 

1267.3961 

1095.2660 

1439.5262 

93.55 

.5759 

.5273 

.6246 

1256.3310 

1064.2116 

1428.4503 

_ 94,35. 

. ,5809 

.5318 

.6300 

_1245.* H 392 

1 jQ 72-*4 B 62. 

1417.9923 

95.16 

.5859 

.5364 

.6354 

1234.1475 

1060.2894 

1408.0054 

. 95.97 

.5908 

.5409 

.6408 

1223.1327 

1048.0887 

1398,1768 

96.77 

.5958 

.5455 

.6461 

1212 • 0997 

1035.5067 

1388.6929 

- 97.58 

.6008 

.5500 

.6515 

1200.4740 

1021.8971 

1379.0508 

98.39 

.6057 

.5546 

.6569 

1185.8283 

1003.6930 

1367.9636 

. 99.19 

.6107 

.5591 

.6623 

1166.4571 

981,6796. 

_1355_._2.347 

100.00 

.6157 

.5636 

.6677 

1145.7223 

952.4611 

1338.9837 


-PCI* , 

tailoff 

-JLLME 


1.33 

.6206 

.5686 

.6725 

1118.7697 

924.0121 

1313.5273 

2.67 

.6254 

.5736 

.6773 

1086.7690 

. .892,7339 . 

1280.8042 

4.00 

.6303 

.5785 

.6321 

1050.3857 

856.7386 

1244.0328 

__ 5,33 

.6352 

.5835 

.6370 

1011.7241 

821.655 2 

1201.7930 

6.67 

.6401 

.5883 

.6919 

971.5473 

786.4144 

1156.6802 

8,00 

.6450 

.5932 

.6968 

931.6782 

752.8625 

1110.4938 

9.33 

.6499 

.5980 

.7017 

892.6076 

721.8643 

1063.3509 

10.67 

.6548 

.6029 

.7067 

855,3248 

. 693,0366 _ 

1017,6131 

12.00 

.6597 

.6077 

.7117 

819.4540 

665.4108 

973.4971 

13,33 

.6646 

.6124 

.7167 

785.5828 

639.5456 

931.6200 

14.67 

.6694 

.6171 

.7217 

753.2597 

615.1995 

891.3198 

. 16,00 

.6743 

.6219 

.7268 

722.6703 

592.2510 

853.0896 

17.33 

.6792 

.6265 

.7319 

693.6676 

570.4025 

816.9326 

18.67 

.6841 

.6312 

.7370 

666.1639 

549,3637... 

.702.9741. 

20.00 

.6890 

.6359 

.7422 

640.3914 

529.4356 

751.3472 

21.33 

.6939 

.6405 

.7473 

615.8021 

510.9543 

720.6499 

22.67 

.6988 

.6451 

.7525 

592.5075 

493.9920 

691.0230 

_ . 24 * 0 0 

.7037 

.6496 . 

.7577 

570.0886 

477.3924 

662.7849 

25.33 

.7066 

.6542 

• 7 t >29 

548.8185 

461.4724 

636.1645 

26.67 

.7134 

.6587 

.7682 

528.5498 

446,140.9 

610,9587. 

28.00 

.7183 

.6632 

.7735 

509.2877 

431.4305 

587.1449 

29.33 

.7232 

.6677 

.7788 

491.0974 

417.5962 

564.5985 

30.67 

.7281 

.6721 

.7841 

473.4347 

403.6292 

543.2402 

32.00 

.7330 

• 6766 

.7894 

456.6112 

390,2160 

523.0063 

33.33 

.7379 

.6810 

.7948 

440.9169 

377.3766 

504.4572 

34.67 

.7428 

.6854 

.8002 

425.7773 

365.3657 

.. _ 486,-1889 

36.00 

.7477 

.6898 

.3056 

411.5136 

354.3595 

468.6677 

37.33 

.7526 

.6942 

.8110 

398,0980 

343.7363 

452.4596 

38.67 

.7575 

.6985 

,8164 

384.9482 

333.3091 

436.5872 

40.00 

.7623 

.7028 

.8219 

372.6211 

323.3540 

421.8882 

41.33 

.7672 

.7071 

.8273 

360.9189 

314.0585 

407.7792 

42.67 

.7721 

.7114 

.8328 

349.3216 

. 304,5711 . 

_._394,0722 

44.00 

.7770 

.7157 

.8383 

338.4454 

295.3758 

381.5149 

... 45 • 33 

.7819 

.7200 

,8438 

327.8011 

286.6873 

368.9149 

46.67 

. 78 t >8 

.7242 

.8494 

317.3519 

278.4996 

356.2041 

48.00 

.7917 

.7285 

.3549 

307.4221 

269.9771 

344,8672 

49.33 

.7966 

.7327 

.8605 

297.9488 

261.7621 

334.1354 

50.67 

.8015 

.7369 

• 8660 

288.7212 

254.5062 

.322.936 2 

52.00 

.8063 

.7411 

.8716 

279.9151 

247.4768 

312.3534 

53.33 

.8112 

.7453 

.8772 

271.7372 

240.9243 

302,5501 

54 . 67 

.8161 

.7494 

. 8328 

263.6080 

234.3644 

292.8517 

56.00 

.8210 

.7536 

.888 4 

255.9788 

228.0174 

283.9402 


TABLE VII. - TRANSFORMED TIME AND CHAMBER- PRESSURE STATISTICAL 
DATA OUTPUT FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR 


PERFORMANCE VERSUS TIME COMPUTER PROGRAM - Concluded 


57.33 

.8259 

.7577 

.8941 

248.3934 

221.6398 

275.1469 

50.67 

.8308 

.7619 

.8997 

241.0648 

215.5795 

266 . 5502 - 

60.00 

.8357 

.7660 

.9054 

234.3061 

209.6224 

258.9898 

61.33 

.8406 

.7701 

.9110 

227.6651 

203.7331 

251 . 5970 - 

62.67 

.8455 

.7742 

.9167 

221.2060 

198.5070 

243.9050 

64.00 

.8503 

.7783 

.9224 

214.8457 

193.2729 

236.4185 

65 • 33 

♦ 8552 

.7824 

.9281 

208.5310 

187.7487 

229.3132 

66.67 

.8601 

.7865 

.9338 

202.1990 

182.8354 

221.5626 

68.00 

.8650 

.7905 

.9395 

196.4787 

178.2372 

214.7202 

69.33 

.8699 

.7946 

.9452 

191.1092 

173.6381 

208.5804 

70 . 67 

.8748 

.7986 

.9510 

105.6251 

169.1642 

202.0860 

72.00 

.8797 

.8027 

. 9567 

180 . 279 a 

164.7641 

195.7957 

73.33 

.8846 

.8067 

.9625 

175.4276 

160.8767 

189.9785 

74.67 

.8895 

.8107 

.9682 

170.6171 

156.4750 

184 . 759 i - 

76.00 

.8944 

.8147 

.9740 

165.6604 

152.0464 

179.2745 

77.33 

.8992 

♦ 8187 

.9797 

160.6167 

147.5117 

173.7218 

78.67 

.9041 

.8227 

.9855 

156.2894 

143.6445 

168.9342 

ao.oo 

.9090 

.8267 

.9913 

152.1756 

139 . 9957 - 

164.3556 

61. 33 

.9139 

.8307 

.9971 

147.9813 

136.9321 

159.0305 

82.67 

.918 a 

.8347 

1.0029 

144.0087 

133.8443 

154.1731 

84.00 

.9237 

.8387 

1.0087 

139.9193 

130.7487 

149.0900 

85.33 

.9286 

.3427 

1.0145 

136.0162 

128.7390 

143.2933 

86.67 

.9335 

.8466 

1.0203 

132.1325 

125.1921 

139.0730 

86.00 

.9384 

.8506 

1.0261 

128.4227 

121.7876 

135.0579 

89.33 

.9432 

.8545 

1.0320 

124.8258 

118.3964 

131.2551 

90.67 

.9481 

.8585 

1.0378 

121.6686 

115.0492 

128.2881 

92.00 

.9530 

.8624 

1.0436 

118.5244 

111.5678 

125.4810 

93.33 

.9579 

.8664 

1.0495 

115.4463 

108.7479 

122.1447 

94.67 

.9628 

.8703 

1 . 0553 

112.2210 

106.0667 

118.3753 

96.00 

.9677 

.8742 

1.0611 

108.9699 

103.3730 

114.5669 

97.33 

.9726 

.8782 

1.0670 

106.0028 

100.9064 

111.0992 

98.67 

.9775 

.8821 

1.0723 

103.1398 

98.1133 

108.1662 

100.00 

• 9824 

.8860 

1.0787 

100.3185 

95.4782 

105.1589 


lj-7 



TABLE m - TRANSFORMED TIME-AND- THRUST STATISTICAL DATA OUTPUT 
FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
VERSUS TIME COMPUTER PROGRAM 


transformed times transformed thrust at pa = 


PCI. 


MEANS 

WITH TWO 

SIDED TOLERANCE LIMITS 


WEB 

MEAN 

MIN, 

MAX. 

MEAN 

MIN. 

MAX. 

TIME 







.00 

.0000 

-.0000 

.0000 

29.9100 

-85.4260 

145.2459 

.81 

* 0050 

* Q045 

• 0054 

128,6946 

-81,6653 

. 339,0546 

1.61 

.0099 

.0091 

.0108 

281.9090 

1.7245 

562.0935 

2.42 

,0149 

.0136 

,0162 

471,5607 

140,4345 

802,6369 

3.23 

,0199 

.0162 

.0215 

694.4515 

321.5112 

1067.3919 

4.03 

.0248 

.0227 

. U2fa9 

924.5868 

547.2773 

1301.8962 

4.84 

.0298 

.0273 

.0323 

1137.4147 

732.9413 

1541.8882 

5.65 

.0348 

.0316 

.0377 

1333.2352 

896.6773 

1769.7933 

6.45 

.0397 

.0364 

.0431 

1509.5691 

1058.4309 

1960.7075 

7,26 

.0447 

.0409 

.0485 

1663. 7580 

1204.7103 

2122. 8Q56 

8.06 

.0497 

.0455 . 

.0538 

1801.8282 

1336.0550 

2267.6016 

8.67 

,0546 

.0500 

.0592 

1920.1355 

1451.1768 

2389.0942 

9.68 

.0596 

.0545 

.0646 

2020.9526 

1542.3207 

2499.5845 

10.48 

.0645 

.0591 

.0700 

2111.9834 

1620.7030 

2603,2638 

11.29 

.0695 

• 063b 

.0754 

2194.1727 

1700.4281 

2687.9172 

12.10 

.0745. 

.0682 

.0808 

2267.3524 

1766.0181 

2768,6867 

12.90 

.0794 

.0727 

.0862 

2332.6759 

1817.9741 

2847.3777 

13,71 

.0844 

.0773 

.0915 

2391.3051 

1873.9083 

£908,7020 

14.52 

.0894 

.0818 

.0969 

2445.1623 

1932.0522 

2958.2724 

15.32 

.0943 

.0864 

.1023 

2495.1170 

1983.4714 

3006,7626 

16,13 

.0993 

.0909 

.1077 

2539.3426 

2032.2972 

3046.3880 

16,94 

. .1043 

.0955 

.1131 

2577.5675 

2078.8729 

3076.2620 

17.74 

.1092 

,1000 

.1185 

2613.2377 

2122.4005 

3104.0750 

18,35 _ 

.1142 

.1045 

,1238 

2645.6056 

2165,5149 

3125.6963 

19:35 

.1192 

.1091 

.1292 

2673.7763 

2205.3074 

3142.2453 

20.16 

.1241 

.1136 

.1346 

2699. S305 

2236.7407.. 

3162.9203 

20.97 

.1291 

.1182 

.1400 

2722.2712 

2265.7028 

3178.8398 

22,77 

.1341 

.1227 

.1454 

2740.8897 

2299.9138 

3181.8657 

22.58 

.1390 

.1273 

. 1508 

2758.4704 

2328.6655 

3188.2754 

23.39 

• 144Q 

.1318 

.1562 

2774,7448 

2352.4293 

3197.0602 

24.19 

.1490 

• l3t>4 

.1615 

2788.7643 

2373.2424 

3204.2862 

25.00 

.1539 

.1409 

.1669 

2800.5063 

2394.6269 

3206.3857 

25.81 

.1589 

.1455 

.1723 

2811.7792 

2415.0017 

3208.5566 

26,61 

.1638 . 

.1500 

.1777- 

2821.9678 

2429.0447 

3214.8910 

27.42 

.1688 

.1545 

.1831 

2830.4479 

2444.2039 

3216.6918 

28.21 

.1738 

,1591 

.1885 

2837.2719 

2456.9333 

3217.6106 

29.03 

.1787 

.1636 

.1938 

2842.8841 

2466.1893 

3219.5789 

29.84 

.1837 

.1682 

.1992 

2849.0198 

2475.7554 . 

3222.2841 

30.65 

.1887 

.1727 

.2046 

2854.2417 

2486.5305 

3221.9529 

31,45 

.1936 

,1773 

.2100 

2858,9887 

2494,0985 

3223,6790 

32.26 

.1986 

.1818 

.2154 

2863,4712 

2497.6341 

3229.3004 

33.06 

.2036 

.1864 

.2208 

2866.8565 

2503.7991 

3229.9140 

33.87 

.2085 

.1909 

.2261 

2869.7965 

2512.1626 

3227.4305 

34.68 

,2135 

.1955 

-,2315 

2871 • 5731 

.2518,4520. _ 

.3224.6342 

35.48 

.2185 

.2000 

.2369 

2873.1184 

2522.2716 

3223.9653 

36.29 

.2234 

.2045 

.2423 

2873.8475 

2524.1828 

3223.5122 

37.10 

.2284 

.2091 

.2477 

2874.6448 

2525.6334 

3223.6563 

37.90 

.2334 

.2136 

.2531 

2875.0168 

2526,4722 

3223.5613 

38.71 

.2383 

.2182 

.2585 

2874.4876 

2524.0527 

3224.9225 

39.52 

.2433 

.2227 

.2638 

2874,7769 

2524.4575 

3225.0962 

40.32 

.2483 

.2273 

.2692 

2874.8362 

2524.5221 

3225.1503 

41*13 

,2532 

.2318 

.2746 

2874.4308 

2526.1292 

3222,7324 

41.94 

.2582 

.2364 

.2800 

2873.1526 

2526.2619 

3220.0433 

42,74 

.2631 

.2409 

.2854 

2871.0146 

2525.0269 

3217.0023 



TABLE vm. - TRANSFORMED TIME-AND-THRUST STATISTICAL DATA OUTPUT 
FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
VERSUS TIME COMPUTER PROGRAM - Continued 


43.55 

.2681 

.2455 

.2908 

2868.3248 

2523.1066 

3213.5431 

44.35 

.2731 

.2500 

.2961 

2865.2769 

2520.5601 

3209*9938- 

45.16 

.2780 

.2545 

.3015 

2862.4629 

2517.4233 

3207.5026 

45.97 

.2830 

.2591 

.3069 

2859.7900 

2513.0028 

3206*5773- 

4b. 77 

.2380 

.2636 

.3123 

2856.6862 

2510.6798 

3202.6925 

47.58 

.2929 

.2682 

.3177 

2854.2220 

2508.1865 

3200.257-6 

48.39 

.2979 

.2727 

.3231 

2851.2232 

2504.6908 

3197.7557 

49.19 

.3029 

.2773 

.3285 

2848.1022 

2502.6069 

3193.5976 

50.00 

.3078 

.2818 

.3338 

2844.7198 

2498.3733 

3191.0663 

50.81 

.3128 

.2864 

.3392 

2839.4213 

2494.3550 

3184.4877 

51.61 

.3178 

.2909 

.3446 

2834.2778 

2491.1850 

3177.3707 

5a . 42 

.3227 

.2955 

.3500 

2829.4647 

2487.6159 

3171.3136- 

53.23 

.3277 

.3000 

.3554 

2824.1272 

2462.0980 

3166.1563 

54. 03 

.3327 

.3045 

.3608 

2818.4720 

2475.7414 

3161*2025- 

54.84 

.3376 

.3091 

.3661 

2812.7021 

2470.3959 

3155.0082 

55.65 

.3426 

.3136 

.3715 

2807.0796 

2464.7331 

3149*4260- 

56.45 

.3476 

.3182 

.3769 

2799.8079 

2456.7768 

3142.8389 

57.26 

.3525 

.3227 

.3823 

2791.6342 

2444.7939 

3138.4X45- 

58.06 

.3575 

.3273 

.3877 

2783.1130 

2433.6220 

3132.6041 

58.87 

.3624 

.3318 

.3931 

2773.3599 

2425.3489 

3121 • 3707, 

59.68 

.3674 

.3364 

.3985 

2762.9403 

2413.6415 

3112.2390 

60.48 

.3724 

.3409 

.4038 

2751.2329 

2399.6761 

3102.7897 

61.29 

.3773 

.3455 

.4092 

2740.2078 

2391.2279 

3089.1877 

62.10 

.3823 

.3500 

.4146 

2728.8394 

2382.7666 

3074.9121 

62.90 

.3873 

.3545 

.4200 

2715.7695 

2372.2278 

3059.3112 

63.71 

.3922 

.3591 

.4254 

2702.5398 

2357.6389 

3047.4407- 

64.52 

.3972 

.3636 

.4308 

2688.8404 

2343.0470 

3034.6339 

65.32 

.4022 

.3682 

.4361 

2675.3877 

2330.1481 

3020.6273 

66.13 

.4071 

.3727 

.4415 

2662.1713 

2319.3293 

3005.0132 

66.94 

.4121 

.3773 

• 44b9 

2649.3055 

2310.287a 

- 2988*3231- 

67.74 

.4171 

.3818 

.4523 

2637.0740 

2309.5849 

2964.5630 

68.55 

.4220 

.3864 

.4577 

2623.4088 

2303.8965 

2942.9211 

69.35 

.4270 

.3909 

.4631 

2608.4831 

2289.6410 

2927.3251 

70.16 

.4320 

.3955 

.4685 

2593.8872 

2280.3134 

2907.4610 

70.97 

.4369 

.4000 

.4738 

2577.2821 

2259.7080 

2894.8563 

71.77 

.4419 

.4046 

.4792 

2561.1299 

2240.5366- 

2881 *-72 32 

72.58 

.4469 

.4091 

.4846 

2545.1835 

2220.8183 

2869.5489 

73.39 

.4518 

.4136 

.4900 

2529.2288 

2200.8713 

2857.5862 

74.19 

.4568 

.4182 

.4954 

2514.2934 

2190.9509 

2837.6359 

75.00 

.4617 

.4227 

.5008 

2499.9293 

2180.8189 

2819.0397 

75.81 

.4667 

.4273 

.5061 

2483.7794 

2163.4445 

2804.1142 

76.61 

.4717 

.4318 

.5115 

2467.0346 

2144.9338,, 

2789..-1355 

77.42 

.4766 

.4364 

.5169 

2450.3193 

2127.9415 

2772,6971 

78.23 

.4816 

.4409 

.5223 

2435.9662 

2118.0287 

2753*4837 

79.03 

.4866 

.4455 

.5277 

2421.6738 

2106.7662 

2736.5813 

79.84 

.4915 

.4500 

,5331 

2406.6572 

2094.8505 

2718*4648 

80.65 

.4965 

.4546 

.5384 

2391.4948 

2084.5418 

2698.4479 

81.45 

.5015 

.4591 

.5438 

2376.3506 

2072.0121 

— 26BQ .6892 

82.26 

.5064 

.4636 

.5492 

2359.9619 

2058.2972 

2661.6267 

83.06 

.5114 

.4682 

.5546 

2343.4240 

2046.6704 

2640*1778 

83.87 

.5164 

.4727 

.5600 

2328.0088 

2034.0621 

2621.9555 

84.68 

.5213 

.4773 

.5654 

2312.4452 

2020.0016 

?604.ftflflfl 

85.48 

.5263 

.4818 

.5708 

2296.8338 

2006.6184 

2587.0491 

86.29 

.5313 

.4864 

■ 5761. 

2280.9606, 

--1994.8274 

2567.0937 

87.10 

.5362 

.4909 

.5815 

2265.0679 

1981.7576 

2548.3783 

87.90 

.5412 

.4955 

.5869 

2249.1552 

1965.7215 

2532.5889 

88.71 

.5462 

.5000 

.5923 

2232.3401 

1949.2994 

2515.3807 

89.52 

.5511 

.5046 

.5977 

2215.3879, 

1934.6026 

2496.1732 



TABLE Vm. - TRANSFORMED TIME-AND- THRUST STATISTICAL DATA OUTPUT 
FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
VERSUS TIME COMPUTER PROGRAM - Continued 


90.32 

.5561 

.5091 

.6031 

2198.1009 

1919.2804 

2476.9214 

91.13 

.5610 

.5136 

.6084 

2181. 0201 

. 1901-8137 

2460*2263 

91.94 

.5660 

.5182 

.6138 

2163.9528 

1883.1708 

2444.7348 

92.74 

.5710 

.5227 

.6192 

2145,8561 

1865.0740 

2426,6382 

93.55 

.5759 

.5273 

.6246 

2127,6139 

1847.1024 

2408.1255 

94.35 

.5809 

.5318 

.6300 

2108.6294 

1827.6688 

2389,5900 

95.16 

.5859 

.5364 

.6354 

2089.5475 

1807.7857 

2371.3093 

95.97 

.5908 

.5409 

.6408 

2070,9862 

1787.6872 

2354,2851 

96.77 

.5958 

.5455 

.6461 

2052.5068 

1769.6109 

2335.4027 

97.58 

.6008 

.5500 

.6515 

2033.8807 

1750.0885 

2317,6729 

98.39 

.6057 

.5546 

.6569 

2010.3788 

1718.7040 

2302.0537 

99.19 

.6107 

.5591 

.6623 

1983.3887 

1682.3773 

2284,3999 

100.00 

PCI . 

tailoff 

TIME 

.6157 

.5636 

.6677 

1947.5539 

1638.6585 

2256.4492 

1.33 

.6206 

.5686 

.6725 

1904.0982 

1589.7141 

2218.4823 

2.67 

.6254 

.5736 

.6773 

1852.2381 

1536.7176 

2167.7586 

4.00 

.6303 

.5785 

.6821 

1793.0064 

1474.8251 

2111.1877 

5.33 

.6352 

.5835 

.6870 

1727,7396 

1415.1121 

2040*3671 

6.67 

.6401 

.5883 

.6919 

1659.3158 

1359.2429 

1959.3887 

8.00 

.6450 

.5932 

.6968 

1591.0564 

1300.3271 

1861.7856 

9.33 

.6499 

.5980 

.7017 

1523.3241 

1242.0927 

1804.5556 

10.67 

.6548 

.6029 

.7067 

1459.6777 

1187.9408 

1731.4146 

12.00 

.6597 

.6077 

.7117 

1399.1327 

1141,0040 

1657.2614 

13.33 

.6646 

.6124 

.7167 

1341.8144 

1097.6074 

1586,0214 

14.67 

.6694 

.6171 

.7217 

1288.3407 

1049.1997 

1527.4816 

16.00 

.6743 

.6219 

.7268 

1236.0442 

10 05 . 2416 _ 

1466.8468 

17.33 

.6792 

.6265 

.7319 

1186.3722 

964.8985 

1407.8460 

18.67 

.6841 

.6312 

.7370 

1139.1027 

927.2142 

. 1350.9913 

20.00 

.6890 

.6359 

.7422 

1095.0792 

892.5298 

1297.6286 

21.33 

.6939 

.6405 

.7473 

1053.0400 

856,7447 .. 

1249.3352 

22.67 

.6988 

.6451 

.7525 

1013.0901 

824.4401 

1201.7401 

24.00 

.7037 

.6496 

.7577 

974.8540 

796.1760 _ 

1153.5320 

25.33 

.7086 

.6542 

.7629 

938.6829 

768.1775 

1109.1883 

26.67 

♦ 7134 

.6587 

.7682 

904.5105 

739.2651 

1069.7559 

28.00 

.7183 

.6632 

.7735 

872.0610 

711.4493 

1032.6727 

29.33 

.7232 

.6677 

.7788 

840,9072 

686.0063 

. 995.8081 

30.67 

.7281 

.6721 

.7841 

810.0386 

659.9553 

960.1218 

32.00 

.7330 

.6766 

.7694 

780.9014 

634.7568 

_ 927,0461 

33.33 

.7379 

.6810 

.7948 

753.9982 

611.2459 

896.7505 

34.67 

.7428 

.6854 

.8002 

728.2765 

589.3771 

_ 867*1760 

36.00 

.7477 

.6898 

.8056 

704.2560 

570.6066 

837.9053 

37.33 

.7526 

.6942 

.8110 

681,1219 

. 550.2586 

. __ 812,0052 

38.67 

.7575 

.6985 

.8164 

658.7540 

528.6794 

788.8285 

40.00 

.7623 

.7028 

.8219 

637.8936 

510.2441 

765,5431 

41.33 

.7672 

.7071 

.8273 

617.7568 

493.5347 

741.9790 

42.67 

♦ 7721 

.7114 

.8328 

598.1981 

_ 476.5672 

._ 7i9 . 889Q „ 

44.00 

.7770 

.7157 

.8383 

580.0054 

460.5784 

699.4325 

45.33 

.7819 

.7200 

.8436 

561.9808 

444.6524 _ 

„ 679.3093 

46 . 67 

.7868 

.7242 

.8494 

544,4940 

430.1607 

658.8073 

48.00 

.7917 

.7285 

,8549 

527,5411 

416.6257 

638,4565 

49.33 

.7966 

.7327 

,8605 

511.1165 

402.8734 

619.3596 

50.67 

.8015 

.7369 

,8660 

495.4716 

387.9626 

. _ 602 - 9B07 

52.00 

.6063 

.7411 

.8716 

480.7791 

375.7061 

585,8521 

53.33 

.8112 

.7453 

.8772 

465.7836 

363.2112 

568.3561 

54.67 

.8161 

.7494 

.8828 

451.2083 

348.3974 

554.0193 

56.00 

.8210 

.7536 

,8884 

438.0007 

336.8681 

539.1334 


50 



TABLE vm. - TRANSFORMED TIME -AND -THRUST STATISTICAL DATA OUTPUT 
FORMAT FOR THE SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 


VERSUS TIME COMPUTER PROGRAM - Concluded 


57.33 

.8259 

.7577 

.8941 

424.4465 

324.1052 

524.7877 

58.67 

• 8308 

.7619 

.8997 

412.4682 

314.0708 

510.8655 

60.00 

.8357 

.7660 

.9054 

401.7756 

307.3071 

496.2441 

61.33 

• 8406 

.7701 

.9110 

390.5586 

298.3997 

482.7176 

62.67 

.8455 

.7742 

.9167 

380.7095 

290.9875 

470.4315 

64.00 

.8503 

.7783 

.9224 

370.5793 

279.2780 

461.8806 

65.33 

.8552 

.7824 

.9281 

359.7040 

265.4061 

454.0018 

66.67 

• 8601 

.7865 

.9338 

348.3283 

253.0975 

443.5590 

68.00 

.8650 

.7905 

.9395 

337.3324 

242.0264 

432.6385 

69.33 

.8699 

.7946 

.9452 

328.3641 

234.3927 

422.3355 

70.67 

.8748 

.7986 

.9510 

319.0463 

227.0833 

411.0094 

72.00 

.8797 

.8027 

.9567 

309.4565 

219.1186 

399.7944 

73.33 

.8846 

.8067 

.9625 

300.9247 

211.2617 

390.5877 

74.67 

.8895 

.8107 

.9682 

292.8332 

204.5929 

381.0734 

76.00 

.8944 

.8147 

.9740 

284.5126 

196.5410 

372.4843 

77.33 

.8992 

.8187 

.9797 

276.4818 

168.9881 

363.9755 

78.67 

.9041 

.8227 

.9855 

269.1266 

182.8596 

355.3936 

80.00 

.9090 

.8267 

.9913 

262.2100 

177.3512 

347.0687 

81.33 

.9139 

.8307 

.9971 

255.5590 

173.4801 

337.6380 

82.67 

.9188 

.8347 

1.0029 

248.8581 

166.6553 

331.0609 

84.00 

.9237 

.8387 

1.0087 

241.7264 

159.6450 

323.8078 

85.33 

.9286 

.8427 

1.0145 

235.0206 

155.6245 

314.4167 

86.67 

.9335 

.8466 

1.0203 

228.3924 

150.5706 

306.2142 

88.00 

.9384 

.8506 

1.0261 

222.2223 

146.2521 

298.1924 

89.33 

.9432 

.8545 

1.0320 

216.5944 

141.0506 

292.1381 

90.67 

.9481 

.8585 

1.0378 

211.6711 

136.8468 

286.4954 

92.00 

.9530 

.8624 

1.0436 

206.2429 

131.1978 

281 .2881 

93.33 

.9579 

.8664 

1.0495 

200.4277 

123.9730 

276.8824 

94.67 

.9628 

.8703 

1.0553 

195.1405 

117.3585 

272.9225 

96.00 

.9677 

.8742 

1.0611 

190.1864 

113.0641 

267.3087 

97.33 

.9726 

.8782 

1.0670 

185.5576 

110.6363 

260.4790 

98.67 

.9775 

.8821 

1.0728 

180.3647 

108.3940 

252.3354 

100.00 

.9824 

.8860 

1.0787 

175.2670 

105.3755 

245.1585 
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(a) Longitudinal cross section. 


NASA-S-66-2025 FEB 24 



(b) Typical internal-burning six-point star propellant 
grain configuration. 

Figure 1. - Typical solid propellant rocket motor. 


NASA-S-66-2022 FEB 24 



(c) Burning surface geometry of a typical internal -burning six-point 
star propellant grain configuration versus regressed distances 
normal to the original propellant surface. 


Figure 1. - Concluded 





Propellant burning area per nozzle throat area 


NASA-S-66-2033 FEB 24 







NASA-S-66-2032 FEB 24 



Figure 3. - Typical performance variation of a 
solid propellant rocket motor of fixed geometry 
and given propellant that is tested at differing 
prefire propellant temperatures. 
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Increasing distance burned normal to the 
original propellant surface, in. 


Figure 4. - Propellant burning surface area versus 
increasing distance burned normal to the original 
propellant surface for the six-point star propellant 
grain configuration shown in figures 1(b) and 1(c). 
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Figure 5. - Area ratio of propellant burning surface area 
to nozzle throat area versus increasing distance burned 
normal to the original propellant surface for the solid 
propellant rocket motor shown in figure 1. 
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Chamber pressure, psia 
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Prefire propellant temperature, °F 
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Prefire propellant temperature, °F 


Figure 7.- Second- order, least- squares curve 
fit of experimentally determined chamber 
pressures at the individual regressed distance 
of interest normal to the original propellant 
surface versus prefire propellant temperature. 


Figure 8.- Second-order, least-squares curve 
fit of experimentally determined burning rates 
at the individual regressed distance of interest 
normal to the original propellant surface versus 
prefire propellant temperature. 


Motor operating time, sec 
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u « X ± sK 2 = upper and lower limits 

in 



Figure 10. - Example of two-sided tolerance limits 
for a normal distribution. 


Figure 9.- Second-order, least-squares curve fit 
of experimentally determined motor operation 
times to the individual regressed distance of 
interest normal to the original propellant 
surface versus prefire propellant temperature. 
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u u * X + s ^i = upper limit 


to 
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(a) Upper tolerance limit 
* X - sK^ s= lower limit 


to 

c 

<U 



Figure 11. - Example of one-sided tolerance limits for 

normal distribution. 
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Figure 12. - IBM 7094 deck setup for the general solid propellant rocket 
motor performance computer program using the group transformation 
method. 
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13. - Example of choosing 
variation of chamber p: 


Figure 

05 




CT> 

to 



Figure 14. - Uni vac 1107 deck setup for the solid propellant rocket motor 
performance versus time computer program using the group transformation 
method. 
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>f choosing the input points that will best define the variation 
tamber pressure and thrust versus time. 
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(a) Typical transformed time and chamber pressure graphed outputs. 

Figure 16. - Example of graphical format from the solid propellant 
rocket motor performance versus time computer program. 
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APPENDIX A 

LISTING OF THE GENERAL SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
COMPUTER PROGRAM USING THE GROUP TRANSFORMATION METHOD 
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o n n 


SIB JOB DRIVFR GO.SOURCE 
SIBFTC DRIVFR 


DRIVER PROGRAM FOR MOTORS 
DIMENSION XX (4R60 ) 

DIMENSION EXPA (4) » C0UNT(4>» WP (4.15). WE8(4.15). ATT0T(4»16). 

1 AETOT (4.15) . P(9.4.15). M0TN0(4.1b). FVAC1 (9.4. 15) 

COMMON G» CPHI . NT. EXPA. COUNT. WP. WEB. ATTOT. AETOT. 

1 P. MOTNO. FVAC1 . NTIMES. XX. NTOT . NBURN. NFIRST . NP 

COMMON CPI. CP2. CT1 . CT2. VEC1. VEC2. VEC3 

1 READ (5.10) CASE 
WRITE (5.11) CASE 
WRITE (6.21) 

READ (5. BO) 

WRITE ( 6 ► BO ) 

READ (5.79) NTIMES 
WRITE (6.14) NTIMES 
DO 2 1 = 1. NTIMES 

READ (5, BO) 

2 WRITE (6. BO) 

10 FORMAT (A6) 

11 FORMAT (1H1. A6/> 

14 FORMnT (IX //I2 * 33H TIMES WERE USED IN THIS ANALYSIS ) 

21 FORMAT ( 46H0GENERAL SOLID-FUFL ROCKeT PERFORMANCE PROGRAM / 

1 120H ROCKET TEST FIRINGS CONDUCTED AT SEVERAL TFMPERATURES 

2 ARE EFFECTIVELY TRANSFORMED TO SPECIFIFD TEMPERATURES OF INTEREST 

3 / B9H STATISTICAL ANALYSES ARE PERFORMED ON THE TEST FTRING 

4 DATA AFTER IT HAS BEEN TRANSFORMED //) 

79 FORMAT (10X. 12) 

80 FORMAT (BOH 

1 ) 

CALL MOTORS 
CALL MOT 2 
CALL M0T3 
60 TO 1 
END 


SIBFTC STATS 

SUBROUTINE STATS (X. N. XM. SD) 
DIMENSION X(100) 

SIJM1 = 0. 

SUM2 = 0. 

XNO = N 

XNOl = XNO - 1. 

DENOM = XN0*XN01 
DO 10 I = 1 »N 
SUM1 = SUM1 + X ( I ) 

10 SUM2 = SUM2 + X ( I ) *X ( I ) 

XM = SUM 1 / XNO 

XNUM = (XN0*SUM2 - SUM1+SUM1) 

IF (XNUM .LE. 0.) GO TO 14 
SD = SORT ( XNUM/DENOM ). 

IF ((SD/XM) .LT. .00025) GO TO 14 
60 TO 15 

14 SD = 0. 

15 RETURN 
END 


TO 



*h c\] ro * h cvi fO HC\jfO 


ALPHA 


I 


SORIGIN 
SIBFTC MOTORS 

SUBROUTINE MOTORS 
C 

C SOLID PROPELLANT ROCKET MOTOR TEST PROGRAM - FIRST HALF 

C 

DIMENSION ZZ (648) » XXC4860) 

COMMON G» CPHIf NT » EXPA » COUNT. WP. WEB. ATTOT. AETOT. 

1 P. MOTNO. FVAC1 » NTIMES. XX. NTOT. NBURN. NFIRST . NP 

COMMON CPI. CP2. CT1» CT2. VEC1. VEC2. VEC3 

DIMENSION EXTLMt 4). EXPA( 4). EPS( 4). 

COUNT ( 4). SET1<45). SET2C45) » SET3(45). SET4(45)» 

SETS (45) » SET6(45) . SET7(45) . TGRP(3). 

TENTRY (60 ) . TENT2<60). ARRAY(46.4). ENTRY(9»60). 

C0EFK9.4). C0EF2(9.4) . IFLAG1 (9) . IFLAG2(9) 

DIMENSION X (4 ) » TIME (9.4 » 15) . INTP(9.4»15) . P(9.4.15). 

INTF (9.4.15) . F(9»4.15). PK9.4.15). P2(9»4»1S). 

P3(9.4.15> . TK9.4.15). T2(9.4.1S). T3(9.4»1S). 

INTPT1 (9.4. 15) » INTPT2 (9.4.15) . INTPT3 (9.4. IS) 

DIMENSION ATK4.15). AT2(4.15)» AT3(4.15). AT4(4»15). AEK4.15). 

AE2 (4.15) » AE3 (4.15) . AE4(4.15). WP(4.15)» WEB(4.15>. 

DENS (4.15) . PAF (4.15) . LNP(9.4). LNT(9»4>» 

M0TN0(4.15) .TEMP(4.15) . RU (4.15). KN (4.15) 

C 

DIMENSION MP (9.4) . SDP(9*4>. MF(9.4). SDF(9.4). 

1 MTIME (9.4) . SDTIME (9.4) . MPTK9.4). SUPT1 (9.4) . MPT2(9.4). 

2 SDPT2 (9.4) . MPT3(9.4>» SDPT3(9»4>» MTTK9.4). SDTT1 (9.4) . 

3 MTT2 (9.4) . S0TT219.4). MTT3(9.4)» SDTT3(9.4). MPTT(9.4). 

4 SDPTT (9.4) . MTTT (9.4) . SDTTT(9.4). MINTK9.4). SDINT1 (9.4) . 

5 MI NT 2 (9.4) » S0INT2(9.4). MINT3(9»4). SDINT3(9»4)» MlNTT(9»4)» 

6 SDINTT (9.4) » MTEMP < 15) . SDTEMP ( 15) 

C 

C 

DIMENSION R 1 (4.15). R2 (4.15). R3 (4.15). AFT0T(4»15). 

1 ATTOT (4.15) »EA (4.15). EB (4.15). EC (4.15). ED (4. 15). 

2 ETOT (4.15) 

3 »C1MIN(4). C1MAX (4 ) . C2MIN(4). C2MAX(4)» C1TMN(4). C1TMX(4)» 

4 C2TMN ( 4 ) . C2TMX ( 4 ) 

DIMENSION C1MN ( 5 ) . C1MX(5)» C2MN(5). C2MX(5) 

DIMENSION 

1 FVAC1 (9.4.15) . WPU (9.4.15). ISP (9.4.15). IVAC1 (9.4. 15) » 

2 ISP1 (9.4.15). 112 (9.4.15). F12 ( 9 . 4 » 15) . 1SPU12 (9, 4 . 15) » 

3 113 (9.4.15). F13 (9.4. 15) . ISPU13 (9.4.15). 114 (9.4.15). 

4 F14 (9.4.15) .ISPU14(9. 4.15) . CSTU (9.4.15) 

DIMENSION MAT (5) . SDAT<5)» MAE (5 ) . SDAE(5>» ME(5)» SDE(5) 
DIMENSION MINTP (9.4) . SDINTP<9»4)» MINTF19.4). SDINTF(9.4)» 

1 MWPU (9.4) . SUWPU (9.4) . MISP(9.4). SDISP(9.4). MCSTU(9.4)» 

2 SDCSTU (9.4) . MRU ( 4). SDRU( 4). MKN( 4). SDKN( 4). 

3 MIVC1 (9.4) . SDIV/C1 (9.4) . MFVC1 (9.4) » SDFVC1 (9.4) ► MIU1K9.4). 

4 SDIU11 (9.4) . MI12(9.4). SDI12(9.4)» MF12(9.4). SDF12(9.4). 

5 MIU12 (9.4) . SDIU12(9. 4) » MI13(9.4). SDI13(9.4)» MF13(9.4>. 

6 SDF13 (9.4) . MIU13(9»4). S0IU13(9.4). MI14(9.4). SDI14(9.4). 

7 MF14 (9.4) . SDF14 (9.4) . MIU14(9»4)» SDIU14(9»4). IL(4) 

DIMENSION SIGP12 (9) . SIGP22(9)» SIGP32(9). SIGP11 (9) * 

1 SIGP2K9). SIGP3K9). SIGP13(9). SIGP23(9). SIGP33(9>» 

2 PIK11 (9) . PIK2K9). PIK31(Q), PIK12(9). PIK2?(9). PIK32(9>. 

3 PIK13 (9) . PIK23 (9) . PIK33(9) 

C 

EQUIVALENCE 
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n n n non 


(XX< 1) 
(XX( 541) 

(xxaoai) 

(XX (1621) 
(XX (2161 ) 
( XX (2701 ) 
( XX ( 3241 ) 
(XX ( 3781 ) 
( XX (4321 ) 
EQUIVALENCE 


1 

2 

3 

4 

5 

6 

7 

8 
9 


1 

2 

3 

4 

5 

6 

7 

8 
9 


(ZZ( 

(ZZ( 

(ZZ( 


1) 

37) 

73) 


(ZZ< 109) 
<ZZ( 145) 
(ZZ( 181) 
(ZZ< 217) 
(ZZ( 253) 
(ZZ( 289) 


EQUIVALENCE 
(ZZ( 325) 
(ZZ( 361) 
(ZZ< 397) 
(ZZ( 433) 
(ZZC 469) 
(ZZ( 505) 
<ZZ( 541) 
(ZZ( 577) 
(ZZ( 613) 


1 

2 

3 

4 

5 

6 

7 

8 
9 


I12H) * 113(1) * 114(1) * Pl(l))* 

WPU (1)* P2 ( 1 ) ) * 
ISP (1)* P3 ( 1 ) ) *■ 
F12 ( 1 ) * F13(l) * F14 ( 1 ) f T1 ( 1 ) ) * 

INTF (1 ) * T2 ( 1 ) ) # 
INTP ( 1 ) * T3( 1 ) ) * 

ISPlIllr ISPU12 ( 1 ) * ISPU13 ( 1 ) * ISPU14(1). INTPTl(l) 

CSTU ( 1 ) » INTPT2 ( 1 ) 
INTPT3 ( 1 ) 

MIVC1 (1 ) • MI 12 ( 1 ) * MI 13 ( 1 ) * MI14 ( 1 ) * MPTt(l))* 

MPT? ( 1 ) ) * 

MPT3 ( 1 ) ) * 

MFVC1 ( 1 ) » MF12 ( 1 ) r MF13(1)* MF14(1)» MTTl(l))* 

MTT? ( 1 ) ) * 

MTT3 ( 1 ) ) . 

MIU11 ( 1 ) *MIU12 ( 1 ) *MI(J13 (1 ) *MlU14( 1 ) fMINTl (1) ) * 

MINT? ( 1 ) ) * 

MINT3 ( 1 ) ) 

smvci(i>» sdi 12 ( i ) * sni 13 ( i ) * snii4(i)f snPTi ( 1 ) ) » 

SDPT2 < 1 > ) . 
S0PT3 ( 1 ) ) . 

SDFVC1 (l)r SDF12(1). SDF13(l)r SDF14(1). SDTT1 ( 1 > ) * 

S0TT2U) ) . 
SDTT3 < 1 ) ) * 

SOIUll(l) * SDIU12 ( 1 ) * SDIU13 ( 1 ) * SDIU14 ( 1 ) * SDPTT < 1 ) ) . 

SOTTT(l) ) • 
S0INT3 ( 1 ) ) 



INTEGER 

COUNT 9 

ENDMOT. TEMSEL 









REAL 

INTP 

* INTF * INTPT1 


lNTPT?f 

INTPT3 


INTPTT 




1 

MTEMP f 

MTIME 

* SDTIME* MP 


ME * 

MPT1 


MPT2 


MPT3 


2 

MPTT # 

MTT1 

* MTT2 * MTT3 


MTTT r 

MINT1 


MINT? 


MI NT 3 


3 

MINTT t 

LNT 

* LNP » MOTNO 


ISP r 

IVAC1 


ISP1 


112 


4 

I SPU1 2 r 

113 

* ISPU13 . 114 


I SPU1 4 9 

KN 


MINTP 


MINTF 


5 

MWPU r 

MISP 

* MCSTU * MRU 


MKN 9 

MIVC1 


MEVC1 


MIU11 


6 

MI 12 t 

MF12 

* MI 13 * MF13 


MIU13 f 

MI 14 


MF 14 


MHJ14 



REAL MAT t MAE r ME 


) t 
) » 
) 


INPUT FORMATS 


81 

FORMAT 

(112. 4F12.4) 

82 

FORMAT 

(2F12.5* 312) 

83 

FORMAT 

(4F12.4) 

84 

FORMAT 

( A6 * P'12. 4. 12) 

85 

FORMAT 

(4E12.5) 

86 

FORMAT 

(8F9.4* F8.4) 


OUTPUT FORMATS 

16 FORMAT ( 10H0THERE ARE 12* 46H AMdltNT PRESSURES BEING CONSIDERED* 

1 THEY ARE / 4F12.4) 

18 FORMAT (10H0THERE ARE 12* 46H EXACT TFMPERATURES BEING CONSIDERED* 
1 THEY ARE / 4F12.4) 

20 FORMAT ( 41H G AND PHI ARE CONSTANT FOR A RUN AND ARE / 2F12.5) 

2001 FORMAT ( 62H1THR0AT AREA ATA ATB ATC ATP 

1 AT TOTAL / 10H MOTOR NO.//) 

2002 FORMAT (11H0TEMP GROUP 12 // (IX* A6*5 X* 5F10.4)) 

2003 FORMAT ( 62H1EXIT AREA AEA AEB AEC AED 

1AE TOTAL / 10H MOTOR NO.//) 

2004 FORMAT ( 62H1EXPANSI0N RATIO EA EB EC En 

1 E TOTAL / 10H MOTOR NO.//) 
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2005 FORMAT (101H1TIMES Tl T2 T3 T4 

1 T5 T6 T7 T8 T9 / 10H MOTOR NO.) 

2006 FORMAT (11H0TEMP GROUP 12 // (IX. A6.5X* 9F10.4)) 

2007 FORMAT ( 101H1PRESSURES PT1 PT2 PT3 PT4 

1 PT5 PT6 PT7 PT8 PT9 / 10H MOTOR NO.) 


2008 FORMAT C119H1 INTEGRATED 


INTPT1 


INTPT2 


INTPT3 


1NTPT4 INTPT5 INTPT6 INTPT7 

2 9H PRESSURE / 10H MOTOR NO.) 

2009 FORMAT (102H1 CSTAR Tl CSTAR 

1CSTAR T5 CSTAR T6 CSTAR T7 CSTAR T8 CS 

2010 FORMAT ( 52H1 AVERAGE R BASE! 

1 48H ON BURN TIME 

2 10H MOTOR NO.) 

2012 FORMAT ( 101H1PROPELLANT WT. WPTl WP1 

1 WPT5 WPT6 WPT7 WPT8 

2014 FORMAT (1H1. 30Xr 13H PA = AMBIENT //) 


INTPT8 


INTPT9/ 


CSTAR Tl CSTAR T2 CSTAR T3 CSTAR T4 
CSTAR T7 CSTAR T8 CSTAR T9 /10H MOTOR NO.) 

AVERAGE r BASED average kn based 

ON BURN TIME ON BURN TIME / 

0 .) 

ANT WT. WPTl WPT2 WPT3 WPT4 

WPT7 WPT8 WPT9 / 10H MOTOR NO.) 


2015 FORMAT 

( 1 1 9H THRUST 


FT1 


FT2 


FT3 

1FT4 

FT5 

FT6 


FT7 


FT8 


2 1DH MOTOR NO.) 







2016 FORMAT 

(119H TOTAL 

IMPULSE 

IT1 


IT2 


IT3 

* IT4 

ITS 

IT6 


IT7 


IT8 


2 10H MOTOR NO.) 







2017 FORMAT 

( 1 19H SPEC. 

IMPULSE 

ISPT1 


ISPT2 


TSPT3 

1SPT4 

ISPT5 

ISPT6 


ISPT7 


ISPT8 


2 10H MOTOR NO.) 







2024 FORMAT 

(1H1. 30X . 

5H PA = F 

9.3) 






ISPT9 


2030 FORMAT <1H1» 30X. 18H TRANSFORMED TIMF / 9H TIME NO. I?) 

2031 FORMAT (1H1. 30X. 22H TRANSFORMED PRESSURES / 9H TIME NO. 12) 

2040 FORMAT (11H0TEMP GROUP 12 //(IX. A6.18X. 4F10.4)) 

2041 FORMAT (11H0TEMP GROUP 12 // (IX. A6. 18X. 4F12.2. 10X. 4F12.2)) 

2042 FORMAT ( 30X. 18H PRESSURE INTEGRAL 40X. 17H AVERAGF PRESSURE) 

2044 FORMAT (11H0TEMP GROUP 12 // (1X.A6.7X. F10.4. 12X. F10.4)) 

2045 FORMAT (11H0TEMP GROUP 12 // (IX. A6r 5X. 9F12.2)) 

2046 FORMAT (11H0TEMP GROUP 12 // (IX. A6.5X. 9F10.2)) 

2047 FORMAT ( 27X.36H TEMPI TEMP2 TEMP3 TEMP4 / 


2047 FORMAT ( 27X.36H TEMPI 1 

1 10H MOTOR NO.) 

2048 FORMAT ( 27X.43H TEMPI 

1 43H TEMPI 

2 10H MOTOR NO.) 


TEMP2 

TEMP2 


TEMP3 

TEMP3 


2050 FORMAT (1H1. 30X. 23H TRANSFORMATION FACTORS ) 


2051 FORMAT ( 7H PIK11 9F1P.6 / 7H 
1 7H PIK12 9F10.6 / 7H 


PIK21 9F10.6 
PIK22 9F10.6 


PIK13 9F1D.6 / 7H PIK23 9F10.6 / 7H 


TEMP4 

TEMP4 


PIK31 9F10.6 / 
PIK32 9F10.6 / 
PIK33 9F10.6//) 


2052 FORMAT ( 7H SIGP11 9F10.6 / 7H SIGP21 9F10.6 / 7H SIGP31 9F10.6 / 

1 7H SIGP12 9F10.6 / 7H SIGP22 9F10.6 / 7H SIGP32 9F10.6 / 

2 7H SIGP13 9F10.6 / 7H SIGP23 9F10.6 / 7H SIGP33 9F10.6 ) 


2055 FORMAT (12H0MEAN 
2060 FORMAT (12H0MEAN 


12H STAND. DEV 
12 Ft STAND. DEV. 


9F10.4) 

9F10.2) 


2062 FORMAT ( 34H0 CONFIDENCE ON NORMAL DISTRIBUTION / 


1 

25H 

(ONE 

MIN 

4F10.4 

/ 

2 

25H 

SIDED 

MAX 

4F10.4 

// 

3 

25H 

(TWO 

MIN 

4F10.4 

/ 

4 

25H 

SIDED 

MAX 

4F10.4 

) 


2063 FORMAT ( 14H0MEAN 


1 14H STANDARD DEV. F 

2064 FORMAT (25H0MEAM 

1 25H STANDARD DEV. 

2065 FORMAT ( 25HOTOT AL MEAN 

1 25H TOTAL STANDARD DEV 

2066 FORMAT (25H0MEAN 

1 25H STANDARD DEV. 

2067 FORMAT (25HOTOTAL MEAN 


F10. 4. 12X . F10 .4 / 

F10. 4. 12X » F10.4 ) 

4F10.4 / 

4F10.4) 

4F10.4 / 

V. 4F10.4 ) 

4F12.2. 10X> 4F12.2 / 
4F12.2. 10X. 4F12.2 ) 
4F12.2. 10X. 4F12.2 / 
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1 25H TOTAL STANDARD DEV. 4F12.2* 10X* 4F12.2 ) 

2068 FORMAT < 34HOCONFIDENCE ON NORMAL DISTRIBUTION / 

1 25H (ONE MIN 4F12.2* 10X* 4F12.2 / 

2 25H SIDED) MAX 4F12.2 * 10X* 4F12.2 // 

3 25H (TWO MIN 4F12.2* 10X*4F12.2 / 

4 25H SIDED) MAX 4F12.2* 10X * 4F12.2 ) 

2070 FORMAT ( 12HOTOTAL MEAN 5F10.4 / 12H TOTAL S.D. 5F10.4) 

2071 FORMAT ( 34HOCONFIDENCE ON NORMAL DISTRIBUTION / 

1 12H (ONE MIN 5F10.4 / 12H SIDED) MAX 5F10.4// 

2 12H (TWO MIN 5F10.4 / 12H SIDED). MAX 5F10.4) 

2080 FORMAT ( 12H0MEAN 9F12.2 / 12H STAND. DEV. 9F12.2) 

C 

COUNT (1) = 0 
COUNT (2) = 0 
COUNT (3) = 0 
BLANK = 0. 

C 

C READ DATA AND GROUP IT INTO TEMPERATURE GROUPS 

C 

5 READ (5*81) NP r (EXPA(I)* I = 1*NP) 

WRITE (6*16) NP* ( EXPA ( I ) * I = 1*NP) 

READ (5*81) NT* (EXTEM ( I ) * I=1*NT) 

READ (5*81) NT* ( TGRP (I)* I=1*NT) 

TGRP (NT + 1) = TGRP (NT) * 2. 

WRITE (6*18) NT* (EXTEM ( I ) * I = 1*NT) 

READ (5*82) G* PHI* NTQT* NBURN* NFIRST 
WRITE (6*20) G* PHI 
READ (5*83) CT1 * CT2* CPI* CP2 
C 

C ASSIGN DATA FOR A MOTOR INTO PROPER TEMPERATURE GROUP 

C 

9 READ (5*84) XMOT * FIRETP* ENDMOT 
DO 19 J = 1*3 

EPS ( J ) = (TGRP (J+l) - TGRP ( J) ) *.5 
IF ((TGRP (J) + EPS ( J ) ) .GT. FIRETP) GO TO 25 
19 CONTINUE 

25 COUNT ( J ) = COUNT ( J ) + 1 
K = COUNT ( J) 

45 TEMP (J*K) = FIRETP 
MOTNO ( J * K ) = XMOT 
C 

C READ REMAINDER OF DATA FOR A MOTOR 

C 

READ (5*85) AT1(J*K)* AT2(J*K)* AT3(J*K)* AT4(J*K)* 

1 AEKJ*K)* AE2 ( J*K ) * AE3 ( J * K ) * AE4(J*K)* 

2 WP ( J * K ) * WEB ( J * K ) * DENS ( J * K ) * PAF(J*K) 

READ (5*86) (TIME ( I * J*K ) * I = 1*NTIMES) 

READ (5*86) < INTP ( I * J*K ) * I = 1*NTIMES) 

READ (5*86) ( P(I*J*K)* I s 1*NTIMES) 

READ (5*86) ( INTF ( I * J * K ) * I = 1*NT1MLS) 

READ (5*86) ( F(I*J*K>* I = l*NTlMtS> 

IF (ENDMOT .LT. 1) GO TO 9 
C 

C PROCEED WITH MAIN BODY OF PROGRAM 

C 

Jl = COUNT (1) 

J2 = COUNT ( 2 ) 

J3 = COUNT (3) 

TEMSEL = 1 
CPHI = COS (PHI) 

C 

C COMPUTE P(I*J*K) * F ( I * J*K ) AND STATISTICS ON P*F AND TIME 
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non o o n 


DO 115 I = 1*NTIMES 
ISW = 1 

IF (P(I*J*K) .GT. 0.) ISW = ISW + 1 
IF (F ( I t J*K ) .GT. 0.) ISW = ISW + 2 
DO 115 J = 1 *NT 
L = 0 

JJ = COUNT (J) 

DO 108 K = 1*JJ 

GO TO ( 104 * 105 * 106 * 107)* ISW 

104 P ( I * J*K ) = INTP ( I * J*K ) / TIME<I*J*K) 

105 F(I*J*K) = INTF(I*J*K) / TIME(I*J*K) 

GO TO 107 

106 P(I*J*K) = INTP(I*J*K> / TIME ( I * J*K ) 

107 L = L + 1 

SET! (L ) = P(I*J*K) 

SET2 (L ) = F(I*J*K) 

108 SET3 (L ) = TIME ( I * J*K ) 

CALL STATS (SET1 * L* MP (I*J)* SUP <I*J)) 

CALL STATS (SET2* L* MF CI*J)* SDF (I*J)) 

115 CALL STATS (SET3* L* MTIME < I * J) * SUTIMF(IrJ)) 

WRITE ( 6 * 2005 ) 

DO 55 J = If NT 
JJ = COUNT ( J ) 



WRITE 

(6*2006) 

J* 

( MOTNO ( J * K ) * 

(TIME(I.J.K) .1=1.9) . 

K = 

1*JJ) 

55 

WRITE 

(6*2055) 

(MTIME ( I * J) * 1=1 

*9)* (SDTIME(I.J) *1= 

l*o) 


116 

WRITE 

(6*2007) 








DO 118 

J = 1 r NT 







JJ = COUNT (J) 








WRITE 

(6*2046) 

J 9 

(MOTNO ( J *K ) * 

(P(I.J.K). I 

= 1*9) * 

K = 

1*JJ) 

118 

WRITE 

(6*2060) 

(MP(I.J)* 1=1.9) 

* (SOP ( I * J) * 

1=1*9) 




WRITE 

(6*2014) 







216 

WRITE 

(6.2015) 








DO 218 

J = 1 *NT 







JJ = COUNT (J) 








WRITE 

(6*2045) 

J 9 

( MOTNO (J.K ) * 

(F(I*J*K)» I 

= 1*9)* 

K = 

1*JJ) 

218 

WRITE 

(6*2080) 

(MF(IrJ)r 1=1*9) 

* (SDF ( I * J) * 

1=1*9) 




COMPUTE STATISTICS FOR 1 1 MPER ATURE 


DO 270 J = 1 r NT 
L = 0 

JJ = COUNT ( J ) 

DO 260 K = lrJJ 
L = L + 1 

260 SETKL) = TEMP(J*K) 

270 CALL STATS (SET1 * L * MTEMP ( J) * SUTLMP(J)) 

COMPUTE AETOTAL* ATTOTAL AND AREA RATIOS FOR FUTURE USE 


L = 0 

DO 704 J = 1*NT 
JJ = COUNT ( J ) 

DO 704 K = 1 * JJ 
L = L + 1 

AETOT ( Jr K ) = AEKJfK) + AE2<J*K) + AE3 ( Jr K ) + AE4 ( J* K ) 
ATTOT(JfK) = AT1(J»K) + AT2(J*K) + AT3<J*K) + AT4(J*K) 
EA ( J * K > = AE1(J*K) / AT1(J*K) 

EB ( J * K ) = AE2 ( J# K ) / AT2 ( J*K ) 

EC(JrK) = AE3(J*K) / AT3 ( J*K ) 

ED(JrK) = AE4 ( J * K ) / AT4 ( J*K ) 

ETOT(JrK) = AETOT ( J r K ) / ATTOT ( J*K ) 
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n n 


SETKL) = ATI ( J»K) 

SET2(L) = AT2<J»K> 

SET3 (L ) = AT3 ( J»K ) 

SET4 (L ) = AT4(J*K> 

SET5 (L ) = ATTOT(J.K) 

SET6 (L) = AE1 ( J»K) 

SET7 ( L ) = AE2(J»K> 

704 TENT2(L>= AE3(J»K> 

CALL STATS (SETl. L. MAT(l), SDAT(l) ) 

CALL STATS (SET2* L» MAT (2) * S0AT<2>> 

CALL STATS (SET3* L» MAT (3) r SDAT<3>> 

CALL STATS CSET4* L» MAT(4)» SDATt4)) 

CALL STATS <SET5» L» MAT(5>» SDAT<5>) 

CALL STATS (SETS* L» MAE ( 1 ) * SDAEtl) ) 

CALL STATS (SET7» L» MAE<2)» SDAE(2)) 

CALL STATS (TENT2*L» MAE (3). SDAE(3)) 

L = 0 

DO 1704 J = 1 »NT 
JJ = COUNT (J) 

DO 1704 K = 1*JJ 
L = L + 1 
SETl(L) □ EA(JrK) 

SET21L) = EB ( J .K ) 

SET31L) = ECIJ.K) 

SET4IL) = ED ( J.K ) 

SETS (L ) = ETOT(J.K) 

SET6 (L ) = AE4 ( J»K) 

1704 SET7 (L ) = AFTOT ( J»K ) 

CALL STATS (SET1 t Lr ME(1)> SDE(1>) 

CALL STATS (SET2f Lr ME(2)> SDE(2)) 

CALL STATS <SET3» L» ME(3)» SDE(3)) 

CALL STATS (SET4* L» ME(4)r SDE(4)) 

CALL STATS <SET5» L* ME(5>» SDE<5>> 

CALL STATS (SET6» L» MAE(4>» SDAE(4)) 

CALL STATS (SET7 • L. MAE<5>> SDAE(5)) 

DO 1706 K = 1 f 5 

C1MN(K) = MAT (K) - CP1*S0AT(K> 

C1MX(K> = MAT (K ) + CP1*SDAT(K> 

C2MN ( K ) = MAT t K ) - CP2*SDAT(K) 

1706 C2MX (K ) r MAT t K ) + CP2*SDATtK) 

PERFORM CALCULATIONS ON NON-TRANSFORMED DATA 
C 

DO 1100 I = NFIRST rNTIMES 
DO 1100 J = 1 »NT 
L = 0 

JJ = COUNT (J) 

DO 1076 K = 1»JJ 
L = L + 1 

IF ( INTP ( I r J r K ) .NE. 0.) GO TO 1070 
INTP ( I » J*K) = P(I»J»K) * TIME ( I r J*K ) 

1070 IF ( INTF ( I r J.K ) .NE. 0.) GO TO 1075 
INTF ( I » J»K ) = F(I.d#K) * TIME ( I » J» K ) 

1075 WPU ( I ► J»K ) = WP ( J » K ) * I NTP ( I . J » K ) / INTP (NTOT • J.K ) 

CSTU ( I * J.K ) = INTP ( I * J.K) *G*ATTOT ( J.K ) / WPU(I.J.K) 
ISP ( I . J » K ) = INTF ( I . J.K ) / WPU(I.J.K) 

FVAC1(I.J.K)= F(I.J.K) + PAF ( J .K ) *ALTOT (J.K) *CPHI 
IVAC1(I.J.K>= FVACl(I.J.K) * TIMt(I.J.K) 

ISP1 ( I . J.K ) = IVAC1 ( I . J.K ) / WPU ( I . J . K ) 

TENT2(L)= CSTU ( I » J » K ) 

SET 1 (L ) = IVAC1 ( I » J.K) 

SET2 (L ) = INTP ( I . J »K ) 
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"SET3 (L ) = INTF ( I , J ,K ) 

SET4 (L ) = FVAC1 C I , J»K) 

SET5(L) = ISP(I,J,K> 

SET6CL) = WPU(I,J,K) 

1076 SET7 (L ) = ISP1 (I,J,K) 

C ALL STATS ( TENT2 , L t MCSTU ( I , J ) , 


CALL STATS (SET1, L, 
CALL STATS (SET2, L, 
CALL STATS (SET3, L, 
CALL STATS (SET4, L, 
CALL STATS (SET5, L, 
CALL STATS (SET6, L* 
CALL STATS (SET7, Lr 
1100 CONTINUE 


MIVC1(I,J), 
HINTP(I ,J>, 
MINTF ( I , J) , 
MFVC1 ( I , J) , 
MISP (I,J), 
MWPU (I,J)t 
MIU11(I,J ) 9 


C 


SDCSTU ( I , J) ) 
SDIVC1 ( I , J) ) 
SDINTP( I , J) ) 
SOINTF ( I , J) ) 
SOFVC1 ( I , J) ) 
SDISP (I,J>> 
SDWPU (I,J)) 
SDIU11 ( I f J) ) 


WRITE (6,2008) 

DO 1080 J = 1,NT 
JJ = COUNT (J) 

WRITE (6,2045) J* (MOTNO(J,K), ( INTP( I * J,K) , I = 1,9)# K 
1080 WRITE (6,2080) (MINTP(I,J), 1=1,9), (SDINTP ( I , J) , 1=1,9) 
WRITE (6,2014) 


1,JJ> 


WRITE (6,2016) 

DO 1085 J = It NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), ( INTF ( I , J,K ) , I = 1,9), K = 1,JJ> 
1085 WRITE (6,2080) (MINTF ( I r J) , 1=1,9), (SDINTF ( I , J) , 1=1*9) 

WRITE (6,2014) 

WRITE (6,2017) 

DO 1090 J = 1,NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), (ISP(I,J,K>, I = 1,9), K = 1,JJ> 
1090 WRITE (6,2080) (MISP (I,J>, 1=1,9), (SDISP (I,J), 1=1,9) 

WRITE (6,2024) EXPA(l) 

WRITE (6,2015) 

DO 1061 J = 1 ,NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), (FVAC1 ( I , J»K ) , 1=1,9), K=1,JJ) 

1061 WRITE (6,2080) (MFVC1(I,J), 1=1,9), (SDFVC1 ( I , J) , 1=1,9) 

WRITE (6,2024) EXPA(l) 

WRITE (6,2016) 

DO 1062 J = 1,NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), ( IVAC1 ( I , J,K) , 1=1,9), K=1»JJ) 

1062 WRITE (6,2080) (MIVC1(I,J), 1=1,9), (SDIVC1 ( I , J) , 1=1,9) 

WRITE (6,2024) EXPA(l) 

WRITE (6,2017) 

DO 1063 J = 1,NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), (ISP1 ( I , J,K) , 1=1 ,9) , K=1,JJ) 

1063 WRITE (6,2080) <MIU11(I,J), 1=1,9), (SDIU11 ( I , J) , 1=1,9) 

C 

IF (NP .LE. 1) GO TO 1099 
DO 1101 I = NFIRST ,NTIMES 
DO 1101 J = 1 ,NT 
L = 0 

JJ = COUNT (J) 

DO 1160 K = 1,JJ 
L = L + 1 

112(1, J,K) = I VAC1 ( I , J ,K ) - (EXPA(2)*AET0T(J,K)*CPHI)*TIME(I,J,K) 
F12 ( I , J,K) = 112(1, J,K) / TIME ( I , J,K ) 

ISPU12 ( I , J,K) = 112(1, J,K) / WPU( I , J,K ) 

SETKL) = 112(1 * J,K) 



SET4 (L ) = F12(I*J*K) 

1160 SET7(L) = ISPU12 (I * J*K) 

CALL STATS (SET1* L* MI12d.J). SDI12(I*J)) 

CALL STATS (SET4* L. MF12(I*J). SDF12(I*J)) 

1101 CALL STATS (SET7* L. MIU12(I*J)* S0IU12d*J>) 

C 

WRITE (6*2024) EXPA(2> 

WRITE (6*2015) 

DO 1071 J r 1,NT 
JJ = COUNT (J) 

WRITE (6*2045) J* (MOTNO(J*K)* (F12(I*J*K). 1=1*9)* K=1*JJ) 

1071 WRITE (6*2080) <MF12(I.J). 1=1*9)* (SDF12(I*J)* 1=1*9) 

WRITE (6*2024) EXPA(2> 

WRITE (6*2016) 

DO 1072 J = 1*NT 
JJ = COUNT (J) 

WRITE (6*2045) J* (MOTNO(J*K)* (I12(I*J*K>* 1=1*9)* K=1*JJ) 

1072 WRITE (6*2080) (MI12(I*J)* 1=1*9)* (SDI12(I*J)* 1=1*9) 

WRITE (6*2024) EXPA<2) 

WRITE (6*2017) 

DO 1073 J = 1 *NT 
JJ = COUNT (J) 

WRITE (6*2045) J* (MOTNO(J*K>* ( ISPU12 < I * J* K ) * 1 = 1*9)* K=1*JJ> 

1073 WRITE (6*2080) ( MIU12 ( I * J) * 1=1 * 9 ) * (SDIU12( I * J) * 1=1*9) 

IF (NP .LE. 2) GO TO 1099 

C 

DO 110? I = NFIRST *NTIMES 
DO 1102 J = 1 *NT 
L = 0 

JJ = COUNT (J) 

DO 1180 K = 1*JJ 
L = L + 1 

113 ( I * J*K ) = I VAC1 ( I * J*K ) - (EXPA (3)*AET0T(J.K)*CPHI)*TIME<I. J*K) 

F13(I*J*K) = 113(1. J*K) / TIME ( I * J*K ) 

ISPU13 ( I * J *K ) = 1 13 ( I * J*K ) / WPU ( I * J*K ) 

SETKL) = 113(1 *J*K) 

SET4 (L ) = F13 ( I . J.K ) 

1180 SET7 (L ) = ISPU13 ( I . J »K ) 

CALL STATS (SET1. L* MI13(I*J)* SDI13(I*J)) 

CALL STATS (SET4. L* MF13(I*J )* SDF13(I*J)) 

1102 CALL STATS (SET7* L* MIU13(I*J)* SUlU13(I*J)) 

C 

WRITE (6*2024) EXPA(3) 

WRITt (6.2015) 

DO 1081 J = 1 *NT 
JJ = COUNT (J) 

WRITE (6*2045) J* (MOTNO ( J* K ) * (F13(I*J*K)* 1=1*9)* K=1*JJ) 

1081 WRITE (6*2080) (MF13(I*J)* 1=1*9). (SDF13(I*J)* 1=1*9) 

WRITE (6.2024) EXPA(3> 

WRITE (6*2016) 

DO 108? J = 1 *NT 
JJ = COUNT (J) 

WRITE (6*2045) J. (MOTNO(J.K). (I13(I*J*K>* 1=1*9). K=1.JJ) 

1082 WRITE (6*2080) (MI13(I*J)* 1=1*9)* (SDI13(I*J)* 1=1.9) 

WRITE (6.2024) EXPA(3) 

WRITE (6*2017) 

DO 1083 J = 1*NT 
JJ = COUNT (J) 

WRITE (6.2045) J* (MOTNO(J*K). ( ISPU13 ( I * J *K ) * I = 1.Q). K = 1.JJ) 

1083 WRITE (6*2080) (MIU13 (I . J) * 1 = 1 *9) . (SDIU13 ( I * J) .1 = 1*9) 

IF (NP .LE. 3) GO TO 1099 
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DO 1103 I = NFIRST ,NTIMES 
DO 1103 J r 1,NT 
L = 0 

JJ = COUNT ( J) 

DO 1200 K = 1 , J J 
L = L + 1 

I14< I * J*K) = I VAC 1(1, J,K) - ( EXPA (4 ) * AETOT ( J, K ) *CPHI ) *TTME < I , J , K ) 
F14 ( I » J.K ) = I 14(1, J.K) / TIME ( I » J,K ) 

I SPU14 ( I , J, K ) = I14(I»J,K) / WPIJ(I»J,K) 

SETKL) = 114(1, J,K> 

SET4 t L > = F14 ( I » J»K) 

1200 SET7 (L ) = ISPU14 ( I » J,K ) 

CALL STATS (SET 1 , L, MI14(I,J)» Sr>I14(I,J)) 

CALL STATS (SET4 , L, MF14(I,J), SDF14(I,J>) 

1103 CALL STATS (SET7. L, MIU14(I.J), SUIU14(I»J)) 

WRITE (6*2024) FXPA(4) 

WRITE (6.2015) 

DO 1091 J = 1 . NT 
JJ = COUNT (J) 

WRITE (6.2045) J. (MOTNO(J.K), (H14U.J.K). 1=1.9). K=1.JJ) 

1091 WRITE (6.2060) (MF14(I.J). 1=1.9), (SDF14(I.J). 1=1.9) 

WRITE (6,2024) EXPA(4) 

WRITE (6,201b) 

DO 1092 J = l.NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J.K), (114(1, J.K), 1=1,9), K=1.JJ) 

1092 WRITE (6,2060) (MI14(I,J)» 1=1,9), (bnil4(I»J), 1=1,9) 

WRITE (6,2024) FXPA(4) 

WRITE (6,2017) 

DO 1093 J = 1 ,NT 
JJ = COUNT (J) 

WRITE (6,2045) J, (MOTNO(J,K), ( 1SPU14 ( I , J, K ) , 1=1,9), K=1»JJ) 

1093 WRITE (6,2060) (MIU14 ( I , J) , 1=1 ,9) , ( SPIU14 ( I , J) , 1=1 » 9 ) 

C 

C COMPUTE RU. KM AND THEIR STATISTICS 

C 

1099 DO 1104 J = l.NT 
L = 0 

JJ = COUNT (J) 

DO 1120 K = 1,JJ 
L = L + 1 

RU ( J.K ) = WEb(J.K) / TIME ( NPURN , J ,K ) 

KN ( J , K ) = 6 * I NTP ( Ndl JRN » J » K ) / (UFNS ( J »K ) *CSTU (NRIIKN , J, K ) *RIJ ( J , K ) ) 
SET6 ( L ) = RU ( J,K ) 

SET7 ( L ) = KN (J.K) 

1120 CONTINUE 

CALL STATS (SET6, L, MRU ( J) , SURU ( J)) 

CALL STATS (SET7, L, MKN ( J) , SuKN ( J) ) 

1104 CONTINUE 
WRITE (6,2009) 

DO 1095 J = l.NT 
JJ = COUNT (J) 

WRITE (6,2046) J, (MOTNO(J.K), ( CSTU (I, J.K), I = 1,9). K = 1,JJ) 
1095 WRITE (6,2060) (MCSTU( I , J) , 1 = 1 ,9) , (SOCSTIM I . J) , 1 = 1 . 9) 

WRITE (6,2012) 

DO 1077 J = l.NT 
JJ = COUNT (J) 

IftRITE (6,2046) J. (MOTNO( J.K ) , (WPU(I,J,K)» 1 = 1,9), K = l.JJ) 

1077 WRITE (6,2060) (MWPU(I.J), 1=1,9). (SOWPU(I.J), 1=1,9) 

WRITE (6,2010) 

DO 1096 J = l.NT 
JJ = COUNT (J) 



non 


WRITE <6r2044) Jr (MOTNO(JrK)r RU(JrK)r KN( JrK) r KrlrJJ) 

1098 WRITE <8r2063) MRU(J)r MKN ( J ) r SDRUlJJr SDKN(J) 

PRINTOUT OF THROAT AREAS r EXIT AREAS AND AREA RATIOS 

WRITE (6 >2001 ) 

DO 701 J = lrNT 
JJ = COUNT ( JJ 

701 WRITE ( 8 r 2002 ) Jr (MOTNO(JrK)r ATI ( JrK ) r AT2( JrK) r AT3 ( Jr K ) r 

1 AT4 ( JrK ) r ATTOT ( JrK ) r K = lrJJ) 

WRITE ( 8 r 2070 ) MAT r SDAT 
WRITE (6 r 2071 ) CIMNr CIMXr C2MNr C2MX 
UO 1707 K = IrS 
C1MN(K) = MAE(K) - CP1*SDAE(K) 

Cl MX (K ) = MAE (K ) + CPl*SOAE(K) 

C2MN ( K ) = MAE ( K ) - CP2*SDAE(K) 

1707 C2MX(K) = MAE ( K ) + CP2*S0AF(K) 

WRITE (6r2003) 

DO 70? J = lrNT 
JJ = COUNT (J) 

70? WRITE <8r?002) Jr (MOTNO(JrK)r AEI l JrK ) r AE2(JrK)> AE3( JrK ) r 
1 AF4(JrK) r AETOT (JrK ) t K = lrJJ) 

WRITt (8r?070) MmE> SDAE 

WRITE ( 8 >2071 ) CIMNr CIMXr C2MN> C2MX 

DO 1700 K = IrS 

C1MN(K) = ME < K ) - CP1*SDE (K) 

C1MX(K) = MF (K) + CP1*SDF (K) 

C?MN(K) = ME (K) - CP2*S0E (K) 

1708 C2MX ( K ) = mf (K) + CP?*Sl)E (K) 

WRITE ( 8 r 2004 ) 

DO 703 J = lrNT 
JJ = COUNT (J) 

703 WRITE (8r?002) Jr (MOTNO(JrK)r EA(J,K)> ER(JrK)r EC(JrK)r 
1 ED ( JrK ) r ETOT ( JrK ) r K = lrJJ) 

WRITE ( 8 r 2070 ) ME r SHE 
WRITE < 8 r 2071 ) CIMNr CIMXr C2MNr C2MX 
C 

C COMPUTE CONSTITUENTS OF LNP USING LEAST SQUARES 

C 

LL = O 

DO 28? J = lrNT 
JJ = COUNT (J) 

DO ?8? K = lrJJ 
LL = LL + 1 
TENTRY(LL) = TEMP ( JrK ) 

262 1 ENT? (LL) = TENTRY (LL ) **? 

C 

300 DO 311 I = lrNTIMES 
IFLAGl(I) = 0 
IFLAG2 ( I ) = 0 

IF (P(lrlrD) 308 r 308r 309 

308 IFLAGl(I) = 1 
GO TO 311 

309 LL = 0 

DO 310 J = lrNT 
JJ = COUNT ( J) 

DO 310 K = lrJJ 
LL = LL + 1 

tNTRY(IrLL) = Al OG (P(IrJrK)) 

310 CONTINUE. 

311 CONTINUE 

NM = COUNT ( 1 ) + COUNT (? ) + C0UNT(3) 



non n non 


DO 320 I = 1 r NTIMES 
IF ( IFLA61 ( I ) - 1) 312 * 320 r 312 
-312 DO 315 LL = 1 r NM 
ARRAY ( LL # 1 ) = 1. 

ARRAY (LL» 2 ) = TFNTRY(Lt-) 

ARRAY (LL » 3 ) = TENT2 (LL ) 

315 ARRAY (LLr4) = ENTRY ( I rLL) 

CALL GLS1 ( ARRAY r Xr ILr NMr 3r ALPHA r O.r 0.) 

COEFl (Irl) = X<1) 

COEFl (Ir2> = X (2 ) 

COEFl ( I r 3) = X ( 3) 

320 CONTINUE 

COMPUTE CONSTITUENTS OF LNT USING LEAST SQUARES 

400 DO 411 I = lrNTIMES 

IF ( TIME (Iflrl) ) 408 r 408 * 409 

408 IFLAG? ( I ) = 1 
GO TO 411 

409 LL = 0 

DO 410 J = 1 p NT 
JJ = COUNT (J) 

DO 410 K = lrJJ 
LL = LL + 1 

ENTRY (IrLL) = ALOG ( TIME ( I r J r K ) ) 

41T) CONTINUE 

411 CONTINUE 

DO 420 I = 1 p NT IMES 

IF ( IFLAG2 ( I ) - 1) 412 p 420 * 412 

412 DO 415 LL = lrNM 
ARRAY ( LL r 1 ) = 1. 

ARRAY ( LL r 2 ) = TENTRY(LL) 

ARRAY ( LL p 3 ) = TENT2 ( LL ) 

415 ARRAY ( LL r 4 > □ ENTRY ( I rLL) 

CALL GLS1 ( ARRAY p Xr ILr MMr 3r ALPHAr O.r 0.) 

COEF2 ( I r 1 ) = X(l) 

COEF2 ( I r 2 ) = X ( 2 ) 

COEF2 ( I r 3 ) = X(3) 

420 CONTINUE 


COMPUTE ACTUAL LNP AND LNT 


C 


450 DO 460 I = lrNTIMES 
DO 460 J = 1 r NT 

LNP(IrJ) = COEF 1 ( I r 1 ) + COEFl ( I r 2 > *EXTEM ( J ) + COEFl ( I r 3 ) *EXTEM ( J ) **2 
460 LNT(IrJ) = COEF2<I*l)+ COEF2 ( I r 2 ) *EX1 EM ( J) + COEF2 ( I r 3 > *EXTEM ( J ) **2 


DO 462 I = lrNTIMES 
SETl(l) = 0* 

IF (MP(Irl) .GT. 0.) SETl(l) 

1 

SET 1(2) = 0. 

IF (MP(Ir2) .GT. 0.) SETK2) 

1 

SET1(3) = 0. 


COEF 1 ( I r 1 ) + COEF 1 ( I r 2 ) *MTEMP ( 1 ) + 
COEFl ( X r 3) *MTEMP ( 1 ) **2 

COEF 1 ( I r 1 ) + COEFl ( I r 2 ) *MTEMP fi C2 ) + 
COEFl ( I r 3 ) *MTEMP ( 2 ) **2 


IF (MP ( I r 3) .GT. 0.) SET1<3)= COEFl(Irl) + COEF 1 ( I r 2 ) *MTEMP ( 3 ) + 
1 COEFl ( I r 3 ) *MTEMP ( 3 ) **2 

PIKll(I) = (LNP(Irl) - SET1 ( 1 ) ) / (EXTEMU) - MTEMP ( 1 ) ) 


PIK12 ( I ) 
PIK13 ( I ) 


( LNP ( I r 2 ) 
( LNP ( I r 3 ) 


SET1 ( 1 ) ) / (EXTEM (2 ) 
SET1 ( 1 ) ) / (EXTEM (3) 


MTEMP ( 1 ) ) 
MTEMP ( 1 ) ) 
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u u o o 


•FIK2KX) r (LNP(Irl) - SET1 (2) ) / (EXTFM(l) - MTEMP (£) ) 

PIK22 ( I ) = (LNP (1*2) - SETl(2)) / (FXTEM (2) - MTEMP(2>) 

PIK23 ( I ) = ( LNP (1*3) - SET1 (2) ) / (FXTEM (3) - MTEMP ( 2 ) ) 

PIK3KI) = ( LNP (I * 1 ) - 5ET1 (3) ) / (EXTEM(l) - MTEMP ( 3) ) 

PIK32 ( I ) = ( LNP (1*2) - SET1 ( 3) ) / (EXTEM (2 ) - MTEMP ( 3 ) ) 

PIK33 ( I ) = ( LNP (1*3) - SET1 (3) ) / (EXTEM (3) - MTEMP<3>> 

SET2(1) = 0. 

IF (MTIME(Ifl) .GT. 0.) SET2(1)= COEF2(I*l) +COEF2 (I * 2 ) *MTEMP ( 1 ) + 


I 

SET2 ( 2 ) = 0. 

IF (MTIME ( I * 2 ) .GT. 0.) SET2(2) = 

1 

SET2 ( 3 ) = 0. 

IF (MTIME (1*3) .GT. 0.) SET2(3)= 

1 


SIGP11 ( I 
SIGP12 ( I 
S IGP13 ( I 
' SIGP21 ( I 
SIGP22 ( I 
SIGP23 ( I 
SIGP3I ( I 
5 IGP32 ( 1 
462 SIGP33 ( I 
WRITE (6 


) = ( SET2 ( 1 ) 
) = ( SFT2 ( 1 ) 
) = ( SF T2 ( 1 ) 
) = (S^TRfP) 
) = ( SRT2 ( 2 ) 
) = <S^T?(2) 
) = ( SFT2 ( 3 ) 
) r (SfT2(3) 
) = ( SET 2(3) 
* 2050 ) 


LNTdrl) ) 
LNT ( J * 2 ) ) 
LNT (1*3) ) 
LNT ( I r 1 > ) 
LNT (1*2) ) 
LNT (1*3) ) 
LNT (1*1)) 
LNT (I *2) ) 
LNT( I *3) ) 


COFF 2(1*3) *MTtMP ( 1 ) **2 

COFF2 (1*1) +C0EF2(I*2)*MTEMP(2) + 
COFF 2 ( I * 3 ) *MTEMP ( 2 ) **2 


C0EF2 (1*1) +COEF2 (1*2) +MTEMP ( 3 ) 
COFh 2(1*3) *MTEMP ( 3 ) **2 


(tXTEM(l) 
( LXTFM ( 2 ) 
( lXTEM ( 3 ) 
(tXTFM(l) 
(EXTEM (2) 
(EXTEM ( 3 ) 
(tXTFM(l) 
(EXTEM ( 2 > 
( c-XTEM ( 3 ) 


MTEMP ( 1 ) ) 
MTFMP U ) ) 
MTFMP ( 1 ) ) 
MTFMP ( 2 ) ) 
MTEMP (2 ) ) 
MTFMP (2) ) 
MTFMP (3) ) 
MTEMP (3) ) 
MTFMP ( 3 ) ) 


WRITE (6*2051) PIK1J * PIK21* PIK31* PIK12* PIK22* PIK32* PIK13* 
1 PIK23 * PIK33 

toR I TE ( 6 * 2052 ) SI GP1 1 * SI GP21 *SI GP31 *SIGP12 * SIGP2? * SIGP32 * SIGP13 * 
1 S IGP23 * S IGP33 


COMPUTE P TT /'NO TIME TT TRmNSFORMAT IONS - SELECTOR 
CHOOSES WHCIH TEMPFRATURF DATA IS TO bF TRANSFORMED TO 


11 = 1 

12 = ? 

13 = 3 

480 GO ro (500* 520* 540)*TFMSEL 

500 DO 510 I = 1 * NT I MFS 
JJ = COUNT (1) 

DO 502 K = 1*J1 

IF ( IFLAG1 ( I ) .EG. 1) GO TO 501 

PI ( I * 1 * K ) = P ( I * 1 * K ) * EXP (PIKll ( I)*(EXTFM(1) —TEMP ( 1 * K ) ) ) 

501 IF ( IFLAG2 ( I ) .EG. 1) Go TO 502 

T1(I*1*K) = TIME ( I * 1 * K ) / (FXP ( SlGPl 1 ( I ) * ( FXTEM ( 1 ) - TEMP(1*K)>)> 

502 CONTINUE 

DO 504 K =1 *J2 

IF (IFLAGl(I) .EG. 1) GO TO 503 

P2(I*1*K) = P(I*2*K) * EXP (PIK21 ( I)*(EXTEM(1) -TEMP ( 2 * K ) ) ) 

503 IF (IFL/'G2(I) .EG. 1) GO TO 504 

T2(I*1*K) = TIME ( I * 2 * K ) / (FXP ( SIGP21 ( I ) * ( EXTEM ( 1 ) - TEMP(2*K)))) 

504 CONTINUE 

DO 510 K =1*J3 

IF (IELAGl(I) .EQ. 1) GO TO 508 

P3(I*)*K) r P ( T * 3 * K ) * EXP ( PIK31 ( I ) * (EXTEM ( 1 ) -TEMP ( 3 * K ) ) ) 

508 IF ( IFLAG2 ( I ) .EG. 1) GO TO 510 

T3(I*1*K> = T IMF ( I * 3 * K ) / (EXP (SIGP31 < I) * (EXTEM < 1 >- TEMP(3*K)))) 
510 CONTINUE 
JSAVE - 1 
GO TO 580 

520 DO 530 I = IrMTIMES 
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JJ = COUNT (2) 

DO 522 K = ltJl 

IF ( IFLAGl (I) .EQ. 1) GO TO 521 

PI ( I * 2 »K ) = P(ItlfK) * EXP (PIK12 ( I ) * (EXTEM (2 )— TEMP (1 »K) ) ) 

521 IF ( IFLAG2 ( I ) *EQ. 1) GO TO 522 

T1 ( I #2* K ) = TIME C I # 1 #K ) / (EXP (SIGP12 < I ) * (EXTEM (2) - TEMP(1#K)>>> 

522 CONTINUE 

DO 524 K = lfJ2 

IF ( IFLAGl ( I ) .EQ. 1) GO TO 523 

P2 ( I * 2 * K ) = P(I»2»K) * EXP (PIK22 ( I ) * ( EXTEM ( 2 )— TEMP (2»K ) ) ) 

523 IF ( IFLAG2 ( I ) .EQ. 1) GO TO 524 

T2(I.2#K) = TIME (1 f 2 f K > / (EXP (SIGP22 ( I > * (EXTEM (2)- TEMP(2»K)))) 

524 CONTINUE 

DO 530 K =1#J3 

IF (IFLAGl (I) .EQ. 1> GO TO 528 

P3(I»2rK) = P(I*3fK> * EXP (PIK32(I)*(EXTEM(2)-TEMP(3»K) ) ) 

528 IF ( IFLAG2 ( I ) .EQ. 1) GO TO 530 

T3 ( I r 2 » K ) = TIME(If3fK) / (EXP (SIGP32( I ) * (EXTEM (2 )- TEMP(3»K)))) 
530 CONTINUE 
JSAVE = 2 
GO TO 580 

540 DO 550 I = lrNTIMES 
JJ □ COUNT (3) 

DO 542 K = 1#J1 

IF ( IFLAGl ( I ) .EQ. 1) GO TO 541 

PI ( I * 3 r K ) = P(IflrK) * EXP ( P IK 13 ( I ) * ( EXTFM ( 3 ) -TEMP ( 1 ♦ K ) ) ) 

541 IF ( IFLAG2 ( I ) .EQ. 1 ) GO TO 542 

Tl(I«3rK) = TIME ( I f 1 r K ) / (FXP ( SIGP1 3 ( I ) * ( EXTEM ( 3 ) - TEMP ( 1 r K ) ) ) ) 

542 CONTINUE 

DO 544 K -1 » J2 

IF ( IFLAGl ( I ) .EQ. 1) GO TO 543 

P2(Ir3#K) = P(I.2*K> * EXP ( PIK23 ( I ) * ( FXTEM ( 3 ) —TEMP ( 2 * K ) ) ) 

543 IF ( IFLAG2 ( I ) .EQ. 1) GO TO 544 

T2(I#3*K> = TIME(X»2’K) / (FXP (SIGP23 ( I ) * ( EXTEM ( 3) - TEMP (2 f K ) ) ) ) 

544 CONTINUE 

DO 550 K “If J3 

IF ( IFLAGl ( I ) .EQ. 1> GO TO 548 

P3 ( I f 3 r K ) = P(I.3fK) * EXP ( PIK33 (I)*(EXTEM(3) —TEMP ( 3 r K ) ) ) 

548 IF ( IFLAG2 ( I ) .EQ. 1) GO TO 550 

T3 ( I » 3 ► K ) = TlME(Ir3#K) / (FXP (SIGP33 ( I ) * (EXTEM (3)— TEMP(3»K) ) ) ) 
550 CONTINUE 
JSAVE = 3 
GO TO 580 

COMPUTE STATISTICS FOR PTT AND TIME TT 


580 J = JSAVE 

DO 640 I = 1 9 NTIMES 
L = 0 
M = 0 

DO 620 K = ltJl 
L = L + 1 
M = M + 1 

TENTRY (L ) = PKIrJrK) 
TENT2 (L ) = Tl(IfJfK) 
SET1(M ) = Pl(IrJ»K) 
620 SET4 ( M ) = Tl(IrJrK) 

M = 0 

DO 625 K = If J2 
L = L + 1 
M = M + 1 
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TENTRY (L) = P2(IfJfK) 
TENT2(L) = T2(IfJfK) 
SET2(M) = P2 (IfJfK) 
625 SET5 ( M ) .= T2(IfJ#K) 
M = 0 
DO 630 K 
L = L + 1 
M = M + 1 
TENTRY (L ) = 

TENT2CL) = 

SET3 ( M) = 

630 SET6 CM) = 

CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
640 CALL STATS 
IF (TEMSEL 
TEMSEL = TEMSEL 
GO TO 480 


= lr J3 


P3 (I » J»K ) 

T3 (IfJfK) 

P3 (IfJfK) 
T3(Ir JrK) 

( SET 1 w J1 
(SET2 f 
(SET3f 
( SET4 * 

( SETS r 
(SET6# 

( TENTRY f L f 
( TENT2 f L f 
.GE. NT) GO 
+ 1 


J2r 
J3 r 
Jlf 
J2f 
J3f 


MPT1(IfJ)f 
MPT2<I»J) * 
MPT3CX»J ) -f 
MTTKIfJ) f 
MTT2CI# J) f 
KTT 3 ( I f J ) f 
MPTT ( I r J) r 
MTTT(IfJ ) f 
TO 6^5 


SuPTl (If J) ) 
SuPT? ( I f J) ) 
SUPT3 ( I r J) ) 
SOTT1 ( I f J) ) 
SUTT2 ( I f J ) ) 
Sl>TT 3 ( I r J) ) 
SuPTT ( I f J) ) 
SdTTT(IfJ) ) 


COMPUTE INTEGRALS OF P TRANSFORMED 

655 DO 700 I = 1 fNTIMES 
LO 700 J = 1 f NT 
L. = 0 
M = 0 

DO 660 K = lfJl 

1NTPT1 (IfJfK) = PKIfJfK) * TIUfJfK) 

L = L + 1 
M = M + 1 

TENTRY ( L ) = INTPT1 (IfJfK) 

660 SET 1 (M ) = INTPTKIfJfK) 

M = 0 

DO 670 K = 1fJ2 
L = L + 1 
M = M + 1 

INTPT? (IfJfK) = P2 ( I f J f K ) * T2(IfJfK) 

TEN TRY (L ) = INTPT2 (IfJfK) 

670 SET2 ( M.) = INTPT2 (IfJfK) 

M = 0 

DO 680 K = If J3 
L = L + 1 
M = M + 1 

INTPT3 (IfJfK) = P3 ( I f J f K ) * T3(IfJfK) 

TENTRY (L ) = INTPT3 ( IfJfK) 

680 SET3(M) = INTPT3 (IfJfK) 

COMPUTE STATS ON INTEGRALS OF P TRANSf-ORMtiD 

CALL STATS (SETIf J1f MINT1(IfJ)f SDINTIUfJ)) 

CALL STATS (SET?f J2f MINT2(IfJ)f SDINTPUfJ)) 

CALL STATS (SET3f J3f MINT3(IfJ)f SPif 'T3 ( I f J ) ) 

700 CALL STATS (TENTRYfLf MINTT(IfJ)f sniNTT(IfJ)) 

675 DO 565 I = IfNTIMES 
WRITE ( 6 f2030 ) I 
WRITE ( 6 f 2047 ) 

WRITE ( 6 f 2040 ) 1 1 f (MOTNO(IfK)f T1(1f1fK)f T1(Tf2fK)f T1(If3fK)f 
1 BLANK r K = IfJJ.) 
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WRITE (6*2064) (MTT1(I*J>* J=1.4)* (SDTTld.J). J=l*4> 

IF (NT .LE. 1) 60 TO 661 

WRITE (6*2040)12. (MOTNO(2,K)» T2d»l»K), T2(I*2*K). T2d,3,K), 
1 BLANK . K = 1 * J2) 

WRITE (6.2064) ( MTT2 ( I * J) * J=l,4), (&nTT2(I.J). J=l*4) 

IF (NT .LE. 2) GO TO 561 

WRITE (6.2040)13* (MOTNO(3»K). T3(1»1*K>. T3(I,2,K). T3(I.3*K). 
1 RLANK . K = 1.J3) 

WRITE (6.2064) (MTT3(I.J>. J=l,4), (SDTT3d,J), J=1.4) 

561 WRITE (6.2065) (MTTTd.J). J=1.4). (SDTTTd.J). J=1.4) 

00 562 J = 1 .NT 

CIMIN(J) = MTTTd.J) - CT1*SDTTT ( I . J) 

CIMAX(J) = MTTTd.J) + CT1*SDTTT ( I . J) 

C2MIN ( J ) = MTTTd.J) - CT2*SDTTT ( I . J) 

562 C2MAX ( J ) = MTTTd.J) + CT2*SDTTT ( I . J) 

WRITE (6.2062) C1MIN. C1MAX. C2MIN, C2MAX 

565 CONTINUE 

DO 570 I = l.NTIMES 
WRITE (6.2031) I 
WRITE (6.2042) 


WRITE (6.2048) 

WRITE (6.2041)11 . (MOTNO(l.K). INTPTl ( 1 . 1 . K ) . INTPT1 ( 1 . 2 .K ) . 

1 INTPT 1 ( 1 . 3 . K ) . BLANK. Pl(I.l.K). PKI.2.K). PKI.3.K). 

2 BLANK. K = 1.J1) 

WRITE (6.2066) (MINTld.J*. J=l,4). (MPTl(I.J). J=l,4).' 

1 (SO INTI ( I . J) . J =1.4). (SOPTld.J). J= 1.4) 

IF (NT .LE. 1) GO TO 564 

WRITE (6.2041)12. <MOTNO(2.K). INTPT2 ( 1 . 1 . K ) . 1NTPT2 ( 1 . 2 . K ) * 

1 INTPT2 ( I . 3 . K ) . BLANK. P2(I.1.K). P2(I.2.K). P2d.3.K). 

2 BLANK. K = 1.J2) 

WRITE (6.2066) (MINT2d.J). J =1.4). (MPT2(I.J)» J=1.4). 

1 (S 0 INT 2 ( I . J ) . J =1,4). (SnPT2(I,J), J=l,4) 

IF (NT . LE • 2) GO TO 564 

WRITE (6,2041)13, (M0TN0(3,K), INTPT3 ( 1 . 1 ,K ) » 1NTPT3 ( T , 2 »K ) , 

1 INTPT3 ( I » 3 ,K ) . BLANK, P3(I,1,K). P3(I.2,K). P3d,3,K), 

2 BLANK, K = 1.J3) 

WRITE (6,2066) ( MI NT 3 ( I » J ) , J =1,4), <MPT3(I»J), J=1,A), 


( SOINT3 ( I » J ) . J =1,4) , (SOPT3(I,J) » 
( MI NTT ( I , J ) , J =1,4), (MPTT(I.J), 
(SOINTT(I.J), J =1,4) . (SOPTTd.J) . 


1 

564 WRITE ( 6 . 2067 ) 

1 

00 566 J= 1 ,NT 

CIMIN(J) = MINTT(l.J) 
CIMAX(J) B MINTT(I.J) 
C 2 M I N ( J ) = Ml NTT ( I , J) 
C 2 MAX ( J ) = MINTT(I.J) 
LlTMN(J) = MPTT (I.J) 
CITMX(J) = MPTT (I.J) 
C 2 TMN ( J ) = MPTT (I.J) 
566 C 2 TMX ( J ) = MPTT (I.J) 
WRITE ( 6 , 2068 ) C 1 MIN, 
1 C 2 MIN , 

570 CONTINUE 
RETURN 
END 


- CT 1 *SDINTT( I , J) 

+ CTl*SniMTT(I, J) 

- CT 2 *SniNTT(I.J) 

+ CT 2 *SniNfT (I.J) 

- CT 1 *S 0 PTT d.J) 

+ CT 1 *S 0 PTT (I.J) 

- CT 2 *S 0 PTT (I.J) 

+ CT 2 *SDPTT (I.J) 
C 1 TMN. Cl MAX , C 1 TMX, 
C 2 TMM, C 2 MAX. C 2 TMX 


J=1 ,4) 
J=1 , 4 ) , 
J=1 , 4 ) 
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SltiFIC GLS1 



5U5HOIJTIMF GL51 (ArXrILiNrMf ALPHA *F1 * t? ) 

Gl-b^onf’l 


LIMbNSIOu A(46’4)» X(4>* IL<4) 

toLA^Ono? 


KM=W+1 

bLSOono^ 


LL=1 

OLSOOOOti 


DO 60J=lrMM 

faLS^0 n 0 r ' 

f>a 

IL(J>=0 

oLbO:infiA 


1 = 1 

Gl-b^(jn t 7 


DO 3K=1 t MM 

GLbG,)' 1 (’f’ 


11 = 1 + 1 

GLSGOOfia 


DO 4 J=I I #N 

GLS00010 


IF (APS(A(J#K) ) — £1 ) 4 f 4 r6 

GLSQ001 1 

6 

T1=SQRT( (A < J#K) ) **2+ ( A ( I » K ) )**?.) 

GLSOoni? 


S=A(JrK)/Tl 

GLSQ0013 


C=A(I#K)/T1 

GLSG0014 


DO 5L=K r MM 

GLSQ0015 


T 2=C*A ( I # L ) +S*A ( J r L ) 

GLSQ0016 


A ( J r L ) =— S*A ( I r L ) +C^A ( J f L ) 

GLSG0017 

5 

A ( I f L ) =T2 

GLSG001R 


LL=LL+1 

GLSQ0019 

4 

CONTINUE 

GLSG0020 


IF (AP'SCACI'K) )-t£2)3#3r8 

GLSQ0021 

8 

IL(K)=I 

GLSQ0022 


1 = 1 + 1 

GLSQ0023 

3 

CONTINUE 

GLSQ0024 


X (MM) =-1.0 

GLS0002S 


II=M 

GLS0002G 


LO 351= 1 »M 

GLSG0027 

35 

X t 1 )=n. 

GLSQ002P 


1 0 30 J=1 9 M 

GLS00n29 


IF (IL(II) )30*30»31 

GLSG0030 

31 

5=0. 

GLSG0031 


LL=I 1+ 1 

GLS00032 


I=ILCIT) 

GLS00033 


DO 32K=LL>MNi 

GLSG0034 

3? 

5 =S + A ( I r K ) *X ( K ) 

GLS0003B 


X(II>=-S/A<1, II) 

GLSQ0036 

30 

11=11-1 

GLSO0037 


IF (IL(MM) ) 50 f 51 > 50 

GLSQ003R 

51 

ALPHA=0 • 

GLS00039 


GO TO 52 

GLSG0040 

50 

I =IL (MM ) 

bLS<30n41 


ALPHA=A(I»MM) 

GLS00042 

5? 

return 

GL500n43 


END 

GLSQ0044 
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ALPHA 


SORIGIN 
SIBFTC M0T2 

SUBROUTINE M0T2 

SECOND LINK OF MOTORS PROGRAM 

DIMENSION XX (4860 ) » 2Z< 2124) 

COMMON G. CPHI. NT. EXPA. COUNT. WP» WEB. ATTOT. A£TOT» 

1 P. MOTNO. FVAC1. NTIMES. XX. NTOT . NBURN. NFIRST » NP 

COMMON CPI. CP2. CT1. CT2. VEC1. VEC2. VEC3 

COMMON SET! . SET2. SET3. SET4. SET5. SET6. SET7-* C1MIN. CIMAX. 

1 C2MIN. C2MAX. C1TMN. C1TMX. C2TMN. C2TMX 

DIMENSION SAVE1 (9*4.15) * SAVE2 (9.4 . 15) * SAVE3 (9*4* 15) 

DIMENSION MFT11 (9.4). S0FT1K9.4). MIT11 (9.4). SGIT1K9.4). 

1 MFT21 (9.4). SDFT21 (9.4) . MIT21 (9.4) . SDIT21 (9.4) . MFT31 (9.4). 

2 S0FT31 (9.4) . MIT31 (9.4). S0IT3K9.4). MFTT1 (9.4). SUFTT1 (9.4) » 

3 MITT1 (9.4). S0ITTK9.4). MFT12 (9.4). S0FT12(9.4). MIT12 (9.4). 

4 SDIT12 (9.4) . MFT22 (9.4). SDFT22<9»4). MIT22 (9.4). SDIT22<9»4>» 

5 MFT32 (9.4). SDFT32(9.4). MIT32 (9.4). S0IT32(9.4). MFTT2 (9.4). 

6 SOFTT? ( 9.4 ) . MITT2 (9.4). SDITT2(9.4) 

DIMENSION MIS11 (9.4). SDIS1K9.4). MIS21 (8.4). SDIS2K9.4). 

4 MIS31 (9.4). SDIS31 (9.4) . MISTI (9.4). SDIST1 (9.4) » MIS12 (9.4). 

5 SDIS12 (9.4) . MIS22 (9.4). SDIS22 (9 .4 > . MIS32 (9,4). SDIS32<9»4)» 

6 MIST2 (9.4). SDIST2 (9.4) 

DIMENSION 

1 FT11 (9.4.15). FT21 (9.4.15). FT31 (9.4.15). FT12 (9.4.15). 

2 FT22 (9.4.15). FT32 (9.4.15). 

4 IT11 (9.4.15). IT21 (9.4.15). IT31 (9.4.15). IT12 (9.4.15). 

5 IT22 (9.4. 5). IT32 (9.4.15). 

7 ISP11 (9.4.15) . ISP21 (9.4.15) . ISP31 (9 .4. 15) . ISP12 (9.4. 15) . 

8 ISP22 (9.4.15) . ISP32 (9.4.15) 

DIMENSION RATI0K9.4.15) . RATI02 (9 .4 . 15 ) . RATI03(9.4. 15) 

DIMENSION WP1 (9.4.15). WP2 (9.4.15). WP3 (9.4.15) 

VARIABLES OBTAINED FROM CHAIN 1 

DIMENSION PI (9.4.15). P2 (9.4.15). P3 (9.4.15). 

1 EXPA ( 4 ) * P (9.4.15). FVAC1 (9.4.15). 

2 AETOT (4,15) . INTPT1 ( 9 . 4 . 15 ) * INTPT2 ( 9 . 4 . 15 ) . 

3 ATTOT (4.15) » INTPT3 (9.4.15) » T1 (9.4.15). T2 (9.4.15). 

4 T3 (9,4.15). WP (4,15). VECK4.15). VEC2 (4.15), VFC3(4,15> 

DIMENSION C1MIN(4), C1MAX(4)» C2KIN(4>» C2MAX(4) 

DIMENSION SET1 (45). SFT2 (45), SFT3 (45). SET4 (45). SFT5 (45), 

1 SET6 (45 ). TENT2 (60 ) , MWPK9.4), MwP2(9,4>. MWP3(9,4). 

2 snwpl(9,4), SPWP2 ( 9,4 ) » SDWP3 (9.4) » TENTRY ( 60 ) , 

3 MOTNO (4,15) , C0UNT(4). MWPT(9,4)» SDWPT(9,4)» WER(4,15) 

REAL MFT11 , MIT11 » MFT21 » MIT21 , MFT31 » MIT31 » MFTT1 , 

1 MITT1 , MFT12 » MIT12 » MFT22 » MIT22 » MFT32 , MIT3? , MFTT? » 

2 IT11 , MOTNO , MWP1 , MWP2 » MWP3 , MWPT » 

4 MITT? , MIS11 . MIS21 , MIS31 » MISTI , MIS12 » MIS2? » MIS3? * 

7 MIST2 , IT21 » IT31 * IT12 » IT22 » IT32 , ISP11 » 

9 ISP21 . ISP31 , ISP12 » ISP22 » ISP32 , INTPT1, INTPT2. INTPT3 
INTEGER COUNT 

EQUIVALENCE 

1 (XX( 1) » Pl(l) » WP1 ( 1) ) » 

2 (XX( 541). P2 ( 1 ) » WP2 ( 1 ) ) » 

3 (XX(lOftl), P3 ( 1 ) . WP3 ( 1 ) ) » 

4 (XX(1621) , . Till) . IT 11(1). IT12 ( 1 ) » ISPlKl), ISP12(1)). 


i 
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5 (XX (2161 ) » T2(l>. IT21U). IT22(1)» ISP2K1). ISP22(1)). 

6 ( XX (2701 ) » T3 ( 1 ) » IT31 ( 1 ) ► IT32(1>. ISP31 (1) . ISP32(1))» 

7 ( XX ( 3241 ) » INTPT1 ( 1 ) ) » 

8 (XX (3781 >» INTPT2 ( 1 ) ) * 

9 (XX (4321 ) * INTPT3 ( 1 ) ) 

EQUIVALENCE 

1 (ZZ( l)r MFT1K1). MFT12(1)» MFT13(1)» MFT14(1). MWP1 (1)>» 

2 (ZZ( 37 ) * MFT21 ( 1 ) » MFT22(1)» MFT23(1)» MFT24(1>» MWP2 (1))» 

3 (ZZ( 73) r MFT31 ( 1 ) » MFT32 ( 1 ) » MFT33()>» MFT34 ( 1 ) » MWP3 (1))» 

4 (ZZ< 109). MITll(l). MIT12 ( 1 ) . MlSll(l). MIS12(1)). 

5 (ZZ( 145). MIT2K1). MIT22(1). MIS2K1). MIS22(1))» 

6 (ZZ( 181). MIT3K1). MIT32(1). MIS3K1). MIS32(1)) 

EQUIVALENCE 

1 (ZZ( 217). SOFT 1 1 ( 1 ) » SDFT12 ( 1 ) .5UFT13 ( 1 > .SDFT14 ( 1 ) » SDWP1 (1)). 

2 (ZZ( 253 ) » SDFT21 ( 1 ) .SDFT22 ( 1 ) .5UFT23 ( 1 ) .SDFT24 ( 1 ) » SOWP2 (D). 

3 (ZZ( 289). SOFT 31 ( 1 ) .SOFT32 ( 1 ) .SUFT33 (1) .SOFT 34(1) » SDWP3 (1)). 

4 <ZZ( 325). SOIT11 ( 1 ) . SniT12(l)» SDISll(l). SUIS12(1)). 

5 (ZZ( 361). SDIT21 ( 1 ) . SDIT22(1). SDIS21 ( 1 ) . SDIS22(1))» 

6 (ZZ( 397). S0IT31 ( 1 ) * SDIT32(1). S0IS3K1). SDIS32(1)) 
EQUIVALENCE 

1 (ZZ( 433). MFTTl(i). MFTT2 ( 1 ) » PFTT3(1). MFTT4(1))» 

2 (ZZ( 469 ) » SDFTT1 ( 1 ) . 5DFTT2 ( 1 ) .SDFTT3 ( 1 ) .SDFTT4 ( 1 ) ) . 

3 (ZZ( 505). FTll(l). FT12 ( 1 ) . FT13(1). FT14(1)). 

4 (ZZ ( 1045) . FT21(1). FT22(1)» FT23(1). FT24 ( 1 ) ) . 

5 (ZZ ( 1585) . FT31 ( 1 ) » FT32(1)» FT33(1)» FT34 ( 1) ) 

C 

C OUTPUT FORMATS 

C 

2014 FORMAT ( 34X. 5H PA = F5.2) 

2032 FORMAT (1H1. 30X. 20H TRANSFORMED THRUST / 9H TIME NO. 12) 

2033 FORMAT (1H1. 30X. 20H TRANSFORMED IMPULSE/ 9H TIME NO. 12) 

2034 FORMAT (1H1. 30X. 30H TRANSFORMED SPECIFIC IMPULSE /9H TIME NO. 12) 

2037 FORMAT (1H1. 30X. 28H TRANSFORMED PROPFLLANT WT. / 9H TIMF MO. 12) 

2041 FORMAT (11H0TEMP GROUP 12 // (IX. A6. 19X. 4F12.2)) 

2047 FORMAT ( 30X. 43H TEMPI TEMP2 TFMP3 TEMP4 / 

1 10H MOTOR NO.) 

2062 FORMAT ( 34HOCONF IOENCF. ON NORMAL DISTRIBUTION / 


1 

26H (ONE 

MIN 

4F12.2 

/ 

2 

26H sineu 

MAX 

4F12.2 

// 

3 

26H (TWO 

MIN 

4F12.2 

/ 

4 

26H SIHEL) 

MAX 

4F12.2 

) 

2064 FORMAT 

(26H0MEAN 


4F12.2 

/ 

1 

26H STANDARD DEV. 

4F12.2 

) 

2065 FORMAT 

(26HOTOTAL 

MEAN 

4F12.2 

/ 

1 

26H TOTAL 

STANDARD DEV. 

4F 12 • 2 

) 


C 

11 = 1 
12 = 2 
13 = 3 
BLANK = 0. 

J1 = COUNT ( 1 ) 

J2 = COUNT (2) 

J3 = COUNT (3) 

DO 696 I = 1.504 
6°6 ZZ(I) = 0. 

C1MIN (4) = 0. 

C1MAX(4) = 0. 

C2MIN (4) = 0. 

C2MAX (4 ) = 0. 

DO 701 J = 1 » NT 
DO 697 K = 1.J1 
fa97 VECl(J.K) □ T 1 (N3URN . J.K ) 
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DO 698 K = 1#J2 

698 VEC2 ( J#K ) = T2 (N6URN# J#K) 
DO 699 K = 1#J3 

699 VEC3 ( J#K ) = T3 (N6URN# J#K > 
701 CONTINUE 


COMPUTE THRUST # IMPULSE AND THEIR STATISTICS FOR PI 


DO 752 I = 1#NTIMES 
DO 752 J = 1 t NT 
L = 0 
M = 0 

DO 705 K = 1#J1 
L = L + 1 
M = M + 1 

RATIOKI# J'K) = PUIpJ'K) / P (I#1#K) 

FT11(I#J#K> = RATIQ1(I#J#K) *FVAC1 ( I # 1 #K ) -EXPA ( 1 ) *AtTOT ( 1 #K ) *CPHI 
IT11(I#J#K) = FT11<I#J#K> * T1(I#J#K) 

TENTRY ( L ) = FT11(I#J#K) 

TENT2 (L ) = IT11(I#J#K) 

SET 1 (M) = FT 11(I#J#K) 

705 SET4 (M) = IT11(I#J#K) 

M - 0 

DO 708 K = 1#J2 
L = L + 1 
M = M + 1 

RAT 102 ( I # J # K ) = P2<I#J#K> / P (I#2#K) 

FT21 ( I # J # K ) = RATI02(I#J#K)*FVAC1(I#2#K)-EXPA(1>*AET0T(2#K)*CPHI 
IT21 ( I # J # K ) = FT21 ( I # J # K ) * T2(I#J#K) 

TENTRY (L ) = FT2l(I#J#K> 

TENT2 (L ) = IT21(I#J#K) 

SET2 (M) = FT21 ( I # J # K ) 

708 SET5 (M) = IT21(I#J#K) 

M = 0 

DO 710 K = 1 r J3 
L = L + 1 
M = M + 1 

RATI 03 (I # J # K ) = P3(I#J#K) / P (I#3#K) 

FT31 ( I # J #K ) = RATI03(I#J#K) +FVAC1 ( 1 # 3 #K ) -EXPA ( 1 ) *AETOT ( 3 #K ) *CPHI 
IT31 ( I # J #K ) = FT31 ( I # J # K ) * T3(I#J#K> 


710 


752 


TENTRY (L ) = 
TENT2 ( L ) r 
SET3 (M) - 
SET6 CM) = 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CALL STATS 
CONTINUE 


FT31 ( I # J # K 
IT31 ( I # J # K 
FT31(I#J#K 
IT31 ( I # J#K 
(SET1 # Jl# 
( SET4 # J1 # 
(SET2 # J2 # 
(SETS # J2# 
(SET3 # J3# 
(SETS# J3# 
(TENTRY #L# 
(TENT2 #L# 


MFT11 ( I # J) # 
MIT1 1 ( I # J) # 
MFT21 ( I # J ) # 
MIT21 ( I # J ) # 
MFT31 ( I # J) # 
MIT31(IrJ) # 
MFTT1 ( I # J) # 
MITTKIpJ) # 


SDFT11 ( I # J) ) 
SDITIKIpJ) ) 
SDFT21 ( I # J) ) 
SDIT21 ( I # J ) ) 
SDFT31 ( I # J ) ) 
S0IT31 ( I # J) ) 
SDFTTKIpJ) ) 
SDITTKI# J) ) 


PRINT THRUST# IMPULSE AND STATS FOR PI 

DO 1710 I = 1 # NTIMES 
WRITE ( 6 #2032 ) I 
WRITE ( 6 #2014 ) EXPA(l) 

WRITE ( 6 # 2047 ) 

WRITE ( 6# 2041 ) II# (MOTNO ( II #K ) # FT11(I#1#K># FT11(I#P#K># 
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n n n n 


1 FT11 ( I »3»K ) » BLANK . K=1*J1) 

WRITE (6.2064) (MFTll(I.J). J = 1.4). (SUFT11 ( I * J) . d = 1.4) 

IF (NT .LE. 1) GO TO 1708 

WRITE (6.2041) 12. (MOTNO ( I2*K ) * FT21(I»1»K). FT21(I.?*K). 

1 FT21 ( I »3»K) . BLANK. K=1»J2) 

WRITE (6.2064) (MFT2KI.J). J = 1*4). (SuFT21 ( I . J) . J = 1.4) 

IF (NT .LE. 2) GO TO 1708 

WRITE (6.2041) 13. (MOTNO ( 13. K ) . FT3KI.1.K). FT3KI.2.K). 

1 FT31 (I.3.K) » BLANK. K=1.J3) 

WRITE (6.2064) (MFT3KI.J). J = 1.4). (S0FT31 ( I ► J) . J = 1.4) 

1708 WRITE (6.206b) (MFTTl(I.J). J = 1.4). (SDFTT1 ( I » J) . J = 1.4) 

00 1709 J = 1 .NT 

CIMIN(J) = MFTTl(I.J) - CT1*SDFTT1 ( I » J ) 

CIMAX(J) = MFTTl(I.J) + CT1+SDFTT1 ( I » J) 

C2MIN ( J) = MFTTl(I.J) - CT2»SDFTT1 ( I . J) 

1709 C2MAX ( J ) = MFTTl(I.J) + CT2*SDFTT1 ( I. J) 

WRITE (6.2062) C1MIN. C1MAX. C2MIN. C2MAX 

1710 CONTINUE 

DO 1712 I = l.NTIMES 
WRITE (6.2033) I 
WRITE (6.2014) EXPA(l) 

WRITE (6.2047) 

WRITE (6.2041) II. (MOTNO ( 1 1 .K ) » ITll(I.l.K). IT1KI.2.K). 

1 IT1KI.3.K). PLANK. K=1.J1) 

WRITE (6.2064) (MITll(I.J). J = 1.4). (SUlTll ( I . J) . J = 1.4) 

IF (NT .LF. 1) GO TO 1707 

WRITE (6.2041) 12. ( MOTNO < 12 .K ) . IT2KI.1.K). 1T2KI.2.K). 

1 IT2KI.3.K). BLANK. K=1»J2) 

WRITE (6.2064) (MIT2KI.J). J = 1.4). (SUTT21 ( I » J) » J = 1.4) 

IF (NT .LE. 2) GO TO 1707 

WRITE (6.2041) 13. (MOTNOt I3.K ) . IT3KI.1.K). 1T3KI.2.K). 

1 1T3KI.3.K). BLANK. K=1»J3> 

WRITE (6.2064) (MIT31 ( I . J) * J = 1.4). (5UIT31 ( I . J) . J = 1.4) 

1707 WRITE (6.2065) (MITTl(I.J). J = 1.4). (SoITTl ( I . J) . J = 1.4) 

DO 1711 J = 1 . NT 

CIMIN(J) = MITTl(I.J) - CT1*S0ITT1 ( I . J) 

CIMAX(J) = MITT) (I. J) + CT1*SDITT1 (I.J) 

C2MIN ( J) = MITT1 (I.J) - CT2+SDITT1 (I.J) 

1711 C2MAXIJ) = MITT1 (I.J) + CT2*SDITT1 ( I . J) 

WRITt (6.2062) C1MIN. C1MAX. C2MIN, C2MAX 

1712 CONTINUF 

COMPUTE TRANSFORMER PROPFLLANT WFIGHT. SPECIFIC IMPULSE 
FOR PI AND THEIR STATISTICS 

DO 882 I = l.NTIMES 
DO 88? U = 1 .NT 
L = 0 

M = 0 

DO 804 K = 1.J1 
L = L + 1 

M = M + 1 

WPl(I.J.K) = WP(l.K) * INTPT1 ( I . J.K ) / IimTPTI (NTOT. J.K ) 
ISPll(I.J.K) = ITll(I.J.K) / WPl(I.J.K) 

SAVE1 ( I . J.K ) = ISPll(I.J.K) 

TENT2 (L) = ISPll(I.J.K) 

IENTRY(L)=WP1 (I. J.K) 

SET4 ( M ) = WPl(I.J.K) 

804 SFT1 <M) = ISPll(I.J.K) 

M = 0 

UO 806 K = I.J? 

L = L + 1 
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non 


i 

■f 

I 

} 


M = M + 1 

WP2(I*J.K) = WP(2.K) *' INTPT2(I*J*K)-/ INTPT2(NT0T*J*K) 
ISP2KI.J.K) = IT21 ( I * J*K ) / WP2< I * J*K) 

SAVE2 ( I * J*K ) = ISP21 ( I * J* K) 

TENT2(L) = ISP21 ( I * J*K ) 

TENTRYKL ) SWP2 ( I ► J * K ) . . . 

SET51M) = WP2(I.J*K> 

806 SET2 (M) = ISP21tI.J*K) 

M = 0 

DO 808 K s 1# J3 
L = L + 1 

M = M + 1 . 

WP3 ( I * J*K ) = WP (3*K) * INTPT3(I*J*K> / INTPT3(NT0T* J*K) 
ISP3KI.J.K) = IT31 ( I * J*K ) / . WP3( I * J*K ) 

SAVE3 ( I * J*K ) = ISP31 ( I * J»K ) 

TENT2XL) s ISP31(I*J*K) 

TENTRY ( L ) =WP3 ( I . J * K ) 

..SET6(M) = WP3 ( I * J * K ) . . 

808 SET3 <M> = ISP31(I*J*K> 

CALL STATS (SETl* Jl* MISll(I.J). SDIS1K I * J) ). 

CALL STATS (SET2* J2. MIS21(I*J>* SDIS21(I*J)) 

-CALL STATS (SET3. J3. MIS31(I»J)» SDIS311I.J)) 

CALL STATS (SET4* Jl * MWP1(I*J>* SDWP1(I*J)> 

. . CALL STATS. (SET5* J2* MWP2(I*J>* SDWP2(I*J)> 

CALL STATS (SET6* J3* MWP3(I*J>* SDWP3(I.J)> 

CALL STATS (TENT2* L> MIST1(I*J>* SDIST1(I*J)> 

CALL STATS (TENTRY *L* MWPT(I.J)* SDWPT(I*J>> 

882 CONTINUE 

PRINT PROPELLANT WT. AND SPECIFIC IMPULSE FOR PI 

DO 812 I = l'NTIMES 
WRITE (6*2037) I 
WRITE <6*2047) 

WRITE (6*2041) II* (MOTNO ( II *K ) * WPliI*l*K)* WP1(I*2*K>* 

L WP1 ( I*3*K) * BLANK* K = 1»J1) 

WRITE (6*2064) (MWP1(I*J)* J=l*4>* (SDWP1(I*J)* J=l*4) 

IF (NT ,LE. 1) 60 TO 810 

WRITE (6*2041) 12* (MOTNO ( 12 *K > * WP2(I*1*K>* WP2(I*2*K)* 

1 WP2 (I*3*K)* BLANK* K = 1*J2) 

WRITE (6*2064) (MWP2(I*J>* J=l*4)* (SDWP2(I*J)* J=l*4) 

IF (NT .LE. 2) 60 TO 810 

WRITE (6*2041) 13* (MOTNO( I3*K) * WP3(I*1*K)* WP3(I*2*K>* 

1 WP3 (I*3*K># BLANK* K = 1*J3) 

WRITE (6*2064) (MWP3(I*J)* J=l*4)* <SDWP3(I*J>* J=l*4) 

810 WRITE (6*2065) (MWPT(I*J)* J=l*4)* (SDWPT(I*J)* J=l*4) 

DO 811 J = 1 *NT 

CIMIN(J) = MWPT ( I * J) - CT1*SDWPT ( I * J) 

CIMAX(J) = MWPT ( I * J) + CT1*SDWPT ( I * J) 

C2MIN(J) = MWPT ( I * J) - CT2*SDWPT ( I * J) 

811 C2MAX( J) = MWPT ( I * J) + CT2*SDWPT ( I * J) 

WRITE (6*2062) C1MIN* ClMAX* C2MIN* C2MAX 

812 CONTINUE 

DO 1775 I = 1 *NT1MES 
WRITE (6*2034) I 
WRITE (6*2014) EXPA(l) 

WRITE (6*2047) 

WRITE (6*2041)11* (MOTNO(IlrK) * ISP11 ( I * 1 *K) * ISP11 ( I *2*K) * 

1 ISP11 (I *3*K) * BLANK. K=1*J1) 

WRITE (6*2064) <MIS11<I*J>* J=l*4)* (SDIS11 ( I * J) * J=l*4) 

IF (NT .LE. 1) 60 TO 1773 

WRITE (6*2041)12* (M0TN0(I2*K) r ISP21 ( I * 1 *K) r ISP21 (I *2*K) * 
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non 


i.J2> 


1 I5P21 (I . 3.K ) * BLANK . K = 

WRITE (6*2064) (MIS21(I*J>* J=l*4). (SDIS21 ( I * J) * J=l*4) 

IF (NT .LE. 2) GO TO 1773 
WRITE (6.2041)13* (M0TN0(I3*K) * ISP31 (I *1 *K ) * ISP31 ( I *?*K ) * 

1 ISP3KI.3.K) * BLANK. K = 1*J3) 

WRITE (6*2064) (MIS31(I.J). J=1.4). (SDIS31 ( I * J) . J=l»4) 

1773 WRITE (6.2065) (MIST1(I*J). J=1.4)» (S0IST1 ( I * J) . J=1.4) 

DO 1774 J = l.NT 

CIMIN(J) = MISTI ( I * J) - CT1*SDIST1 ( I * J) 

CIMAX(J) = MISTI ( I . J) + CT1*SDIST1 ( I * J) 

C2MIN ( J) = MISTI ( I * J) - CT2+SDIST1 ( I * J) 

1774 C2MAX ( J) = MISTI ( I * J) + CT2*SDIST1 ( I * J) 

WRITE (6*2062) C1MIN* C1MAX* C2MIN* C2MAX 

1775 CONTINUE 

IF (NP .LE. 1) GO TO 835 

COMPUTE THRUST AND IMPULSE FOR P2 ANP THEIR STATISTICS 

DO 1785 N = 2* NP 
DO 767 I □ l.NTIMES 
DO 767 J = 1*NT 
L = 0 

M = 0 

DO 712 K = 1.J1 
L = L + 1 

M = M + 1 

SAVE = FTll(I.J.K) 

FT12(I* J.K) = RATIOl (I,J.K)*FVAC1 ( I * 1 . K ) -FXPA <N ) *AETOT (1 *K > *CPHI 
IT12 ( I * J.K ) = SAVEl(I.J.K) *WP1 (I.J»K)*FT12(I*J*K) / SAVE 
TENTRY (L ) = FTIP(I.J.K) 

TENT2 (L ) = IT12(I.J.K> 

SET1 (M) = FT12 ( I * J* K ) 

712 SFT4 (M) = ITlPd.J.K) 

M = 0 

DO 716 K = 1*J? 

L = L + 1 
M = M + 1 

SAVE = FT21 (I. J.K) 

FT22 ( I * J.K ) = RAT 102 (I. J*K)*FVAC1 (I .2 . K ) -EXPA < N ) *AET0T ( 2 »K ) *CPHI 
IT22 ( I * J.K ) = SAVF,2(I.J.K) *WP2 ( I . J. K ) *FT22 ( I . J. K ) / SAVE 
TENTRY (L>= FT22(I*J.K) 

TENT2 (L ) = IT22 ( I . J.K) 

SET2 (M) = FT22 (I . J. K ) 

715 SET5 (M) = IT22 ( I . J.K) 

M = 0 

DO 718 K = 1.J3 
L = L + 1 
M = M + 1 

SAVE = FT3KI.J.K) 

FT32 (I . J.K ) = RATI03(I.J.K)*FVACld»3.K)-EXPA(N)*AET0T(3.K)*CPHl 
IT32 ( I . J.K ) = SAVE3(I.J.K) *WP3 ( I . J. K ) *< : 'T32 (I . J.K ) / SAVE 
TENTRY (L ) = FT3?(I.J.K) 

TENT2 ( L ) = IT32 ( I . J.K ) 

SET3 (M) = FT32C.J.K) 

718 SET6 (M) = IT32 ( I * J. K ) 

CALL STATS (SET1. Jl. MFT12(I.J). SDFT12(I.J)) 

CALL STATS (SET4 . Jl. MIT12(I.J). SPIT12(I.J)) 

CALL STATS (SET2. J2 . MFT22(I.J). snFT22(I.J)> 

CALL STATS (SETS. J2. MIT22(I.J). S01T22(I.J)) 

CALL STATS ( SET3 » J3. MFT32(I.J). SDFT32(I.J)) 

CALL STATS (SET6. J3. MIT32(I.J). SDIT32d.J)) 

CALL STATS (TENTRY. L. MFTT2(I»J). SPFTT2 ( I . J) ) 
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on non 


1 


CALL STATS (TENT2 ,Lr MITT2(I*J)* SDITT2(I*J)) 

767 CONTINUE 

PRINT THRUST* IMPULSE AND STATS FOR P2 

DO 1716 I = lrNTIMES 
WRITE ( 6 *2032 ) I 
WRITE (6*2014) EXPA(N) 

WRITE (6*2047) 

WRITE (6*2041) II* (MOTNO ( 1 1 * K ) * FT12 ( I * 1 * K ) # FT12(I*2*K)r 
1 FT 12(I*3*K) * RLANK* K=1*J1) 

WRITE (6*2064) (MFT12(I*J)* J = 1*4)* (SUFT12 ( I * J) * J = 1*4) 
IF (NT .LE. 1) GO TO 1716 

WRITE (6*2041) 12* ( MOTNO ( 12 * K ) * FT22(I*1*K)* FT22(I*?*K)* 

1 FT22 ( I * 3 * K ) * BLANK* K=1*J2) 

WRITE (6*2064) (MFT22(I*J)* J = 1*4)* (SuFT22U*J>* J = 1*4) 
IF (NT .LE. 2) GO TO 1716 

WRITE (6*2041) 13* ( MOTNO ( 13 * K ) * FT32(I*1*K>* FT32(I*?*K)* 

1 FT32 ( I * 3 *K ) * RLANK* K=1*J3> 

WRITE (6*2064) (MFT32(I*J)* J = 1*4)* ( SuFT32 ( I * J ) * J = 1*4) 

1716 WRITE (6*2065) (MFTT2(I*J)* J = 1*4)* ( SCFTT2 ( I * J ) * J = 1*4) 
DO 1717 J = 1 * NT 

CIMIN(J) = MFTTJ? ( I * J) - CT1*SDFTT2 1 1 * J) 

CIMAX(J) = MFTT 2 ( I * J ) + CT1*S0FTT2 ( I * J) 

C2MIN ( J ) = MF TT2 ( I * J ) - CT2*SDFTT 2 ( I * J ) 

1717 C2MAX ( J ) - MFTT2 ( I * J ) + CT2*SDFTT? ( I * J ) 

WRITE (6*2062) C1MIN* C1MAX* C2MIN* C2MAX 

1718 CONTINUE 

DO 1720 I = 1 * NTIMES 
WRITE (6*2033) I 
WRITE (6*2014) FXPA(N) 

WRITE (6*2047) 

WRITE (6*2041) II* (MOTNOdl *K) * IT12(I*1*K)* IT12tT,2,K)* 

1 IT12 ( I *3*K ) * PLANK* K=1*J1) 

WRITE (6*2064) (MIT12(J*J)* J = 1*4)* (SuIT12 ( I * J) * J = 1*4) 

IF (NT .LE. 1) GO TO 1715 

WRITE (6*2041) 12* ( MOTNO ( I 2 *K ) * IT?2(T*1*K>* IT22(I*2*K>* 

1 1T2? ( I * 3 * K ) * RLANK* K = 1*J2) 

WRITE (6*2064) (MIT22(I*J)* J = 1*4)* (SuIT22 ( I * J) * J = 1*4) 

IF (NT .LE. 2) GO TO 1715 

WRITE (6*2041) T3* ( MOTNO ( T 3 * K ) * IT32(I*1*K)* IT321I*2*K>* 

1 IT32(I*3*K)* PLANK* K=1*J3> 

WRlTe (6*2064) (MIT32(I*J)r J = 1*4)* ( SDI T32 ( I * J ) * J = 1*4) 

1715 WRITE (6*2065) (MITT2(I*J)* J = 1*4)* ( SDI TT2 ( I * J ) * J = 1*4) 

DO 171 Q J = 1 * NT 

C1M1N(J) = MITT2(I*J) - CT1 *SOITT? ( I * J ) 

CIMAX(J) = MITT 2(1 * J ) + CT1*SP ITT2 ( I * J ) 

C2MIN ( J ) = MITT2(I*J) - CT2*SPITT2 ( I * J) 

1719 C2MAX ( J ) = MITT 2 ( I * J ) + CT2*SDITT2 ( I * J ) 

WRITE X 6 * 2062 ) C1MIN* C1MAX* C2MZM* C2MAX 

1720 CONTIMUF 

COMPUTE SPECIFIC IMPULSE AND It’S STATISTICS FOR P2 

DO 897 I = 1 * NTIMES 
DO 897 J = 1 * NT 
L = 0 
M = 0 

DO 1810 K = 1*J1 
L = L + 1 
M = M + 1 

ISP12 ( I * J*K ) = IT12(I*J*K) / WP l ( I * J * K ) 
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n on 


SAVE1 ( I # J # K ) = ISP1 ? (I # J#K ) 
TENT2CL) = ISP32(I#J#K> 

1810 SET1 (M) = ISP J 2 d # vJ# K ) 


813 


814 


897 


M = 0 


DO 813 K = 1#J? 

L = L -f 1 
M = M + 1 

ISP22 (I # J# K ) = I T22 { I » Jt K ) / WP?d#J*K> 
SAVE2(I*J*K) = ISP2?d*J#K) 

TENT2 (L> = ISP?2(I*J#K) 

SET2 (M) = ISP?? (I # J # K ) 


M = 0 

DO 814 K = 1#J3 
L = L + 1 
M i M + 1 

1SP32 (I # J# K ) = IT32 ( I 
SAVE3d#J#K> = ISP32C 
TENT2(L) = ISP32 (I # J# ! 
SET3 (M) = ISP32d*Jr! 
CALL STATS ( SET 1 # Jl# 
CALL STATS (SET? r J2# 
CALL STATS (SET3# J3* 
CALL STATS (TENT2 # L# 
CONTINUE 


*J#K) / WP3 ( I # J#K ) 
IpJpK) 

<) 

<) 

MI SI 2 ( I # J) p SD1S12 ( I 
MIS22 ( I # J ) # SHIS??(I 
MIS32 (I # J ) t SD1S32 ( I 
MIST2 ( I # U ) p sniST2(I 


# J) 

p J ) 

# J> 
pJ) 


) 

) 

) 

) 


PRINT SPECIFIC IMPULSE FOR P2 


DO 1780 I = ltNTIMES 
WRITE <6r?034) I 
WRITE ( 6 r 2014 ) FXPA(N) 

WRITE ( 6 p 2047 > 

WRITE (6*2041)11# (MOTNO (II # K ) # 
L 

WRITE (6*2 064) (MIS12(I#J)# J= 
IF (NT ,LE. 1) GO TO 1778 
WRITE (6*2041)12# ( MOTNO (12 # K ) # 
L 

WRITE ( 6 p 2064 ) (MIS22(I#J)* J= 
IF (NT .LE. 2) GO TO 1778 
WRITE <6*2041)13# <MOTNO<I3#K)# 


ISP12(I#1#K> * ISPl?(I*?rK) 
ISP12 (I # 3 # K ) 9 PLANK # K = 1 
1*4)# (SDIS12 ( I # J ) 9 J=l*4) 

ISP22 (I r 1 # K ) * ISP?? (I # ? # K ) 
ISP2?(T*3#K) # BLANK, K = 1 
1 #4 ) # (SHIb?2( I # J ) 9 J=l#4) 


WRITE ( 6 #2064 ) <MIS32(I#J># 

1778 WRITE ( 6 # 2065 ) (MIST2(I#J)# 

DO 1779 J = 1 #NT 
CIMIN(J) = MIST2 ( I # J ) 

CIMAX(J) = MIST2 ( I # J) 

C2MIN ( J ) = MIST? ( I # J) 

1779 C2MAX ( J ) = MIST2<I*U) 

WRITE (6*2062) C1MIN* C1MAX* C2M1N# C?MAX 


ISP3?(I#1#K) # ISP32 ( I # ? # K ) 
ISP3? ( J # 3 # K ) 9 PLANK# K = 1 
J=l#4>* (SniS32(I#J> # J=l#4) 
J=l#4)# (SDIST2 ( I # J) # J=l*4) 


CTl*SDIST?d#J) 
CTl*SniST2(I#J) 
CT2*SDIST? ( I # J ) 
CT2*SDIST?(I#J) 


1780 

1785 

835 


CONTINUE 

IF (NP .LF. 

CONTINUE 

CONTINUE 

RETURN 

END 


2) GO TO 836 


#J1) 

* 

# J?) 

# 

# J3) 



SIBFTC M0T3 

SUBROUTINE M0T3 
DIMENSION XX (4860) 

COMMON G* CPHI * NT * EXPA* COUNT* WP* WEB* ATTOT * AETOT* 

1 P* MOTNO* FVAC1 * NTIMES* XX* NTOT * NBURNt NFIRST * NP 

COMMON CPI * C P2* CT1* CT2* VEC1* VEC2* VEC3 

WP1 (9*4*15)* WP2 (9*4*15) * 

CST1 (9*4* 15 ) * CST2 (9*4*15)* 

INTPT1 (9*4*15) * INTPT2 (9 * 4* 15) * 

ATTOT (4* 15) * AETOT (4*15) * 

VEC1 (4 * 15 ) * VEC2 (4*15) * VEC3 
R2 (4*15)* R3 (4*15)* MOTNO 

( 4 ) * 


DIMENSION 

DIMENSION 

1 FVAC1 (9*4*15>* 

2 P (9*4*15) * 

3 WEB (4 * 15 ) 9 

4 R1 (4*15>* 

DIMENSION 

1 C1TMN ( 4 ) * 

2 SET3 (15) » 

3 MRT ( 4 ) * 
DIMENSION 

1 

EQUIVALENCE 


WP3 (9*4*15) 
CST3 (9*4*15)* 
INTPT3(9*4* 15) * 

WP (4*15)* 

(4*15)* 
(4*15) 

C1MIN (4 ) * C1MAX (4 ) * C2MINI 4)* C2MAX( 4)* 


EXPA 

C1TMX ( 4)* C2TMN (4) * C2TMX(4)* SET1 (15)* SET2 ( 15 )* 
TENTRY (60 ) * TENT2 (60 ) *MR1 ( 4)* MR2 (4)* MR3 < 4)* 
SDR1 ( 4)* SDR2 ( 4)* SDR3( 4)* SDRT ( 4)* COUNT (4) 
MCST1 (9*4)* MCST2 (9 *4) • MCST3<9*4)* MCSTT(9*4)* - 
SDCST1 (9*4) *SDCST2 (9*4) * SDCST3(9*4) *SDCSTT(9*4) 


1 

2 

3 

4 

5 

6 


WPl ( 1 ) * CST1 ( 1 ) ) * 
WP2 ( 1 ) « CST2 ( 1 ) ) * 
WP3(l)r CST3 ( 1 ) ) * 
INTPT1 ( 1) ) • 

INTPT2 ( 1 ) ) * 

INTPT3 ( 1 ) ) 


► MCST2 * MCST3 
INTPT2* INTPT3 


* MCSTT * MR1 


* MR2 


MR3 


(XX( l)i 
(XX( 541)* 

( XX ( 1081 ) * 

( XX ( 3241 ) * 

(XX (3781) * 

(XX (4321) * 

REAL 

1 MRT * MCST1 . 

2 MOTNO* INTPT1* 

INTEGER COUNT 

OUTPUT FORMATS 

( 1H1 * 30X * 36H TRANSFORMED CHARACTERISTIC VELOCITY / 

9H TIME NO. 12) 

36H TRANSFORMED CHARACTERISTIC VELOCITY* 14X* 
33H TRANSFORMED AVERAGE BURNING RATE / 

9H TIME NO. 12* 70X* 18H BASED ON WEB TIME ) 
(11H0TEMP GROUP 12 // (1X*A6* 19X* 4F10.2* 10X* 4F10.4)) 
( 11H0TEMP GROUP 12 // (1X*A6* 

( 29X* 36H TEMPI TEMP2 

10H MOTOR NO.) 

( 29X* 36H TEMPI TEMP2 

36H TEMPI TEMP2 

10H MOTOR NO.) 

(34HOCONFIDENCE ON NORMAL DISTRIBUTION 


( 1H1 * 30X* 


2035 FORMAT 
1 

2036 FORMAT 
1 
2 

2040 FORMAT 

2041 FORMAT 

2047 FORMAT 
1 

2048 FORMAT 
1 
2 

2062 FORMAT 

1 26H (ONE 

2 26H SIDED 

3 26H (TWO 

4 26H SIDED 


19X * 4F10.2)) 

TEMP3 TEMP4 / 


TEMP3 

TEMP3 


TFMP4 

TEMP4 


* 17X * 


MIN 

MAX 

MIN 

MAX 


2064 FORMAT (26H0MEAN 

1 26H STANDARD DEV. 

2065 FORMAT (26H0T0TAL MEAN 

1 26H TOTAL STANDARD DEV. 

2066 FORMAT (26H0MEAN 

1 26 H STANDARD DEV. 

2067 FORMAT (26HOTOTAL MEAN 

1 26H TOTAL STANDARD DEV. 


4F10.2/ 

4F10.2 // 

4F10.2 / 

4F10.2 ) 

4F10.2 / 

4F10.2) 

4F10.2 / 

4F10.2 ) 

4F10.2* 10X* 4F10.4 / 

4F10.2* 10X* 4F10 .4 I 

4F10.2* IPX* 4F10(M J 
4F10.2* 10X* 4F10.4 ) 


2068 FORMAT ( 34HOCONFIDENCE ON NORMAL DISTRIBUTION / 


1 26H (ONE MIN 

2 26H SIDED) MAX 

3 26H (TWO MIN 

4 26H SIDED) MAX 


4F10.2* 10X*4F10.4/ 
4F10.2* 10X*4F10.4 
4F10.2*10X*4F10.4/ 
4F10.2*10X*4F10.4 ) 
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non non 


u = i 
12 = 2 
13 = 3 

J1 = COUNT (1) 

J2 = COUNT (2) 

J3 = COUNT (3) 
SDR114) = 0. 
MCSTT (5*4) = 0. 
BLANK = 0. 


COMPUTE R TRANSFORMED AND ITS STATISTICS 


DO 1055 J = If NT 
L = 0 
M = 0 

DO 1010 K = 1»J1 
L = L + 1 
M = M + 1 

Rl(JrK) = WEB(lfK) / VECl(JfK) 

TENT2(L) = RKJrK) 

1010 SET1 (M) = RKJ.K) 

M = 0 

DO 1015 K = ltJ2 
L = L + 1 

M = M + 1 

R2(J»K) = WEB ( 2 » K ) / VEC2(J»K) 

TENT2<L) = R2(J,K) 

1015 SET2 (M) = R2(JrK) 

M = 0 

DO 1020 K = 1 * J3 
L = L + 1 

M = M + 1 

R3(JrK> = WEB(3»K) / VEC3(JrK) 

TENT2(L) = R3 ( J»K ) 

1020 SET3 (M) = R3(JrK) 

CALL STATS (SETl* J1 1 MR1<J>» SDRUJ) ) 
CALL STATS (SET2* J2» MR2(J)» SDRtCJ)) 
CALL STATS (SET3r J3f MR3(J)» SDR31J) ) 
CALL STATS <T£NT2r L t MRT(J)» SDRT(J)) 


1055 CONTINUE 


COMPUTE C-STAR AND ITS STATISTICS 

DO 1000 I = N8URN» NTOT 
DO 1000 J = If NT 
L = 0 
M = 0 

DO 950 K = 1»J1 
L = L + 1 
M = M + 1 

CST1 ( I » J»K) = INTPT1 (NTOT* Jr K) * G * ATTOT(lfK) / Wf*l»K) 
TENT2(L) = CST1 < I » Jr K > 

950 SET1 (M) = CSTlU*JrK> 

M = 0 

DO 955 K = 1 t J2 
L = L + 1 
M = M + 1 

CST2 ( I » J>K ) = INTPT2 (NTOT » J »K ) * G* ATT0T(2rK) / WP (?»K) 
TENT2(L) = CST2(IfJfK) 

955 SET2 (M) = CST2(IrJ»K) 

M = 0 
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DO 960 K = 1.J3 
L = L + 1 
M = M + 1 

CST3 ( I . J.K ) = INTPT3CNTOT .J.K) * G * ATTOT (3.K) / WP(3.K) 

TENT2 (L ) = CST3 ( I . J .K ) 

960 SET3 (M) = CST3(I.J.K) 

CALL STATS (SET1. J1 . MCSTl(I.J). SDCSTl(I.J)) 

CALL STATS (SET2. J2. MCST2(I.J>. SDlST2(I.J)) 

CALL STATS (SET3* J3. MCST3(I.J). SDCST3(I*J)) 

CALL STATS (TENT2. L. MCSTT(I.J). SDCSTT(I.J)) 

1000 CONTINUE 

DO 1795 I = NBURN. NTOT 
DO 1785 J = l.NT 

C1MIN(J> = MCSTT(I.J) - CT1*SDCSTT(I.J) 

CIMAX(J) = MCSTT ( I » J) + CT1*SDCSTT ( I . J) 

C2MINIJ) = MCSTT ( I . J) - CT2*SDCSTT ( I . J ) 

1785 C2MAX ( J ) = MCSTT(I.J) + CT2*SDCSTT l I . J) 

IF (I .EQ. NBURN) GO TO 1786 

WRITE ( 6 . 2035 ) I 
WRITE (6.2047) 

GO TO 1788 

1786 WRITE (6*2036) I 
WRITE (6*2048) 

1788 IF (I .EQ. NBURN) GO TO 1791 

WRITE (6.2041 ) II* (MOTNO( I1*K) . CSTl(I.l.K). CSTKI.2.K)* 

1 CST 1 ( I * 3 * K ) * BLANK. K = 1,J1) 

WRITE (6*2064) (MCSTl(I.J). J=1.4). (SDCST1 ( I * J) . J=1.4) 

WRITE (6.2041)12. ( MOTNO ( 12 *K ) . CST2U.1.K). CST2(I*2*K). 

1 CST2 ( I * 3. K ) * BLANK. K = 1.J2) 

WRITE (6.2064) (MCST2(I.J). J=1.4). (SDCST2( I . J) . J=l»4) 

WRITE (6.2041)13. (M0TN0(I3.K) . CST3(i.l.X). CST3U.2.K). 

1 CST3(I.3.K). BLANK. K = 1.J3) 

WRITE (6.2064) ( MCST3 ( I . J ) . J=1.4). (SDCST3 ( I . J) . J=1.4) 

WRITE (6.2065) (MCSTT(I.J). J=1.4). (SDCSTT ( I . J) . J=1.4> 

WRITE (6.2062) C1MIN. C1MAX. C2MIM. C2MAX 
GO TO 1795 

1791 WRITE (6.2040)11 . (MOTNO(Il.K) . CSTl(l.l.K). CSTKI.2.K). CST 1(1.3. 
1 K). BLANK. R1U.K). R1(2.K). R1(3.K). BLANK. K = 1.J1) 

WRITE (6.2066) (MCSTl(I.J). J=1.4). MR1 . (SDCSTl ( I . J) . J=1,4),S0R1 
WRITE (6.2040)12. (M0TM0(I2.K) . CST2U.1.K). CST?(I.2.K). CST2(I.3. 
1 K). BLANK. R2 ( 1 . K ) . R2(2*K). R2(3.K). BLANK. K = 1.J2) 

WRITE (6.2066) (MCST2(I.J). J=l, 4 ), MR2 . ( SDCST2 ( I . J ) . J=1.4),sr)R2 

WRITE (6.2040) 13. (MOTNOi 13. K) . CST3(I.1.K). CST3(T.2.K). CST3(I»3. 
1 K). BLANK. R3(1»K>. R3(2.K). R3(3.K). BLANK. K = 1.J3) 

WRITE (6.2066) (MCST3U.J). J=1.4>. MR3. ( SDCST3 ( I . J ) . J=1.4).SBR3 

WRITE (6.2067) (MCSTT(I.J). J=1.4). MRT. < SDCSTT ( I . J ) . J=1.4).SPKT 

DO 1794 J = 1 .NT 
CITMN(J) = MRT(J) - CT1*SDRT(J) 

CITMX(J) = MKT(J) + CT1 *SORT ( J ) 

C2TMN ( J ) = MRT(J) - CT2*STtRT(J) 

1794 C2TMX ( J ) = MRT(J) + CT2*STIRT(J) 

WRITE (6.2068) C1MIN. C1TMM. CIMmX. C1TMX. 

1 C2MIN. C2TMM. C2MwX. C2TMX 

1795 CONTINUE 
RETURN 
END 
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APPENDIX B 


LISTING OF THE SOLID PROPELLANT ROCKET MOTOR PERFORMANCE 
VERSUS TIME COMPUTER PROGRAM USING THE 


GROUP TRANSFORMATION METHOD 


uuu u u u 


*P FOR MAIN 

GENERAL SOLID-PROPELLANT ROCKET MOTOR PERFORMANCE VS. TIME 
DIMENSION XX (1700 ) 

DIMENSION XLNP (3.15.200). FSAVE (3.15.200) . PSAVE (3. 15.200 ) 
DIMENSION CC (200*4) * PCO(200)* FOC200J* Y(200)* 

1 WPRNT (200) * TREP(200)» PC (200)* T(200)* F(200) 

DIMENSION AE1 (3*15) * AE2(3*15>* AE3(3.15)» AE4<3.15). WEB(3*15)* 

1 XM0TN0(3*15) * TEMP (3* 15) * PAF (3*15) *TAIL(3*15) * AETOT (3.15) 

DIMENSION TGRP ( 4)* EXTEM ( 3). EXPA ( 3)* K0UNT( 4)* EPS( 4) 
DIMENSION BCDX ( 12) * 8CDK12). BCD2(12)* BC03(12). 0CO4(12) 
DIMENSION PCTW(IO). NW(10). BREAKW ( 10 ) * WINC(IO) 

DIMENSION PCTT ( 5). NT(5>* BREAKT ( 5). TINC( 5) 

DIMENSION NGO ( 10 ) * NSTOP(IO)* XNW(IO)* XNT ( 5). SETK10) 

COMMON XX 

COMMON XLNP. FSAVE. PSAVE 
COMMON BCDX. BCD1 * BCD2* BCD3* BCD4 
EQUIVALENCE 


1 

( XX ( 1). 

NSV2 ) # 



(XX( 3). 

ICON) * 

2 

( XX ( 4). 

TAV) # 

(XX( 5). 

Nl) . 

(XX( 6). 

MPRIME) . 

3 

( XX ( 7). 

NPTS ) # 

<XX( 8). 

Tfl) . 

(XX( 9). 

TT) . 

4 

(XX< 10). 

WPCT ) # 

(XX( 11). 

TPCT), 



5 

<XX< 13). 

CT2 ) » 

(XX( 14). 

TREP) . 

(XX( 214). 

T) . 

6 

( XX < 414). 

CO# 

( XX ( 1214 ) . 

NSV1 ) . 

(XX ( 1215 ) * 

NXNT) . 

7 

(XX ( 1216) * 

NA) # 

( XX ( 1217 ) . 

KOUNT) . 

(XX ( 1221 ) . 

TEMP) . 

8 

( XX ( 1266 ) * 

IPKNT) # 

( XX ( 1267) . 

N2 ) . 

(XX ( 1268) . 

TO) . 

9 

(XX ( 1269) * 

NP) # 

(XX (1270) . 

EXPA) . 

(XX ( 1273) . 

AETOT) 

EQUIVALENCE 






1 

( XX ( 1318 ) * 

PAF) # 

( XX ( 1363) . 

CPHI ) . 

(XX (1364) . 

TO) . 

2 

( XX ( 1365) . 

WPRNT) # 

( XX (1565 ) . 

XMOTNO) . 

(XX ( 1601 ) . 

NXNW) 

EQUIVALENCE 






1 

( XX ( 1802 ) * 

BREAKW) # 

( XX ( 1612 ) . 

BREAKT) . 

( XX ( 1617) . 

WINC) , 

2 

(XX ( 1627 ) . 

TINC) # 

(XX ( 1632 ) . 

NGO) . 

OCX (1642) . 

NSTOP) . 

3 

( XX ( 1652 ) * 

PCTW ) # 

(XX (1862) . 

PCTT) 




INPUT FORMATS 






81 

FORMAT 

(112. 4F12.4) 


82 

FORMAT 

( A6 . 6F12.4. IX. 

ID 

83 

FORMAT 

(216. 3F12.4.I6) 


84 

FORMAT 

( 5 (F10 .4 . 16)) 


85 

FORMAT 

(4E12.5) 


100 

FORMAT 

(316) 


101 

FORMAT 

(2 (E16.8 . E12.5. 

E12.5) ) 

102 

FORMAT 

( 12A6 ) 



OUTPUT FORMATS 


109 FORMAT (38H1THE FOLLOWING OUTPUT IS FOR MOTOR NO. A6> 


110 

FORMAT (54H0 

PCT WEB TIME 

TIME 

PC 

THRUST/) 

111 

FORMAT (F12 

• 4 # 3F14 «4) 




112 

FORMAT (54H0 

PCT TAILOFF TIME 

TIME 

PC 

THRUST/) 


1 READ (5. 83) NXNW. NXNT * TD. PHI* CT2.NSW 
IF (NSW .GT. 0) READ (5.84) RSIA 


WRITE (6 . 

83) 

NXNW. NXNT. TD. 


PHI 

t 

CT2 

READ 

<5# 

84) 

{PCTW ( I ) # 

NW ( I ) 

# 

I 

z 

1 »NXNW ) 

WRITE (6 . 

84) 

( PCTW ( I ) # 

NW(I) 

# 

I 

z 

1 # NXNW ) 

READ 

(5# 

84) 

(PCTT ( I ) # 

NT ( I ) 

# 

I 

z 

1 # NXNT ) 

WRITE ( 6 # 

84) 

(PCTT ( I ) # 

NT ( I ) 

# 

I 

z 

1 # NXNT ) 

READ 

(5# 

81) 

NA r (TGRP 

(I)# 

I 

z 

1 #NA ) 
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WRITE (6# 81) 
READ <5* 81) 
WRITE (6* 81) 
READ (5*102 ) 
WRITE (6* 102 ) 
READ (5*102) 
WRITE (6* 102) 
READ (5*102) 
WRITE (6*102 ) 
READ (5*102) 
WRITE (6*102) 
READ (5*102) 
WRITE (6* 102) 


NA $ (TGRP 
NP* (EXPA 
NP* (EXPA 
(BCDX ( I ) * I 
(RCDX(I) *1 
(BCD1 ( I ) * I 
(BCDl(I) *1 
(BCD2( I ) * I 
( BCD2 ( I ) * I 
(BCD3 ( I ) * I 
<BCD3(I) *1 
( BCD4 ( I ) * I 
(RCD4 ( I ) * I 


(I)* I 
(I)* I 
(I)* I 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 
= 1 * 12 ) 


1 *NA ) 
1 *NP ) 
1 *NP) 


KOUNT(l) = 0 
KOUNT ( 2 ) = 0 
KOUNT ( 3 ) = 0 
NGO ( 1 ) =1 

DO 4 I = 1 * NXNW 

XNW(I) = NW(I) 

WINC(I) = 1. / XNW ( I ) 

NGO ( I +1 ) = NGO(l) + NW(I) + 1/1 

N5T0P(I) = NGO(I+l) r 1 
CONTINUE 

DO 5 I = InMXNT 
XNT(I) = NT ( I ) 

TINC(I) = 1. / XNT(I) 

I SUB = (NXNW +1+1) 

NGO ( I SUM ) = NGO (NXNW + I) + NT ( I ) 
i\lSTOP(fJXrjW + I) = NGO(ISUH) - 1 
CONTINUF 


IPASS = 0 

CPU! = COS (PHI) 

IMSUM = NXNW + MXNT 
NS VO = NSTOP(MXNW) 

N1 = NS VO - 1 
NSVl = NSVO + 1 
NSV2 = NSTOP(NSUM) 

DO 7 I = 1 * NXNW 

SET 1 ( I ) = PCTW ( I ) * ."1 

DO 8 I = 1 * NXNW 

PCTW ( T + 1 ) = SET 1(1) 

PCTWll) = 0. 

DO 9 I = 1 * MXNT 
SETHI) = PCTT(I) * .01 
DO 10 I = 1 *NXNT 
PCTT ( I +1 ) = SET) (I) 
PCTTU ) = n. 


GROUP MOTORS INTO TEMPERA 1 URF GkOUPS 
TGRP(MA+i> = TGRP(NA) * 2. 

READ (5* 82) XMOT* FI RE TP * PAX* IP* IT* TOrPCP* NOMOT 
WRITE (6* 8?) XMOT* FIRFTP* PAX* TF * IT* |0*PCP* NPM01 
DO 19 J = 1* NA 

fc.PS(J) = (TGRP ( J+l ) - TGRP(J)) *.5 

IF ( (TGRP(J) + tPS(J)) - FIRFTP) 19* 19* 25 

CONTINUF 

KOUNT (J) = KOUNT (J) + 1 
K = KOUNT (J) 

READ (5* 85) AE1 ( J * K ) * AE2(J*K)r AE3(J*K)* AE4(J*K) 



WRITE (6* R5) AE1 ( JrK ) r AE2(JrK)r AE3 ( Jr K ) r AE4 1 JrK ) 

TEMP (JrK) = FIRETP 
XMOTNO ( JrK ) = XMOT 
PAF (JrK) = PAX 

WEB (JrK) = TB 

TAIL (JrK) = TT 

AETOT (JrK) = AFl(JrK) + AE? ( Jr K ) + AF3 ( JrK ) + AE4 ( JrK ) 

C 

READ ( 5 r 100 ) NPTSr ICONr IPRNT 
WRITE (6r 100 ) NPTSr ICONr IPRNT 

READ ( 5r 101 ) (TREP(I)r PC ( I ) r F(I)r l=lrNPTS) 

C 

DO 11 1=1 rNPTS 

11 PC(I) = PC (I) + PCP 
PRES1 =' 0.0 
PRES2 = 0.0 

IF ( NSW .GT. 0) GO TO 5005 
5001 EREAKW ( 1 ) = TO 

DO 26 I = lrNXNW 

26 HREAKW ( I +1 ) = TO + PCTW ( 1+1 ) *TB 
UREAKT(l) = TB + TO 
DO 2H I = lrNXNT 

2P BREAKT ( 1+1 ) = TR + PCTT ( 1+1 ) * ( TT-Tb ) + TO 
C 

IF ( PRES1 .NE. 0.0) GO TO 5002 
5005 CALL CURVE (TREPr PC r NPTS) 

IF (NSW .LE. 0) GO TO 5002 
IONEW = TRFP(l) - .00001 - TAV 
DO 5003 1=1 r 100000 

TOMEW = TONEW + .00001 

RRCS1 = ( (CC(lr4)*T0NEW + CC(lr3)> *TONEW + CC ( 1 r2).) *TOMEW+CC ( 1 r 1 ) 
PRES2 = CC ( 1 r 1 ) + CC(lr?)*TONEW + CC ( 1 r 3 ) *T0NEW»*2 
1 + CC ( 1 r 4 ) *T 0MEW**3 

IF (PKFS1 ,GF. PSIA .OR. PRES2 .GE. PSIA) GO TO 5004 
50U3 CONTINUE 
5004 TO = TONEW + TAV 

WRITE ( 6 r 1 0 02) TOrPSIArPRESlr PRES2 

1002 EORMnT ( 1H1 /// 3llTO= FJ 0 . 4 r 6X r 5HPSI A= Flo . 3r 6X r 6HPKES1= F?0.4r 
1 6Xr6HPKFS2= F?0 . 4//////// ) 

GO TO 5001 

6002 CALL SRCHl(PCO) 

DO 30 L = lr NSV2 
aRG = PC 0 ( L ) 

('SAVE (JrK rL ) = ARG 
30 XLNP(JrKrL) = ALOG (ARG) 

C 

CALL CURVE (IRFPr F r NPTS) 

CALL SRCH1 (FO) 

DO 35 L = lr MSV2 
35 FSAVE ( Jr K r L ) = FO(L) 

C 

WRITE ( 6 r 109) XMOTNO (JrK) 

WRITE ( 6 r 110) 

DO 56 IJ = lr MXNW 
A I = 0. 

NG = NGO(IJ) 

NSP= MSTOP(IJ) 

UELPw = PCTW ( I J+l ) - PCTW(IJ) 

DO 55 I = NG t NSP 
IF (IJ .NE. 1) GO TO 53 
IF (I .EIO. 1) GO TO 54 
53 A I = A I + 1 
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54 IF (IPASS .EQ. 1) GO TO 55 

WPRNT(I) =(PCTWCIJ) + WINC ( I J) * A I*UELPW ) * 100. 

55 WRITE (6* 111) WPRNT ( I ) * T(I)* PCO(I>* FO(I) 

56 CONTINUE 
WRITE <6* 112) 

UO 60 IJ = 1* NXNT 
A I = 0. 

NG = NGO (IJ + NXNW ) 

NSP = NSTOPCIJ + NXNW) 

DELPT = PCTTlIJ+l) - PCTT(IJ) 

UO 58 I = NGf N5P 
IF (IPASS .EQ. 1) GO TO 58 
AI = AI + 1. 

WPRNT ( I ) =(PCTT(IJ) + TINC ( I J) *AI*uFLPT) * 100. 
58 WRITE (6 * 111) WPRNT ( I ) * T(I>* PCO(I)* FO(I) 

60 CONTINUE 
IPASS = 1 

IF (NDMOT ) 65* 15* 65 
65 CALL LSTSO ( 1* NSVO* WEB * 1) 

CALL LSTSO (NSV1. NSV2* TAIL* 2) 

GO TO 1 
END 
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•IP FOR SRCH1 

SUBROUTINE SRCH1 (Y) 

DIMENSION XX (1700) 

DIMENSION T(200)# CC<200#4)# TREP1200)# Y(200) 

DIMENSION BREAKW(IO) # BREAKT (5) # WINC ( 10 ) # TIMC ( 5 ) # NGO ( 1 0 ) # NSTOP ( 10 ) 
COMMON XX 


EQUIVALENCE 

1 (XX( 4 ) # 

2 (XX( 14 ) # 

3 (XX ( 1268 ) # 

4 (XXQ601)# 

5 (XX< 11)# 
EQUIVALENCE 

1 (XX<1602)# 

2 ( XX ( 1627 ) # 


TAV) # 
TREP) t 
TO) # 
NXNW) r 
TPCT) # 

BREAKW ) # 
TINC) # 


(XX( 8) # 
(XX( 214 ) 9 
( XX ( 1215 ) f 
(XX< 7 ) 9 
< XX ( 3># 

( XX ( 1612 ) # 
( XX ( 1632 ) # 


Tb) # 
T) r 
NXNT) # 
NPTS) # 
ICON) 


(XX( 
(XX( 
( XX ( 
(XX< 


10 ) # 
414) * 
5 ) # 
9) # 


WPCT ) p 
CO 9 
Nl) * 
TT ) «! 


BREAKT) # 
NGO) 9 


( XX ( 1617 ) 9 
( XX ( 1642 ) # 


WINC) # 
NSTOP) 


NX = ICON - 1 
XI = -1. 

NSUM r NXNW + NXNT 
DO 45 IJ = 1* NSUM 
XI = 0. 

NG = NGO(IJ) 

NSP= NSTOP ( I J) 

DELW = BREAKW ( I J + l ) - BREAKW ( I J) 

IK = IJ - NXNW 

DELT = BREAKT ( IK+1 ) - BREAKT ( IK ) 

DO 40 I = NG 9 NSP 
IF (IJ .GT. NXNW) GO TO 15 
IF (IJ «NE • 1) GO TO 8 
IF ( I *EQ • 1) GO TO 9 

8 XI = XI + 1. 

9 T ( I ) = BREAKW < I J) + XI*WINC ( I J) *DELw 

GO TO 21 

15 XI = XI + 1. 

T ( I ) = BREAKT ( IK ) + X I *TINC ( IK ) *DELT 

21 DO 28 J = 1 r NPTS 

IF (TREP(J) - T ( I ) ) 28 9 28# 22 

22 JS = J - 3 

IF (JS) 23# 23# 25 

23 JS = 1 
GO TO 32 

25 IF (JS - NPTS + NI) 32# 32# 26 

26 JS = NPTS - NI 
GO TO 32 

28 CONTINUE 

32 TQ = T(I) - TAV 

Y ( I ) = ( (CC< JS#4)*TQ + CC ( JS # 3 ) ) *TQ + CC(JS#2))*T© + CC(JS#1) 
40 CONTINUE 
45 CONTINUE 
RETURN 
END 
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•IP FOR LSTSQ 

SUBROUTINE LSTSQ (NGO, NSTOP, TIME, IPASS) 

DIMENSION XX (1700) 

DIMENSION XLNP (3»15»200>» FSAVE (3, 15* 200 ) . FD1 ( 3, 15,200 ) » 

1PSAVE(3» 15,200), F02 (3,15,200), FD3(3, 15,200) , PD(3,15,200) 

DIMENSION C2MAX ( 200,3) , C2MIN(200,3) , XMFD(200»3>» 


1 

PIKD(3»200) , 

XMLNP (3# 

200) # 

XMOTNO ( 

3# 15) # 

AETOT ( 3,15), 

2 

PAF (3, 15), 

TEMP (3# 

15) # 

TIME (3 # 15)# 

ARRAY (46# 

4) # 

3 

XLNPD (200 ) , 

C1MIN 

(200 ) # 

C1MAX 

(200) # 

WPRNT (200) , 

4 

XMTD (200 ) , 

C3MIN 

(200 ) # 

C3MAX 

(200)# 

XMPTD(200) , 

5 

SET1 ( 45), 

SET2 

( 45)# 

XMTEMP 

( 3)# 

OENOM ( 

3)# 

6 

EXPA ( 3), 

SIGP 

( 3)# 

X 

( 4)# 

KOUNf ( 

4) # 

7 

IL( 4), 

LGO 

( 10)# 

PCTW 

( 10)# 

PCTT ( 

5) # 

8 

WINC ( 10), 

TINC 

( 5) # 

SET3(45) » TEMX (45) 


DIMENSION BCDX ( 12 ) , 

0COK12) 

, BCD2 ( 12) , BCD3 ( 12) , BCD4(12> 



COMMON XX 

COMMON XLNP, FSAVE, PSAVE 

COMMON BCDX, BCD1, 8CD2, BC03, BCD4 

EQUIVALENCE 


1 

( XX ( 5), 

Nl) r 

C XX C 11), 

TPCT ), 

(XX( 10)# 

WPCT) # 

2 

(XX ( 1216) , 

NT) # 

( XX ( 1217 ) , 

KOUNT ) , 

( XX ( 1221 ) # 

TEMP) # 

3 

( XX ( 1269 ) , 

NP) # 

< XX ( 1270 ) , 

EXPA) , 

( XX ( 1273) # 

AETOT) # 

4 

(XX ( 1318) , 

PAF) » 

(XX ( 1363 ) , 

CPHI ) , 

( XX ( 1364 ) # 

TD) # 

5 

( XX ( 1 365) , 

WPRNT) # 

(XX ( 1565 ) , 

XMOTNO) , 

<XX( 13)# 

CT2) # 

6 

( XX ( 1), 

NSV2 ) # 

( XX ( 1214) , 

NSV1 ) , 

l XX ( 160 1 ) # 

NXNW) # 

7 

( XX ( 1632 ) , 

LGO) r 

, (XX ( 1652 ) , 

PCTW) , 

( XX <1862 ) # 

PCTT) # 

8 

(XX ( 1617 ) , 

WINC) # 

( XX ( 1627) , 

TINC) 



EQUIVALENCE 

( XLNP ( 1 ) # 

FDl(l), FD2 ( 1 ) , FD3 ( 1 ) , PD(1)) 



C 

C OUTPUT FORMATS 

C 

200 FORMAT ( 1H1 , 13X, 56H TRANSFORMED TIMES TRANSFORMED CHAMB 

1ER PRESSURE ) 

201 FORMAT ( 7H PCT., 16X, 38H MEANS WITH TWO SIDtD TOLERANCE LIMITS/ 

1 36H WEB MEAN MIN. MAX. ,7X, 30H MEAN MIN. 

2 MAX. / 7H TIME ) 

202 FORMAT (1H1, 13X, 54H TRANSFORMED TIMES TRANSFORMED THRUS 

IT AT PA =FB.2) 

205 FORMAT (IX, F6.2, 3F10.4, 3F12.4) 

1201 FORMAT ( 5H PCT. / 8H TAILOFF / 5H TIME ) 

C 

NS VO = N1 + 1 
IF (IPASS - 2) 66, 72, 66 
66 NM = KOllNT(l) + KOUNT(2> + KOUNT(3) 

C 

C COMPUTF STATISTICS FOR TEMPtRATURE 

C 

DO 70 J = 1 »NT 
LL = 0 

JJ = KOUNT(J) 

DO 68 K = 1,JJ 
LL = LL + 1 

68 SETKLL) = TEMP(J»K) 

CALL STATS (SET1, LL, XMTEMP(J), AR6) 

DENOM(J) = TD - XMTEMP(J) 

70 CONTINUE 
C 

C LEAST SQUARE FIT LOG OF BURN TIME VS. TFMP • 

C 

72 LL = 0 

DO 210 J = 1 ,NT 
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non non nnoonoo 


JJ = KOUNT(J) 

DO 210 K = 1 * J J 
LL = LL + 1 
ARRAY (LL* 1 ) = 1. 

ARRAY (LL*2) = TEMP(J*K> 

ARRAY (LL*3) = TEMP(J*K)**2 
ARG = TIME ( J*K ) 

ARRAY (LL. 4) = ALOG(ARG) 

210 CONTINUE 

CALL 6LS1 (ARRAY* X. IL* NM* 3. ALPHA* 0.* 0.) 

XLNPD(l) = ( X ( 3 ) *TD + X(2))*T0 + X(l) 

COMPUTE LOG OF MEANS OF BURN TIMES AMU SIGP*S 

DO 230 J = 1*NT 

XMLNP ( J* 1 ) = ( X ( 3 ) *XMTEMP ( J ) + X (2 )) *XMTEMP ( J) + X(l) 

230 SIGP(J) = ( XMLNP ( J* 1 ) - XLNPD(l)) / uEMOM(J) 

COMPUTE TRANSFORMED TIMES AND IN THE SAME LOOP LFAST SQUARE 
FIT LOG OF PRESSURE VS. TEMP 

IJ = 0 

DO 240 L = NGO* NSTOP 
IQ = IJ + 1 + ( IPASS-1 ) *NXNW 
IF (L .Nt. LGO(IQ)) GO TO 231 
IJ = IJ + 1 
AL = 0. 

I NCR = 1 / IJ 
XINC r I NCR 

231 AL = AL + 1 
LL = 0 

DO 235 J = 1 *NT 
JJ = KOUNT(J) 

DO 235 K = 1*JJ 
LL = LL + l 

ARRAY (LL* 4 ) = XLNP(J*K*L) 

ARG = SIGP ( J) * (TO - TEMP( J*K ) ) 

TEN = TIME ( J* K ) / (EXP (ARG)) 

IF (IPASS - 1) 232* 232* 233 

232 TEMX(LL) = TEM 

DELW = PCTWUJ+1) - PCTW(IJ) 

SETKLL) = TEM * (PCTW(IJ) + ( AL-XINC ) *WINC ( I J) *OELW ) 

GO TO 235 

233 DELT = PCTT(IJ+1> - PCTT(IJ) 

TAM = TEM - TEMX(LL) 

SETKLL) = TEMX(LL) + TAM * (PCTTUJ) + AL*T INC ( I J ) *DELT ) 

235 CONTINUE 

CALL STATS (SET1* LL* XMTD(L) * SUTU) 

CIMIN(L) = XMTn(L) - CT2*SDT0 
CIMAX(L) = XMTP(L) + CT2*SDTQ 

CALL GLS1 (ARRAY* X* IL* NM* 3* ALPHA* 0.* 0.) 

XLNPO(L) = ( X ( 3 ) *TO + X(2))*TD + X(l) 

COMPUTE XMLNP FRQM CURVE FIT AND PIK*S 

DO 240 J = 1* NT 

XMLNP ( J *L ) = (X(3)*XMTFMP(J> + X (2 ) ) *XMTtMP ( J) + X(l) 

PIKO(J*L) = (XLNPD(L) - XMLNRT J*L ) ) / DENOM(J) 

240 CONTINUE 

COMPUTE AND PRINT OUT TRANSFORMED PRESSURES 
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non 


C 


90 


mo 


UO 100 L = NGO 9 NSTOP 
LL = 0 

DO 90 J : 1 r NT 
JJ = KOUNT(J) 


DO 90 K : 
LL = LL + 
ARG = PIKD 
PD(JfKfL) : 
SET3 (LL) = 
CONT1NUF 
CALL ST * TS 
C3M IN ( L ) = 
C3MAX ( L ) = 
CONTINUE 


f JJ 

J f L ) * (TO 
: EXP (XLNP 
PO( JfKfL) 

(SET3 

XMPTD(L) - 
XMPTD(L) + 


- TEMP(JfK) ) 
(JfKfL) + ARG) 


f NMf XMPTD (L ) f SPPTP) 
CT2*SDPTD 
CT2+S0PTD 


IF (IPASS - 2) 182 f 101 9 18? 

101 WRITE ( 6 f 200 ) 

WRITE (6r 201) 

DO 10B L = If NSV2 

IF (L - NSV1) 1 0 5 r 104# 105- 

104 WRITE ( 6 f 1201) 

105 WRITE (6# 205) WPRNT(L)# XMTD(L)# 
1 XMPTD (L ) f C 3M I N ( L ) f C3MAX(L) 

CALL GUIKMV <-1f 46 f RCDXf BCDlr 

CALL GUIKMV ( Of 67 f RCDXf RCDIf 

CALL GUIKMV ( Of 77 f RCDXf BCDIf 


C1MIN(L)f C1MAX(L)f 

— NSV? f C 1MAX f C3MAX) 
-NSV2f XMTDf XMPTD ) 
—NSV? f C1MIN f C3MIN) 


COMPUTE TRANSFORMED THRUSTS 


182 


184 

186 

190 

194 

196 

198 


UO 198 L = NGOf NSTOP 
LL = 0 

UO 190 J = 1 f NT 
JJ = KOUNT(J) 

LO 190 K = 1 f JJ 

LL = LL + 1 

F<AT 10 = PD(JfKfL) 

F VAC = FSAVE(JfKfL) + 
FD1(JfKfL)= RATIO+FVAC 
SET 1 ( LL ) = FD1 (JfKfL) 
IF (NP - 1 ) 190 f 190 f 
FD2(JfKfL)= RATIO*FVAC 
SET2CLL) = FD2 (JfKfL) 
IF (NP - ?) 190 f 1 90 f 
FD3(JfKfL)= RATIO*FVAC 
SET3 ( LL ) = FL)3 (JfKfL) 
CONTINUE 

CALL STATS (SETIf NMf 
C2MIN ( L f 1 ) = XMFU(LfI) 
C2MAX(Lf1) = XMFU(LfI) 
IF (NP - 1) 198 f 198 f 
CALL STATS (SET2f NMf 
C>MIN(Lf2) = XMFD ( L f 2 ) 
C2MAX ( L f 2 ) = XMFU(Lf2) 
IF (NP - 2) 198 f 198 f 
CALL STATS (SET3 
l2MIN(Lf3> = XMFD <Lf 3 ) 
C2MAX ( L f 3 ) = XMFU(Lf3) 
CONTINUE 


/ PSAVE(JfKfL) 

PAF ( JfK ) *AETOT ( JfK ) *CPH1 

- EXPA(1)*AETOT(JfK)*CPMI 

184 

- EXPA(2)+ALTOT(JfK)*CPHI 
186 

- EXPA(3)*AET0T( JfK)*CPHI 


XMFD(Lf1)f SuFu) 

- CT2*SDFU 
+ CT?*SDFU 

194 

XMFD ( L f 2 ) f SuFu) 

- CT2*SDFU 
+ CT?*SDFU 

196 

f NMf XMFD(Lf3)f SuFU) 

- CT?*SDFu 
+ CT?*SDFU 


IF (IPASS - 2) 209# 301 f 209 


301 M = 0 
204 M = M + 1 
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WRITE (6. 202) EXPA(M) 

WRITE (6* 201) 

DO 206 L = 1* NSV2 
IF (L - NSV1 ) 206* 207* 206 
207 WRITE (6* 1201) 

206 WRITE (6* 205) WPRNT(L) * XMTD(L) * CIMIN(L)* CIMAX(L)* 
1 XMFD (L *M) * C2MIN(L*M)* C2MAX(L*M) 



IF (M - NP) 

204 r 

2l8 r 

218 



218 

CALL 

QUIKMV 

( 

-If 

46 r 

BCDXr 

BCD2* -NSV2*C1MAX* 

C2MAX (1*1)) 


CALL 

QUIKMV 

( 

Or 

67 r 

RCDXr 

BCD2* -NSV2* XMTD* 

XMFD (1*1)) 


CALL 

QUIKMV 

( 

Or 

77 r 

RCDXr 

BCD2 * -NSV2 *C1MIN* 

C2MIN ( 3*1)) 


IF (NP - 1) 

209 r 

209 r 

220 



220 

CALL 

QUIKMV 

( 

-lr 

46 r 

RCDXr 

BCD3* -NSV2.C1MAX* 

C2MAX ( J *2) ) 


CALL 

QUIKMV 

( 

Or 

67 r 

RCDXr 

BCD3* -NSV2* XMTD* 

XMFD (1*2) ) 


CALL 

QUIKMV 

( 

Or 

77 t 

RCDXr 

BCD3* -NSV2 * C1MIN * 

C2MIN ( 1 *2) ) 


IF (NP - 2) 

209 r 

209 r 

222 



222 

CALL 

QUIKMV 

( 

-lr 

46 r 

RCDXr 

BCD4 * -NSV2*C1MAX* 

C2MAXO *3) ) 


CALL 

QUIKMV 

( 

Or 

67 r 

BCDXr 

BCD4* -NSV2* XMTD* 

XMFD (1*3) ) 


CALL 

QUIKMV 

( 

0, 

77 r 

RCDXr 

BCD4* -NSV2*C1MIN* 

C2MIN ( 1 *3) ) 


209 RETURN 
END 
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•IP FOR CURVF 

SlMROUTlNF 
C CURVE FIT 

DIMENSION 
DIMENSION 

CURVF (Tr 

XXU 700) 
ALi) ( 40 0 ) r 

Sr N) 

DAL ( 40D ) r 

ALC (400) r 

ALW (1000) » 

TC (200) r 

1 

T<200) r 

s( 2 .on) f 

A3(?a0) r 

SC( 200) r • 

DP ( pno ) f 

2 

COMMON XX 

DtRC200) r 

CC ( ?0n f 4 

) 



EQUIVALENCE 






1 (XX( 3)r ICON)* (XX( 4) r TA V ) r (XX( 414)r CC)r 

2 (XX(1266)r I PR NT ) 

WRITE ( 6 r 1 1 1 ) 

110 FORMAT (2 F15.5) 

111 FORMAT ( 1H3 r 12HINPIJT ARRAYS) 

WRITE ( 6 r 1 1 0 ) (T(I)rS(I)* I=lrN) 

***THE FLFiMENTS OF THE T AND S ARRAYS ARF ARRANGED TIM ASCFMniNG ORUFR 
r A NO ARE THEN NORMALIZED. 

CALL ACCEND ( T »S »N) 

CALL NORMLZ(TrTCrNrTAV) 

CALL NORMLZ(SrSCrNrSAV) 

UO 1 1=1 rN 
AL0(2*I )=TC ( I ) 
aLD (2*1-1) =SC { I ) 

UAL ( 2*1 )=SC ( I ) 

1 UAL(2*I-1)=TC(T) 

HERE BEGINS PIFCEWISE KITS WITH CUolCSr lMoRACING M POINTS Aj TIME 
K = 1 
M = ICON 

IF ( M-2* ( M/2 ) ) 23 9 24 f 23 

23 IS=0 
00 TO 25 

24 IS=1 

25 NH=M/2 
MK = M 
NM=N+1-M 

IF (IPRNT) 6 29 64r 02 

62 WRITE ( 6 r 63) TAVr M 

63 FORMAT ( 76H1 FOLLOWING ARF C OF FF I C I DOTS OF SLIGHTLY uTS^OMT INMO'IS CH 

6311. ICS IN THE QTY ( T IME-F 1? . 7 r 2611 ) r FOR TImF RANOFS SHO'-'h . /A6H0T IMt 
632KANGE COEFF. OF Q**0 COFFF. OF 0**1 COt-FF . OF Q**P C 

6330EFF. OF Q**3rI4r22H POINTS LINKtP PtF Sc.T//) 

♦ ♦♦COMPUTING THE COEFFICIENTS OF THE CUPICS AND STORING THFM IN rHh 
CC ARRAY 

64 UO 52 1 = 1 r MM 

CALL CF2F1 (OrALCrOrA3rALD(2*I-l ) r3rM) 

CALL CF2F2(0.n. ALCr3rDFRf 3r3) 

CC( Irl) = DER(l) + SAV 
H = l. 

UO 52 J=1 f 3 
H = h * FLOAT (J) 

52 CC( If J+l) = rER(2*J+l) / R 
CALL CF?Fl(Of ALCrOf A3»ALDt3rM) 

C ***COMPUT ING THE DISCREPANCIES FOR ThF PIFCE-WISt FITS. 

DO 53 1 = 1 rN 
0=1— ( MH-IS ) 

IF ( J-l ) 57r 57r 56 
56 IF ( J-MM ) 5fl r 58 r 57 

56 CALL CF2Fl(nrALCrOrA3rALO(2*J-l) r3rM) 

57 CALL CF2F2(TC(I) r ALC r 3 r DER r D r 3 ) 

53 DP ( I ) = SC(I) - DFR(l) 

E =T ( 1 ) 
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IZ=MH-IS 

H=T(XZ+1) 

DO 166 I = 1 rMM 
IF ( IPRNT ) 51 » 55» 51 
51 WRITE (6» 54) Rf H» (CC(I»J)» J=l»4) 

54 FORMAT (F7«2»3H T0Ffl.2»4E17.fl) 

55 6=H 
J=I+MH+1-IS 

IF ( 1 + 1— MM) 164 * l65r 164 

165 J=N 

H = T(J) 

GO TO 166 
164 H=T( J) 

166 CONTINUE 

WRITE (6» 124) MK 

124 FORMAT ( 17H PIECEWISE CUFIICS / 21H NO. OF POINTS LINKEr ,1?) 
WRITE 1 6 • 127 ) 

127 FORMAT (lH0r3(30Xr6HDELTAS) ) 

DO 151 J=1»N*1 

151 WRITE (6» 153) T(J) »S(J) rDP(J) 

153 FORMAT (F10 .4 »2X »F10 . 3» F15.6 ) 

CALL STnEV (DP»N»S2) 

WRITE <6» 11) S2 

11 FORMAT <19H0 STD. DEV. Fli.fa) 

RETURN 

END 


•IP FOR STATS 

SUBROUTINE STATS (Xr Nf XM» SD) 
DIMENSION X(100) 

SUM1 = 0. 

SUM2 = 0. 

XNO = N 

XNOl = XNO - 1. 

DENOM = XN0*XN01 
DO 10 I = 1*N 
SUM1 = SUM1 + X ( I ) 

10 SUM2 = SUM2 + X ( I ) *X ( I ) 

XM = SUM1 / XNO 

XNUM □ <XNO*SUM2 - SUM1*SUM1) 

IF (XNUM .LF. 0.) GO TO 14 
SD = SORT ( XNUM/nENOM ) 

IF ((SD/XM) .LT. .00025) GO TO 14 
GO TO 15 

14 SD = 0. 

15 RETURN 
END 
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GLS00001 


*IP FOR GLS1- 

SUBROUT INE GLS1 ( B * X * IL *N * M * ALPHA * El * E2 ) 

DIMENSION A (46*4 ) * B(46*4)f X(4)* IL(4) 

MM=M+1 GLS00003 


2 

DO 2 J = 1 f N 
DO 2 K = 1*MM 
A(JrK) = H(JtK) 



LL=1 

GLSO0004 


DO 60J=lfMM 

glsoooos 

60 

IL(J)=0 

GLS00006 


1=1 

6LSQ0007 


DO 3K=1 * MM 

GUSOOOOA 


11=1+1 

GLS00009 


DO 4 J=I I f N 

GLsoomo 


IF (ARS(A( J»K) )-El)4r4f6 

GLS00011 

6 

Tl=SQRT ( ( A ( J#K ) ) **2+ ( A ( I »K ) ) **2 ) 

GLS00012 


b=A ( J r K ) /T 1 

GLS00013 


C=A(IfK)/Tl 

GLS00014 


DO 5L=KtMM 

GLSQ0015 


T2=C*A(IrL)+S*A(J#L) 

GLS00016 


A ( J » L ) =— S* A ( I » L ) +C*A ( J r L ) 

GLSQ0017 

5 

A ( I rL ) =T2 

GLSQ001P 


LL=LL+1 

GLSQ0019 

4 

CON f INI IF. 

GLS00020 


IF ( ARS (A(IfK) ) — E2 ) 3 # 3 » 8 

GLS^0021 

a 

1 L < K ) = I 

GLS00022 


1=1 + 1 

GLSQ0023 

3 

CONT INI IF 

GLS00024 


X(MM)=-1.D 

GLS00025 


I T=M 

GLS00026 


DO 3bI=l*M 

GLS00027 

3b 

xm=n. 

GLSOOn28 


DO 30 J=1 f M 

GLS00029 


IF (IL(I1)) 30 f 30 » 31 

GLSoonSn 

31 

S = 0 • 

GLSQ0031 


LL= I 1 + 1 

GLS00032 


I=IL( II) 

GLS00033 


DO 32K=LLrMM 

6LS00034 

32 

S=S+A(IrK)+X(K) 

GLSQ0035 


X(II)=-S/A(IrII) 

GLS00036 

30 

11=11-1 

GLS00037 


IF (IL(MM) ) 50 *51 *50 

GLS0003F 

51 

ALPHA=0 • 

GLSO0039 


00 TO 52 

GLSO0D40 

50 

I = I L ( MM ) 

GLSO0P41 


ALPHA= A ( I *MM) 

GLSO0O42 

52 

RETURN 

GLS00043 


END 

GLSO0044 
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•IP FOK ACCFND 


4 


5 

b 

1 

H 


7 

2 


9 

3 

10 


SUBROUTINE ACCENO ( X r Y f N ) ASCN001 n 

THIS SR SORTS (XpY) POINTS INTO A SLGNENCE OF ASCtNOlNG X VALUES • ASCN0020 
N IS THE NO. OF POINTS IN THF SEOUFNCE WHILF X AND Y ARF ASS0CI«-ASCN0n30 
ATEU ARRAYS. THE ARRAYS OCCUPY THE SAMc. STORAGE aFTFR SORTING ASASCMQ040 
THEY HID HEFORF SORTING. N IS fHL NO. OF POINTS. ASCNOOSH 

DIMENSION XC1)#Y(1) ASCN0060 

DIMENSIONS OF ABOVE VARIABLES ARE ACTUALLY EFFFCTtD PY THE ASCMOO 7 O 

HIGHER (CALLING) PROGRAM OR SR. ASCNOOflO 

EOUI V/ 1 LENCE < I p T ) ASCM009P 

J=1 ASCN0100 

J IS THE INDEX OF THE NEXT MEMBER Of- THt SET OF POINTS WHICH WILLASCNOUP 
BE ORDERED BY OPERATIONS IN THE INNER LOOP HO 8 ON T • ASCM0120 

GO TO 3 ASCN013P 

THE ABOVF TRANSFER AVOIDS MIS-OPFKAT TON IF N=1 OR LE^Sp S£F CARD ASCM014O 
MO 305. NORMALLY p PROGRESS TO STMT. 4. ASCNOlbH 

K=J ASCN016P 

K IS THE TENTATIVE INDEX OF THF SMALLEST UN-ORDERtD X VALUE. ASCM0170 

I=J + 1 ASCNOIRO 

GO TO b ASCN019P 

1=1+1 ASCNOPOO 


I F ( X ( J)-X(K) ) 1 p 8 p 8 
K = I 

IF( I-N)Sp7p7 

K IS NO L0N6FR TENTATIVE. IT IS INDEED THF. INDEX OF SMALLFST Xp 
SO FAR UMOROERFD . 

I F ( K — J ) P p 9 r 2 
T=X(K) 

X ( K ) =X ( J ) 

X( J)=T 
T=Y ( K ) 

Y(K)=Y(J) 

Y(J)=T 

BOTH X AND Y HAVE BEEN SWAPPFDp USING T AS A TFMPORARY STORAGE. 
J=J+1 

IF ( J— N ) 4 p 10 p 10 
KFTURN 
a no 


ASCN0F1P 

ASCN0220 

ASCM0P30 

ASCN0P4P 

aSCMOFbP 

ASCN0P6P 

ASCM0P7P 

ASCM0R80 

ASCNOP90 

ASCN0300 

ASCM0310 

ASCN032P 

ASCM0330 

ASCN034P 

ASCNG35P 

ASCM036P 

ASCN0370 


IP FOR NOKMLZ 

SUHRul IT INF MORMLZ ( X p Y » N p AV 

THIS ROUTINE AVERAGES M ELEMENTS IN X ARRAY GETTING AV AS THE 
RESULT. IT THEM PRODUCES ARRAY Y WHICH IS SAMF AS XpMTNUS AV, 
TERM BY TERM 
DIMENSION X(l) pY(1 ) 

AV=X ( 1 ) 

DO 1 T=?,N 

1 AV=AV+XU) 

AV = AV / FLOAT (N) 

DO 2 1=1 pN 

2 Y(X)=X(I)-AV 
RETURN 

END 
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•IP FOR FIXIT 

SUBROUTINE F 1 X I T ( A , n > , N , N2 * N7 ) 

DIMENSION A (10,1) »B(5» 1) #0(50) ,KK(4> #S1 (IS) ,S?(15) 
00 200 J=N?,N 7 
00 100 1 = 1, N 
100 D(I)=A(JtI) 

200 CALL STDEV(n,N,Sl(J) ) 

KK(1)=B 
KK(2)=10 
KK (3)=12 
KK(4)=14 
00 400 J=l,4 
DO 300 1=1, N 
300 U(I>=P(J,I) 

400 CALL ST DEV ( D , N , S2 ( J ) ) 

WRITE <6, 11) (Sl(J)r J=N2,N7>, (S2(J), J=l»4) 

11 FORMAT ( 13H0 STD. DEV • , 6F 11.6, 2X,4F 11.6) 

RETURN 

END 


»IP FOR STuFV 

SUBROUTINE ST DEV ( X , N r S ) 
DIMENSION X(l) 

SUMXrO . 

XSQR=0. 

F.N=N 

DO 1U0 1=1, N 
SUMX=SUMX4*X(I) 

100 XSQR=XSOR+(X(I) )**2 

SP=XSOR-( (SUMX**2)/FN) 

S = SORT (SP / (EM - 1.) ) 

RETURN 

FNU 
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•IP FOR CF2F1 

SUHROUTINE CF2F1 ( J, ALC r IT» ALW » ALD»K »N) CF2 

1/ 1 MENS I ON TMP (?) r ALW ( 1 ) * ALP ( 1 ) » ALC ( 1 ) CF2 

K6=6*K+4 CF2 

FK=K CF2 

fjT = 2 CF2 

IF < IT) 1 1 * 1 » 11 CF2 

11 NT = 3 CF2 

1 IF(J>33,3»33 CF2 

33 TMP(2)=0.0 CF2 

TMP(l)=0.n CF2 

GO TO 2 CF2 

3 ALC(K6+4)=0.0 CF2 

ALC(Kf>+3)=0.0 CF2 

FN=N CF2 

NN=0 CF2 

DO r> 1 = 1 >N CF2 

WN = NN + NT CF2 

ALC ( Kf>+4 ) = ALC ( K6+4 ) + ALD ( NN ) CF2 

5 ALC ( KfS+3 ) =ALC ( K6+3 ) +ALO ( NN— 1 ) CF2 

TMP(l)=ALC(K6+3)/FN CF2 

ALC ( K6+3 ) =TMP ( 1 ) CF2 

TMP (2 ) =ALC (K6+4 ) /FN CF2 

ALC ( KfS+4 ) =TMP ( 2 ) CF2 

2 NN=0 CF2 

NM=-2 CF2 

DO 6 1 = 1 »N CF2 

NM=NM+5 CF2 

NN = NN + NT CF2 

ALW (NM ) =ALO ( NN) -TMP (2 ) CF2 

6 aLW(NM-1)=ALD(NN-1)-TMP(1) CF2 

NM=-2 CF2 

NN=— 4 - CF2 

DO 8 1=1 fM CF2 

NN=NN+5 CF2 

IF ( IT ) 77 f 7 r 77 CF2 

77 NM = NM + NT CF2 

ALW(NN)=ALO(NM) CF2 

GO TO 8 CF2 

7 ALW (NN) =1.0 CF2 

8 CONTINUE CF2 

NN=0 CF2 

DO 9 I = 1 » N CF2 

NNrNN+5 CF2 

9 ALW ( NN ) =1 • 0 CF2 

FJJ=0.0 CF2 

ALC (K6+5 ) =J CF2 

IF(J)17rl5fl7 CF2 

17 t- JJ=FJJ+1 .0 CF2 

NN=0 CF2 

DO 18 1=1 »N CF2 

NN=NN+5 CF2 

18 ALW(NN)=ALW<NN)*aLW(NN-2) CF2 

IF (F JJ-ALC (K6+5) )17»15»15 

15 ALC (K6+1 )=0.0 CF2 

ALC (KA+2 ) =0 . 0 CF2 

ALC ( Kft— 2 ) =0 • 0 CF2 

NN=— 4 CF2 

LO 19 1=1 rN CF2 

NNrNN+5 CF2 

I MP ( 1 ) = ALW ( NN ) * ALW ( NN+4 ) CF2 



ALC (K*+?)=TMP(1 ) *ALW ( NN+1 ) +ALC ( Kb+2 ) 

7 MP (1 ) =TMP < 1 ) * ALW ( NN+4 ) 

ALC (K6+1 ) =TMP ( 1 ) + ALC (KF + 1 ) 

19 ALC(K6-2)=TMP(l)*ALV/(NN+2)+ALC(K6~2) 

ALC ( KF-2 > = *LC ( Kb-2 ) / ALC ( K^+l ) 

11 1 ALC ( K0 + 2 ) = ALC ( K6+2 ) / ALC ( Kf*+1 ) 

112 ALC(K6-3)=0.0 
24? L3=D 

Ll=-1 

L4=0 

113 K6=Kf>-6 

114 lF(FJJ-FK)U5rl?2f 122 
lib FJJ=FJJ+1.0 

mLC <Kb + 2)=fU0 
ALC(Kb+l>=0.0 
ALC (Kb-?)=0.0 
ImM=-4 

00 119 1=1 fN 

i>n=nn+s 

L2=NN+L1 

LS=HU+L4 

1HP(1 )=ALC (K6+3) *ALW(L?+4) 

^Lw ( L? +4 ) = ( ALW ( NW+2 ) - M.C ( K b+4 ) ) * aLw ( L5+4 ) -TMP ( 1 ) 
TMP(l) = ALW ( L?+4 ) * ALW CNN) 

ALC ( Kb + 2 ) = T MP ( 1 ) * ALW (NN+1) + ALC ( Kb+2 > 
TMP(1)=TMP(1 ) * A L/H L2+4 ) 

ALC ( Kb + 1 )rTMP(l )+ALC(Kb+l) 

119 mLC ( K b— 2 ) = TMP(l) * ALWCMM+3) + ALC(Kb-2> 
/>LC(Kb-?)=ALC(Kb-?)/ALC (Kb+t ) 

ALC ( K b+? ) =ALC ( Kb ♦*? ) /mLC ( Kb+1 ) 

/.LC (Kb-3)=ALC <l<b+l)/ALC (Kb+7) 

IF ( L3 ) 94? »l?0r?42 

1?0 L 3 = 1 
L 1 = 0 
L4=-l 
00 ro 113 
1?2 kf TUKN 
1ND 


CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 

CF2 


115 



* IP FOR CF2F2 



SUBROUTINE CF2F2 ( ARG * ALC fKRfOERr IQ#K ) 

CF2 


DIMENSION TMP(6) f ALC(l) ^DER(l) 

CF2 


XRAR=ARG 

CF2 


FKB=KR 

CF2 


K6=6*K+4 

cFa 


IF(ALC<K6+5) >22,222*22 

CF2 

22? 

XB AR=XBAR-ALC ( K6+4 > 

CF2 

22 

F JJ=ALC ( K6+5 ) 

CF2 


MN=K6+6 

CF2 


NQ=IG+1 

CF2 


DO 23 I=1*NG 

CF2 


NNsNN-6 

CF2 

23 

ALC(NN)=0.0 

CF2 


TMP(l)=n.n 

CF2 

24 

IMP (2) =1-0 

CF2 

25 

IF ( ALC (K6+5 ) -TMP ( 1 ) ) 231 1 r 231 ,2311 

CF2 

2311 

TMP ( 1 ) =TMP ( 1 > +1 . 0 

CF2 


TMP(2)sTMP(2)*XRAR 

CF2 


1F<TMP(2) ) 25,24,25 

CF2 

231 

ALCtK6>=TMP<2> 

CF2 


IF ( TMP ( 1 ) ) 2322 , 232 r 232? 

CF2 

2322 

KK6=K6 

CF2 

236 

KK6=KK6-6 

CF2 


IF ( XHAR ) 234,2344 , 234 

CF2 

2344 

TMP ( 2 ) =ALC ( KK6+6 ) 

CF2 


AUC(KK6+6)=0*0 

CF2 


GO TO 235 

CF2 

234 

TMP ( 2 ) =ALC ( KK6+6 ) /XBAR 

CF2 

235 

ALC ( KK6 ) =TMP-( 2 ) *TMP ( 1 ) 

CF2 


TMP(l)=TMPm~1.0 

CF2 


IF (TMP ( 1 ) >236,232,236 

CF2 

232 

NN=-1 

CF2 


TMP(l)=ALC(Kr>+2) 

CF2 


TMP(2)=ALC(K6+1) 

CF2 


KK6S=K6 

CF2 


KK6SK6+A 

CF2 


DO 241 1=1, NO 

CF2 


KK6=KK6-6 

CF2 


NNsNN+2 

CF2 


DER ( NN ) =TMP ( 1 ) * ALC ( KK6 ) 

CF2 

241 

DEH ( NN+1 ) = ( ALC < KK6 ) **2 ) /TMP ( 2 ) 

CF2 

242 

L3=n 

CF2 


L7=0 

CF2 


Ll=-1 

CF2 


L4=6 

CF2 

26 

IF (FJJ-FKB ) 211 1* 211*211 

CF2 

2111 

f JJsFJJ+1.0 

CF2 


KKASKK6S 

CF2 


NN=-1 

CF2 


TMP(l)=ALC(KK6-2) 

CF2 


TMP (2 ) ? ALC ( KK6-3 ) 

CF 2 


TMP ( 3 ) =ALC ( KK6-4 ) 

CF2 


TMP (4 ) =ALC ( KK6— 5 ) 

CF2 


KK6S=KK6-6 

CF2 


TMP (6) =0 « 0 

CF2 


KK6=K6+6 

CF2 


DO 29 1 = 1 , NO 

CF2 


KK&=KK6-6 

CF2 


NN=NN+2 

CF2 


L6=KK6+L7 

CF2 
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L2=KK6+L1 CF2 

TMP(5)=TMP(2)*ALC(L2> CP2 

TMP(5>=<XBAR-TMP<1> ) *>ALC <L6) -TMP (R) CF2 

L5=KKfc+L4 CF2 

ALC ( L2 ) =TMP < 6 > * ALC < L5 > +TMP ( 5 ) CF2 

UER (NN)=AUC <L2> *TMP (3J+DER INN) CF2 

DER(NN+1)=< (ALC(L2)**2)/TMP<4> )+UER(NN+l) LF2 

29 TMP (6) 3TMP (6) +1 »0 CF2 

IF(U3)242»30»242 CF2 

30 L3=l CF2 

L1=0 CF2 

L4=5 

L7=-l CF2 

00 TO 26 ^ F 2 

211 IF (ALC (K6+5) )212»2122»212 CF2 

2122 DER ( 1 ) 2DER ( 1 ) +ALC ( K6+3 ) CF2 

212 RETURN CF2 

tND 

• XOT MAIN 


NASA- Langley, 1966 s-105 
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" The aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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